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ABSTRACT 

In this paper concepts and numerical calculations are reviewed in which baryons 
are considered as many body systems composed of three quarks coupled to a polar- 
ized Dirac sea. The interaction is the one of the Nambu-Jona-Lasinio model and 
of local and biquadratic form in a way which allows for spontaneous chiral symme- 
try breaking. The models are treated explicitely in the 0-boson and 1-quark loop 
approximation (Hartree mean field approach), where the real part of the effective 
Euclidean action is gauge invariantly regularized. In the meson sector the parame- 
ters of the models, including a cut-off for the Dirac sea, are fixed to PCAC, on-shell 
meson masses and decay constants. Furthermore the regularization schemes are 
chosen to get good values for vacuum condensates and current masses leaving the 
constituent mass M as only free parameter. The review is focussed on the evalua- 
tion of baryon properties. Assuming a spherically symmetric hedgehog ansatz for 
the quark fields the solitonic solutions of the system are obtained in a selfconsistent 
way by solving the classical equations of motion. 

After a semiclassical quantization procedure the solitons can be related to the 
physical states of the spin 1/2 and 3/2 multiplets. In this way observables and form 
factors can be calculated. Various types of quark fields and quark-quark couplings 
are investigated and reviewed. These are in SU(2) given by sigma and pion fields 
and in SU(3) they are complemented by kaon and eta-fields. In addition SU(2) 
vector mesons of p, Ai and cu-type are considered as well. For these models relevant 
observables and form factors of the nucleons are calculated. Some results on the 
other members of the baryon octet and decuplet are reviewed as far as they are 
available and of theoretical and experimental interest. For scalar and pseudoscalar 
couplings a clear picture emerges corresponding to localized valence quarks coupled 
to moderately polarized sea quarks and being separated from them by a finite 
energy gap. Using heat kernel and gradient expansion techniques the models of 
the Nambu-Jona-Lasinio type can be related to fully bosonized approaches of the 
Skyrme type, on the one hand, and to the chiral sigma model, representing valence 
quarks and dynamical meson fields, on the other. 



1. Introduction 

The present article deals witii tiie description of nucleons and hyperons by 
means of relativistic chiral effective models. The prominent model discussed in 
detail is the Nambu-Jona-Lasinio (NJL) approach. It incorporates certain features 
of the Quantum Chromodynamics (QCD) and is related to low energy strong 
interaction phenomena. 

Quantum Chromodynamics is generally considered to be the proper theoretical 
framework for the description of structures and reactions being dominated by the 
strong interaction. QCD is a SU(3) colour gauge theory combining quarks of 
Nf different flavours (u,d,s,c,b,t) as fermions and gluons (g) as gauge fields. Its 
Lagrangian reads 

with the covariant derivative 



D^ = di, + igAl-\^ 



1 
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and the field strength tensor 



Gl, = d^Al - d,A- - gr^'^AlAl 

Here A" are the generators of SU{?>)c (Gell-Mann matrices) the /q^^ an the struc- 
ture coefficients of SU{Z) and rrij are the current masses of the quarks. The basic 
features of QCD are (Donoghue et al. 1992) the following being all related to 
symmetry properties: 

Universality: There is only one coupling constant g for all interactions among 
quarks and gluons. This is a direct consequence of the local gauge invariance under 
the 5't/(A^c)-gi'oup, with Nc = number of colours. 

Asymptotic freedom: At very high energies quarks behave as free particles. In 
this region the effective coupling constant becomes small and perturbation theory 
may be applied. This property has been observed experimentally. 

Confinement : Up to now neither free quarks nor gluons have been observed 
outside the volume occupied by a hadron. This feature has not fully been explained 
yet by theoretical reasoning. The idea is that only colourless systems are stable. 



Chiral Symmetry: In the limit of vanishing quark masses, the QCD La- 
grangian is invariant under the chiral SU{Nf)ji (g) SU{Nf)i group of global 
transformations 

q{x) ^ exp (za'*A"/2)Q(x) 



g(x)^exp(z75/?«A72)g(x) 
where a — 1 , . . . , Ni — 1 , so that the corresponding conserved currents read 



Aiix) 
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In the absence of strong interactions the quarks have a current mass of electroweak 
origin which breaks chiral symmetry explicitly. That means 

d^'A'l;^ ^ iq{x){m, X^'jqix) 

where the current quark mass matrix m is defined as ?fi = diag{mu, mj^, rus, rric, m^, mt) 
The attributed values are (Gasser and Leutwyler 1982) given in Tab. 1.1. 

The chiral symmetry is well realized with respect to up- and down-quarks 
because of their small current mass (^lOMeV). The symmetry is spontaneously 
broken resulting in a chiral condensate < uu + dd >^/'^= — (283±31)MeV (Gasser 
and Leutwyler 1992). This gives the up- and down-quarks a dynamically generated 
constituent mass of about 350 — 450MeV. For those quarks apparently the spon- 
taneous breaking of chiral symmetry dominates over the explicit breaking caused 
by finite current masses. For heavy quarks (c,b,t) the situation is rather opposite. 
The strange quark s lies in an intermediate region and its relevance for low-energy 
non-strange physics is at present under debate. 

Finally let us mention that in the presence of external fields chiral symmetry is 
also broken because of quantum corrections leading to the chiral anomaly (Adler 
1969, Bell and Jackiw 1969). Vector current conservation implies an anomaly in 
the divergence of the axial current (Bardeen 1069), which predicts successfully the 
decay rate for the process tv'^ -^ 27 if Nc = 3. 



Scale Invar iance : In the chiral limit {nij = for all flavours) the QCD La- 
grangian does not contain any dimensional parameters. Thus it is invariant under 
the scale transformation q{x) -^ X^ ' q{Xx). The corresponding dilatational cur- 
rent is not conserved because of quantum corrections leading to the trace anomaly. 

The QCD has been treated in a clean way only for high energy processes, be- 
cause there the asymptotic freedom (Politzer 1974, Gross und Wilczek 1973a, 1973b) 
facilitates the solution tremendously (perturbative regime). For low energies, in 
particular the structure of hadrons, the gluon self interaction ~ {gf Aj^^A^) 
creates non-linearities and causes enormous problems (non-perturbative regime). 
Here the Lagrangian as such has only been treated in a regularized way on a 4- 
dimensional lattice in the form of lattice gauge theory (Wegner 1971, Wilson 1974). 
Although these calculations are potentially exact they are practically hampered 
by technical and conceptual problems associated with the choice of large current 
masses, fermion doubling, small size of lattices, etc. In addition calculations of 
this sort require extreme amounts of computer time (Creutz 1983). Because of 
this dilemma in the non-perturbative regime effective theories are very much in 
use. They explicitely treat those degrees of freedom, which are relevant for low 
energy structures, and ignore the others or parametrize them in form of coupling 
constants. This only makes sense if those low energy degrees of freedom are clearly 
separated from the high energy ones. However, such a question can be decided 
only by an inspection of the low energy phenomenology, which to some extend has 
been done already long before the advent of QCD. Of course for such a limited 
theory one must not expect that all properties of QCD, discussed above, are also 
properties of the effective models. Hence the guideline for the effective models is 
the necessity to reproduce low energy phenomena and to reflect some of the QCD 
symmetry properties. 

There are some basic phenomena of low energy hadronic physics which over the 
many years of experience with effective models have been identifled as indispensible 
for the construction of an effective theory for baryon and meson structure (Cheng 
and Li 1984). It is clear that the pion fleld plays a dominant role in all hadron 
physics. This goes back to Yukawa's original idea that nuclear forces are mediated 
by pions and corresponds to the fact that basic feature of nuclear structure at 
low and intermediate energies can successfully be explained in turn of nucleons 



and mesons instead of their constituents quarks and gluons, as they are identified 
at high energies. The pion is closely connected to the chiral condensate. In fact 
it is the most prominent Goldstone-boson (Goldstone 1961, Nambu 1960) of the 
spontaneously broken chiral symmetry. 

The pion decays via the axial current A"(a;) and its decay constant /jr = 
93MeV is defined as (Cheng and Li 1984) 

< I A«(0) I n\p) >= tUp^5^^ 

For hadrons the axial current A^(a;) is partially conserved due to the small but 
finite current masses of the up- and down-quarks. According to the hypothesis of 
the partial conservation of the axial current (PC AC) (Nambu 1960, Chou 1961, 
Gell-Mann and Levy 1960) 

d^'Alix) = Umln^ix) 

with TOtt = 139MeV being the mass of the pion field t:"'. The iso-vector current 
Vn{x) is conserved unless one explicitely attributes different current masses to up- 
and down-quarks. Both currents allow the formulation of chiral charges (Cheng 
and Li 1984) 

3, 



Q^x) = J d-'xVQ^ix) 

Q\{x)= jd^'xAlix) 

and a corresponding charge algebra, which exists even if the currents are not fully 
conserved: 

[Ql,{x),Q\,{x)]=ir'^Q'y{x) 

m{^),Q\{x)] = ir^''Q\{x) 

[Q\{x),Q\{x)]=tr^^Q'y{x) 

The relevance of this commutator algebra lies in the fact that several low energy 
theorems and sum rules involving strong and electro- weak processes can be derived 
and confronted with experimental data. A typical example is the Goldberger- 
Treiman relation (Goldberger and Treiman 1958), which is experimentally fulfilled 
to a precision of 7%: 



Here M^y = 938MeV represents the nucleon mass, gj\ = 1.25 is the axial couphng 
constant of the neutron beta-decay and Qt^nn = 13.6 is the pion-nucleon-nucleon 
couphng constant. In general the violation of the chiral symmetry is less than 10% 
and it is indeed after the isospin symmetry the best known symmetry of strong 
interactions at low energies. It is clear that any effective theory of hadron has 
to incorporate the spontaneously broken chiral symmetry and the corresponding 
Goldstone boson, i.e. pions in SU{2) and in addition the kaons and the eta in 
SU{3). 

Another fruitful idea in the study of electromagnetic properties of hadrons is 
given by the vector meson dominance model, which assumes that photons do not 
interact directly with hadrons but rather through the virtual propagation of an 
intermediate massive vector meson state (Sakurai 1960,1969). In quantum field 
theory such an idea can be formulated through the so-called current field identities 
(KroU et al. 1967) 

9 2 



4^' - ^^-,(.), jf 1 = ^^p«(.) 



9ujNN ' ■ "" 9pNN 

for the isoscalar and isovector components of the electromagnetic currents. Here 
the vector meson masses rriuj = 783MeV and nip = 770MeV appear together with 
the p and uj exchange coupling constants of NN scattering. In most cases, the pre- 
dictions of the vector meson dominance model (VDM) have been very satisfactory 
(see e.g. Gourdin 1974). The most famous example is provided by the prediction 
for the mean squared electromagnetic pion radius 

< r2 >^"^ = -^ = 0.39fm2 (1.1) 

mi 

which compares reasonably well with the corresponding experimental value of 
(0.44 ± 0.01)fm2. 

Finally, let us mention that the union of current algebra with current-field 
identities leads to the field algebra whose most prominent result is the KSFR 
relation (Kawarabayashi and Suzuki 1966, Riazuddin and Fayazuddin 1967) 

with (/pTTTT = 5.48 being the decay constant of the strong process p ^ nn. The 
accuracy of this relation is 2%. 



As soon as strange quarks are involved one has the Gell-Mann Okubo mass 
formula (Gell-Mann, 1962; Okubo 1962), which deviates from the experiment by 
less than 3%: 

M^-Mn = ^(Ms - Mn) + ^(Ms - Ma) 



254MeV 190McV 58MeV 

and the Coleman- Glashow formula (Coleman and Glashow, 1961) 

Mn- Mp + M^- - M^o = M^- - M^+ 
1.3McV 6.4±0.6MeV 8±0.8MeV 

which fits perfectly with experimental data and which was originally attributed to 
the electromagnetic interactions between the quarks. 

An effective theory of low energy QCD degrees of freedom in the baryonic 
sector can only be successful if it reproduces the above phenomenological features. 
Therefore we concentrate in this article on models showing spontaneously broken 
SU{Nf)]^ (g) SU{Nf)ji symmetry with a possible explicit breaking due to finite 
current masses. In some cases we will also consider scale invariance. We consider in 
detail the Nambu-Jona-Lasinio model (Nambu and Jona-Lasinio 1961a, 1961b) and 
will relate it to all others. The theory will be explained for scalar and pseudoscalar 
quark-quark interactions in the light quark sector (up, down). Generalisations 
towards vector interactions and/or SU(3) will be done in separate chapters. Thus 
the NJL-Lagrangian studied in detail is 

G 
C-NJL = qih^d^ - mo)q + — [{qqf + [qi-f^rqf 

with q having an up- and down-component and ttiq = diag{mm 'n^d)- In rnost cases 
we will assume ruu = m^. The NJL-model is distinguished because it allowes to 
build conceptual bridges to most of the currently used effectives models for baryon 
structure. In such a way these models and also the NJL model can be better 
understood. 

There are by now several attempts to relate the NJL to some low energy 
limit of QCD. One has been given by Diakonov and Petrov (1986) assuming the 
QCD vacuum to be an instanton liquid. By variational field theoretic means they 
obtain after various plausible approximations the zero-boson and one-quark-loop 



approximation to the NJL-Lagrangian. This procedure is based on previous works 
by t'Hooft (1976a,1976b), Callan, Dashen and Gross (1978,1979a,1979b), Carhtz 
and Cremer (1979) and Shifman, Vainstein and Zakharov (1980a,1980b), see also 
Shuryak (1986). BaU (1987) on the other hand integrates out the gluon fields 
in some certain approximation and obtaines in the end the NJL-Lagrangian in 
SU{Nf)ji (g) SU{N f)i including vector mesonic couplings. Schaden et al. (1990) 
formulate QCD in terms of path integrals over the field tensor rather than the 
vector potential. At a certain level they perform stationary phase approximations 
and end up, after several further plausible arguments, at the NJL-model. Cahill, 
Roberts and others (Cahill and Roberts, 1985; Roberts et al., 1988; Cahill, 1992) 
introduce bilocal meson and diquark fields into a functional form of QCD. The 
NJL model appears as some pointlike approximation. There are also steps in a 
similar direction by Dhar and Wadia (1984, 1985), McKay and Munczek (1985), 
Adrianov (1985) and Chanfray et al. (1991). One should note that none of these 
theories is really a clean cut derivation of NJL. For example most of them are not 
able to say where in the chain of arguments and approximations on the way from 
QCD to NJL the confinement has been lost, nor can they say how the NJL-model 
can be improved. 

Actually the effective chiral models for baryon structure can roughly be clas- 
sified in three groups: 

i) Chiral quark models: Here the system is dominated by three valence quarks 
coupled to a dynamical pion. In die chiral MIT- model (Chodos and Thorn, 
1975) the quarks are confined by a bag with infinite walls and chiral sym- 
metry is preserved by a dynamical pion field and the continuity of the axial 
current at the bag's surface. The linear versions of these models can be re- 
lated to the cloudy bag model (Theberge et al. 1980; Thomas AW 1983). 
In the chiral a-model (Gell-Mann and Levy 1960) we do not have confining 
walls but a sigma field with finite vacuum value. Here the valence quarks are 
localized by binding forces exerted by the sigma- and pion field (Birse and 
Banerjee 1984,1985, for a review see e.g. Birse 1990 and references therein). 
The baryon number is carried by the three valence quarks. Altogether these 
models can be characterized by valence quarks and meson fields. 



ii) Skyrme models: In these models (Skyrme 1961, Adkins et al 1983, Zahed 
and Brown 1986, Holzwarth and Schwesinger 1987) one does not deal with 
explicit quarks but solely with dynamical meson fields. They are collectively 
quantized in order to get the proper quantum numbers of e.g. a nucleon. 
The formalism is comparatively simple and hence these models enjoy great 
popularity. Conceptually they are based on the large A^c-expansion of the 
QCD ('t Hooft 1973, Witten 1979b). This would require an infinite set of 
boson fields which is in practice replaced by the pion field in the simplest 
form, and by vector mesonic fields in more recent versions (Meissner U G 
1988). The baryon number is hidden in the topology of the Goldstone fields. 
At the center of the baryon the Goldstone fields must show a winding number, 
which is to be identified with the baryon number. These models can be 
characterized by meson fields and topology. 

iii) Nambu-Jona-Lasinio models: As described above one deals solely with quarks 
without dynamic meson fields (Nambu and Jona-Lasinio 1961a, 1961b). Va- 
lence quarks appear in a natural way as bound, discrete and localized states 
(Kahana and Ripka 1984). Mesons enter the theory only as non-dynamic 
auxiliary fields generated by the polarization of the Dirac sea (Eguchi 1975, 
Kleinert 1976). In fact these models can be characterized by valence quarks 
and sea quarks. 

iv) With respect to the above characterization there are some hybrid models 
which combine certain features. This is the chiral bag of the Stony Brook 
group (Brown and Rho 1979). There the pion field outside the bag shows 
some topological properties and carries hence a fractional baryon number. 
The linear chiral sigma model of Kahana and Ripka (1984) involves sea 
quarks and dynamic pion- and sigma fields and requires renormalization: In 
fact these authors are the first to treat the polarization of the Dirac sea 
in this context (Ripka and Kahana 1987, Soni 1987, Li et al. 1988). It 
shows, however, an instability and hence has never been applied to nucleon 
properties. Recently it was shown by Kahana and Ripka (1992) that the 
introduction of vector mesons in this model can cure the vacuum instability, 
though no numerical calculations for nucleons were performed yet. 
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It is interesting to note, that the NJL-Model (valence- and sea-quarks) can 
be directly related to the chiral bag models (valence-quarks and mesons) and the 
Skyrme- models (topology and mesons). In order to understand this, the figures 
1.1-1.4 may be helpful. Fig. 1.1 shows schematically the single particles energies 
of the chirally symmetric plane wave vacuum as it emerges from the NJL-model in 
the 0-boson and 1-quark loop approximation. We concentrate on up- and down- 
quarks and hence a small gap of 2mo = niu+m^ exists between the upper and lower 
continuum. In a plausible Fock-state picture the lower continuum is occupied by 
Nc = 3 quarks for each fiavour in each single particle level, as indicated by crosses. 
The spontaneously broken chiral vacuum is shown in Fig. 1.2. There the gap is 2M 
and the single particle states are again of plane wave nature with the constituent 
mass M. The vacuum in both cases is characterized by the baryon number B = 0. 
The soliton for a 'small' constituent mass M is shown in Fig. 1.3. A single particle 
state from the positive continuum has entered the gap and has become a bound 
state. It is explicitely occupied by A^c = 3 valence quark and hence the baryon 
number of the system has increased from B = to B = 1. The interaction of 
the continuum state with the valence quarks causes a polarization of the single 
particle states of the positive and negative continuum which is indicated by wavy 
lines. For large constituent masses the bound state has moved further down, see 
Fig. 1.4, and approaches the Dirac sea. In such a case both continua are strongly 
polarized and the Dirac sea incorporates another level and hence acquires the 
baryon number B = 1. The scenarios of Fig. 1.3 and Fig. 1.4 are distinguished by 
the values of the constituent mass M. As we will see for sigma and pionic couplings 
in SU(2) as well as in SU(3) the scenario of Fig. 1.3 is the one which reproduces 
the relevant nucleon and hyperon observables. If vector mesons are included the 
picture is not yet settled. 

The formal techniques to relate the NJL model to the bag model and to the 
Skyrme model are the gradient expansion (Aitchison and Frazer 1984, 1985a and 
1985b) the heat kernel expansion (Gasser and Leutwyler 1984, Ebert und Rein- 
hardt 1986), or the expansion of Chan (1985) all applied to the fermion determi- 
nant of NJL. Actually these expansions are well converging for smoothly varying 
and largely extended boson fields or quark wave functions which is not always 
fulfilled for realistic systems. Nevertheless, if one uses only a certain number of 
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expansion coefficients in the expansion of tfie effective action of tfie scenario of 
Fig. 1.3, the Lagrangian of GeU-Mann and Levi (1960) is obtained. There the 
baryon number is carried by the valence quarks and the polarization of the Dirac 
sea creates dynamic meson fields whose possible topology does not carry a baryon 
number. If one considers the scenario of Fig. 1.4 and terminates the expansion of 
the action in an appropriate way and if one ignores certain destabilizing terms then 
a Skyrme type Lagrangian (Skyrnie 1961) is obtained. The highly polarized Dirac 
sea including the additional valence level is described in terms of a dynamic pion 
field. This exhibits a winding number to be identified with the baryon number. 
One has to consider these relations between the models with some care: First, 
the expansions are often badly converging. Second, a gradient expansion of an 
observable of NJL does not necessarily agree with the observable evaluated from 
the gradient expanded NJL Lagrangian. Nevertheless the NJL model is extremely 
helpful in understanding the meaning of valence quarks and topology the concepts 
of which are the basic features of all presently used effective chiral models. 

There are two other branches relevant in the context of the NJL-model, one 
is the meson sector with baryon number B = 0, the other is the evaluation of 
properties of a hot and dense medium and the embedding of a baryon in it. The 
meson sector has been the object of intensive studies in the last years and several 
review articles have been written on it (Vogl and Weise 1991, Klevansky 1992). 
We will discuss the meson sector in this paper only as far it is needed to explain 
the appearance of the chiral condensate and to adjust the parameters of NJL. 
The medium calculations are not discussed in the present article at all. The 
interesting point is that with increasing density and temperature of the medium 
the chiral symmetry of the system is restored. The reader is referred to Hatsuda 
and Kunihiro (1985, 1987a, 1987b, 1988), Meissner U G (1989a, 1989b), Bernard 
V et al. (1986, 1987), Jaminon et al. (1989, 1992), Ruiz Arriola E et al. (1990), 
Christov C V et al. (1990a, 1990b), Lutz and Weise (1991) and Hatsuda and 
Kunihiro (1994). 
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2. The Nambu-Jona-Lasinio-Model With SU(2)-Flavor: 
Vacuum and Mesonic Sector 

The purpose of this chapter is to present the Nambu-Jona-Lasinio-model 
(NJL) (Nambu and Jona-Lasinio 1961) in the simplest realistic case, i.e. with 
scalar and pseudoscalar couplings. To begin with, we study the classical La- 
grangian and its symmetries. The bosonization procedure is presented by means 
of path integrals and the stationary phase approximation is defined. We extract 
the ultraviolet structure of the model and present some different regularization 
schemes. The gradient and heat kernel expansions are discussed as well as the 
formalism for extracting vacuum and mesonic properties. The parameters of the 
model are fixed in order to reproduce known properties of the mesonic sector. Dif- 
ferent ways of parameter fixings are introduced and the corresponding numerical 
results are shown. 

2.1. The Semibosonized NJL - Effective Action 

Classical Lagrangian - Symmetries 

The simplest version of the Nambu-Jona-Lasinio-model is represented by the 
following Lagrangian 

^NJL = ^kin + ^br + ^int (2-1) 

where the kinetic, interaction and chiral breaking mass terms are given by 

^km = q{x)i^q{x) 

^br = -q{x)mq{x) (2.2) 



r -G 

'-'int r) 



2 / . . ^ . .\2 



q{x)q{x)) + [q{x)i'^^Tq{x 



respectively. Here q{x) denotes a quark field with u and d fiavors and Nc col- 
ors. The G is the coupling constant with dimensions of length squared. The m 
represents the current quark mass matrix given by 

( ruu \ 
m= =mil + r3Am (2-3) 

V md) 

where the average and difference masses are defined as 

'^1 = -^{mu + md), Am = mu-md (2.4). 
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Their attributed theoretical values at a scale of IGeV are (Gasser and Leutwyler 
1982) given in Tab. 1.1. In the following we will assume the quark mass degeneracy 
ruu = m^ = mi, corresponding to exact isospin symmetry. 

At the classical level the total Lagrangian is invariant under the global U{1) 
transformation 

q{x) ^' q (x) = exp{id)q{x) (2.5) 

The kinetic and interacting terms Cj^j^j^ and £j^i are also invariant under the 
SU{2)v ® SU{2)a global chiral group 

q(x) -^ q'ix) = expiif ■ d)q(x) 

(2.6) 
q{x) -^ q'{x) = exp(i75f ■ (3)q{x) 

These transformation properties generate the following Noether currents (Cheng 
and Li 1984) 

B^{x) = q{x)'^^q{x)j (baryon current) 

-^ _ f 

V^{x) = q{x)'y^—q{x), (vector current) (2.7) 

-> _ 'f 

^m(^) = ^(^)7m752^(^)' (axial current) 

whose divergences on the classical level are given by: 

c>^S^ = 0, 8^"% = ^, ^^'A^='-q{m,r}q (2.8) 

corresponding to the conservation of the baryon number and isospin and the partial 
conservation of the axial charge respectively. Actually , even for m = the axial 
current is not conserved on the quantum level in the presence of external vector 
and axial vector fields refiecting the chiral anomaly (Wess and Zumino 1971, Wess 
1972). Finally let us mention that the total Lagrangian shows a SU{Nc) global 
invariance, leading to the conservation of the color current. Since the four fermion 
interaction in eq. (2.1) is of color singlet form, this conservation does not have 
dynamical implications. However, there are other versions of the model which 
involve color octet-octet interaction in the vector and axial channels (Klimt et 
al. 1990 ; Takizawa et al. 1990) which are formally equivalent to the Lagrangian 
(2.2) after a Fierz rearrangement. Those have been discussed at length in other 
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reviews ( Vogl and Weise 1991; Klevansky 1992) and will not be considered here 
specifically. 

Collective Boson Fields 

Like the Fermi theory of weak interaction the Nambu-Jona-Lasinio model is 
non renormalizable, because the coupling constant of the 4-fermion interaction G 
has the dimension [G] = mass (Itzykson and Zuber 1980, Cheng and Li 1984). 
This means that with each increasing order in G a new graph with a higher degree 
of ultraviolet divergence appears. In order to get a well defined theory it is therefore 
necessary to specify how the infinities of the model have to be treated. 

in which order the model has to be handled. 

For our later purposes it will be convenient to use a semibosonized form of 
the 4-fermion interaction eq. (2.1). The idea of bosonization has first been formu- 
lated in solid state physics, where it is called Hubbard-Stratononovich transfor- 
mation (Negele and Orland 1987) and has been applied to the present theory by 
Eguchi (1976), Kikkawa (1976) and Kleinert (1978). It consists in resummarizing 
the graphs of the 4-fermion point interaction into a quark-meson interaction of 
the Yukawa type by introducing collective scalar-isoscalar (a) and pseudoscalar- 
isovector (yf) background fields, which carry the quantum numbers of the interac- 
tion channels (qq) and (qir'^^q), respectively. To this end we insert the functional 
identity: 

1 = I VaVnexp!^-t I d^x^[{a + ^{qq - ^))2 + (^ + ^q^r ^^qf]^ (2.9) 
into the generating functional Zjsfjj;^: 

and obtain: 

Z'njl = f VqVqVaVne' ^ '^^'^^NJLi^) (2.11) 



with 



Lat jt = qiDq —[a + tt J H a 

iD = i^ — g{a + inr'-y^) 
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and: 

G = g^/^? (2.13) 

Here g and // are newly introduced parameters which arise due to the present form 
of bosonization. Their values will be fixed in sect. 2.6.. 

We want to stress that on this stage the fields a and tt are non-dynamical 
collective background fields and no kinetic term \ d^ad^^a + d^Tid^Tr appears on 
the classical level. 

The semibosonized Lagrangian C j^jj^ remains invariant under the chiral trans- 
formation of SU{2)ji (g) SU{2)i ® U{1)q, if (cr, yf) transforms under the (^, \) 
representation of this group (Cheng and Li 1984). Hence we demand: 

0- ^ cr + 2J3n 

^ ^ ^ ^ ^ (2.14) 

TT ^ TT — 2pa — 2a X Tx 

eq. (2.14) is consistent with the equations of motion for C' 

Grassmann-Integration - Effective Chiral Action 

The quark contribution qiDq in eq. (2.12) is bilinear in the quark fields q 
and q so that the functional Grassmann integration JVqVq can be performed 
analytically. After transformation to Euclidean space time (cf. appendix A) we 
get: 

Znjl" = jVaVne-^-ff^''^'^^ (2.15) 

The effective chiral action S^ff{a, tt) contains the fermion determinant S^jj (a, n) 
as well as an mesonic mass term Sgff^'^''^{a,n) and the chiral breaking term 

S,ff{a, Tf) = S,ffP{a, Tf) + 5e//"^''(^, ^) + Sef/'i^. Tf) 
Sef/{a,n) = -lndet(-iD) = -Spln(-iD) 

5,;;--(a,7f) = ^/d4^^(a2 + 7f2) 

^eff (c^,7r) = d XECr 

9 •^ 

Sp denotes the total trace in functional as well as matrix space: 

Sp^ = I d'^xTr^rci.x \A\ x)= I d'^kTr^rcik \A\ k) 



(2.16) 
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Generally the fermion determinant in Euclidean space time has a real as well 
as an imaginary part, reading: 

ReS,f/{a,n) = (-)^Spln(Dtl))| (2.17) 

and: 

lmS,f/{a,n) = (-)^Spln (-^) (2.18) 

For time independent meson fields a and if it will be shown in chapter 3, that the 
imaginary part ImS'gff (cr, tt) vanishes, because the 1- particle hamiltonian h is 
hermitian. This changes if time-like vector mesons (e.g. u mesons) are coupled 
to the Lagrangian destroying the hermiticity of h. In some cases the imaginary 
part is related to the chiral anomaly. Details will be discussed in chapter 7. In 
addition imaginary parts of the action occur if the system is considered in the 
rotating frame. 

2.2. Constituent Quark Mass - Stationary Phase Approximation 

The main observation of Nambu and Jona-Lasinio (1961) in their original work 
was that the four fermion interaction generates a dynamical mass for the fermions 
if the coupling constant is bigger than a certain critical value. This was done 
in the canonical formalism using a Bogoliubov-Valantin transformation from bare 
massless quarks to constituent massive quarks. The idea is that if the interaction 
is strong enough, the vacuum lowers its energy creating a mass gap between the 
positive and negative energy continua of the Dirac spectrum. Such a phenomenon 
is called dynamical mass generation. In the following sections we will see explicitly 
how a dynamical mass generation takes place within the path-integral approach 
to the NJL model. 

Spontaneously Broken Chiral Symmetry 

From eq. (2.12) one recognizes at once that for a finite vacuum expectation 
value ay of the a-field the quarks acquire a finite dynamical or constituent mass 
M due to: 



M 



-NJL' 



y = 9^V (2.19) 



dqdq 

which means that chiral symmetry is spontaneously broken (Cheng and Li 1984). 
In order to have a translation and parity invariant vacuum one has to demand that 
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TT = and that ay is independent of x. The value of ay will be determined later 
on. Normal ordering with respect to this vacuum is performed by subtracting the 
vacuum value S^jj{ay, ny) from the effective chiral action eq. (2.16) 

Stationary Phase Approximation 

From now on we will assume classical fields for a and tt or in other words 
we perform a stationary phase approximation in 0th order with respect to the 
stationary points of the effective action 

5Seff{cf,Ti 



6a 
SS^ff{a,Tv) 





(2.20) 





which is the Schwinger-Dyson or gap equation. This means we have now specified 
the approximation in which the model will be treated. It can be shown that the 
stationary phase approximation of the semibosonized version is fully equivalent 
to the Hartree approximation in the original 4-fermion version, which has been 
treated by a number of authors using the Bethe-Salpeter formalism (Bernard et 
al. 1984, Bernard 1986, Ferstl.et al. 1986, Bernard et al.l987, Providencia et al. 
1987, Bernard et al.l988, Klimt et al. 1990) and has been reviewed by Vogl and 
Weise (1991) and Klevansky (1992). Whereas in this chapter we restrict ourselves 
to the vacuum solutions ay and ivy of eq. (2.20), the main aim of this article is 
to study non-translational invariant solutions of this equation corresponding to a 
system with finite baryon number, which will be done in the next chapter. 

Vacuum Expectation Values - Quark Condensate 

In the following we need the expectation value of the bilinear density q{x)ICq{x) 
in the stationary phase approximation, i.e. for classical a and tt, where /C is an 
arbitrary matrix in spin, isospin and color space. This can be expressed by: 

{q{x)ICq{x)} = - 



J VqVqe^ / d*XEqix){-iD)q{a 



5 



-In / VqVqe' 



J d xj^q{x){—iD— K{x)JC)q{x) 



5k{x) J ^ ^ <x)=0 (2.21) 

^-^Spln(-zL> - /t(x)/C)|^(^)^o 

{xE\Tr^Tc[{iD)-^K]\xE) 
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One important special case is the expectation value {qq) , the quark condensate, 
which is an order parameter characterizing the strength of the spontaneous break- 
down of chiral symmetry, like it does M or ay. In the present model it can be 
easily determined from the general expression eq. (2.21) with /C = X by applying 

the Schwinger-Dyson equation (2.20) to eq. (2.16): 

„2 

{qq)v = i^E \Tr[iDv]-^\xE) = (-) — (ay - mo) (2.22) 

2.3. Divergences and Regularization 

UV Divergences 

In order to extract the UV divergent terms in the fermion determinant S^ff {a, if) 
we write its real part with eq. (2.19) in the form: 

ReSef/{a,n) = {-)^S-p\nD^D 

= (-)^Spln[-aij2 + M^ + ig^{a + 2x^75) 

+ (7^((T2 + 7f2-(7y2)] (2.23) 

Introducing the abbreviations: 

G= {-dE^ + M^)-^ 

2 2 ^2 2 (2.24) 

V = ig^{a + iriv'^^) + g {a + tv — ay ) 



F, ^. 1^ - ^ 1. 



we obtain: 

ReSef/{a, n) = -SplnG - -Spln(l + GV) (2.25) 

The first term ^Spln is an infinite constant independent of a and n and vanishes 

after normal ordering. The second term can be expanded in powers of (GV): 

00 -1 

ln(l + GV)=y (-r+^-(GVr (2.26) 

1 n 

747 

In momentum space the functional trace Sp gives an integral / . .^ . Therefore 
we see that all terms with n > 3 are UV convergent. The contributions n = 1 and 
n = 2 can be easily explicitly calculated leading to (Eguchi 1976): 

Re%/(a,ff) = [{-)4NcIi{M)]g^ J d'^XEia^ + if^) 

+ [2iV,//2(M)] [d'^XEia^ + ^2 - 4)2 

^ (2.27) 

+ [(-)4iVe^2^2(M)] f d^XEk[id^a){df,a) + (a^7f)(a^7f)] 



+ Seffi 



0", TT 
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where the integral 

/ (M) = / ^ i r: (2.28) 

is quadratically divergent for k = 1 and logarithmically divergent for k = 2. The 
Sgff{a,n) contains only UV convergent terms and is of higher order in the field 
amplitudes V and their gradients dV. A systematic expansion in {V, dV) can be 
performed by the gradient or heat kernel expansion and will be discussed in the 
next section. 

Regularization Schemes 

The imaginary part ImS^f ^ of the fermionic action in Euclidean space is finite 
because it is in a fundamental way connected with the anomaly structure of the 
theory. The real part ReS'g i- ^ is always divergent corresponding to /i and I2 being 
infinite. 

The divergent integrals /i and I2 have to be regularized consistently by apply- 
ing a regularization scheme to Spjf (a, tv) with a finite UV cutoff A, which enters 
as an additional parameter into our model. This is not considered to be a big 
problem since the model is believed to be an effective low energy approximation 
of QCD, with the cutoff as the relevant scale for low energy hadronic phenomena. 
Physically, the cutoff function models in some sense the gluon cloud around a 
single quark and the correct cutoff function should be derived ideally from QCD 
itself. Because there exists no such derivation up to now one is forced to apply 
general schemes which are used in the literature and fixes the cutoff A in order to 
describe the physics of the mesonic sector as well as possible. In practice there is 
no reason to prefer some scheme in favor of the other provided some constraints 
concerning symmetries are obeyed (Ball 1989). The question of different regular- 
ization schemes has been treated in detail by Meissner Th et al. (1990b) for the 
vacuum and by Blotz et al. (1990) and Doering et al. (1992) in the soliton sectors. 
In this review we will mainly use two of them for practical calculations: Proper 
time ( Schwinger 1951 ) and Pauli-Villars (1949), as it is used in this article, both 
preserving Lorentz-, chiral and vector-gauge invariance (Ball 1989). Actually both 
can be obtained from a generalized proper-time formalism (Ball 1989). 

Gauge Invariant Schemes 
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The generalized proper-time regularization procedure (Ball 1989 ) is based on 
the following identity for the difference of two logarithms 

loga- log/3 = lim / — 0(r, A)(e-™ - e"^^) (2.29) 

A^oo JO r 

with the additional condition ^(r, oo) = 1 and a properly chosen shape for ^(r, A) 
in the limit of small r. Using this representation and assuming that the vacuum 
contribution is subtracted, the generalized proper time regularized effective action 
reads 

logdetDtD - logdetD JDo = |°° — (/)(r, A) [tr e"^°^° - tr e^^^oDo] (2.30) 

The original proper time method of Schwinger corresponds to the choice 

Hr)=0{j^-T) (2.31) 

with ^ as a single step function^^ . The Pauli- Villars scheme is obtained as 

(/)(r) = 1 - (1 + AV)e-^^' (2.32) 

corresponding to usual mass subtractions in Feynman diagrams. 

Some other schemes have been proposed, which however do not show the 
invariances mentioned above and will not be used in this article. Those are e.g. a 
relativistic four dimensional sharp cutoff method 

TrlogiD -^ Tr{^(-a2 - A^)\ogiD} (2.33) 

and also a three dimensional sharp cutoff 

TrlogiD -^ Tr{^(-V2 - h^)\ogiD} (2.34) 

Dynamical Mesonic Terms 

As we can see the logarithmically divergent terms are proportional to the ki- 
netic mesonic part [{d^Gi){d^a) + (9^'7f)(9^7f)] and the mexican hat self interaction 
(cr + yf — (Jy) . Those are exactly the dynamical mesonic terms appearing in 
the linear chiral sigma model (Gell-Mann and Levi 1960). We therefore have gen- 
erated those parts in our model from the fermion determinant S(,ff (a, n). The 
proportionality constants in front of them are dependent on the cutoff A and will 
therefore be determined by the physics of the mesonic sector. 



^ JSee sect. 6 for the double step proper time regularization, which offers some 
more degrees of freedom. 
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2.4. Gradient and Heat Kernel Expansion 

The effective action is a highly non-local functional in the a and tt fields, 
i.e. it involves products of fields at all space-time points. Therefore it is rather 
interesting to investigate its behavior in some limiting cases. In this section we 
will study the case of slowly varying fields, since then the effective action reduces 
to a local Lagrangian. We will also see how kinetic and interacting contributions 
arise in this simple limit as a consequence of the polarization of the Dirac sea. 
This allows also to study some qualitative and quantitative features of the model. 
There are basically two schemes: gradient and heat kernel expansion. 

In essence, the gradient expansion is an expansion in the number of Lorentz 
indices. The problem is merely technical and many methods have been suggested 
(Aitchison and Frazer 1984, 1985a, 1985b). A very elegant and powerful one has 
been proposed by Chan (1985), so that it will be exposed here. The method 
has been worked out assuming the cyclic property to be valid so that the results 
are automatically vector gauge invariant. Chan (1985) has applied it to second 
order elliptic operators up to the fourth order and more recently Caro and Salcedo 
(1993) up to sixth order. Hence it makes sense to apply it to the real part of the 
regularized effective action with 

D^D = -a^ + i^E + E^E 

with E = (/((T + zf7f75). The main idea is to realize that the integral representing 
a one loop Feynman diagram is invariant if all external momenta are shifted by 
the same amount, so that one can as well average over all possible shifts. In the 
proper-time scheme this can be expressed as follows 

ReS^jj = / —(/){t, A)Spexp[-r(-a2 + i^E + ETe) 

The next step is to expand the exponential in powers of the derivative operator d 
up to the desired order and to make use of the cyclic property. This corresponds 
to a derivative or gradient expansion. If one further expands around the vacuum 
configuration S = M one gets the heat kernel expansion. In other words, to get 
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a given order of the derivative expansion one needs an infinite number of terms 
from the heat kernel expansion. The final expressions for the known orders are 
rather lengthy and will not be quoted here. For our purpose it is sufficient and 
instructive to restrict ourselves to the heat kernel expansion up to second order. 
The result is (Kleinert 1976, Ebert and Reinhardt 1986) 

2 

C = Ncg'^I2{M)txr\^J:^^^'t+{t^t-M'^f^+{-^-2Ncg^Il{M))trr{t^t-M'^) 

L ^ J 4^^ 

(2.35) 
The /rj-integrals are given by (2.28). The important point of this expression is that 
even for slightly space-time dependent fields the whole effect can be summarized 
in a kinetic and interaction terms for the fields a and tt. Thus although those 
did not appear explicitly in the original lagrangian, they are indeed present due 
to the polarization of the Dirac sea. In chapter 3 we will show how one can 
go beyond the limit of slowly varying fields by computing the effective action in 
an exact manner. The similarity of eq. (2.35) with the Gell- Mann-Levy (1960) 
sigma model is also very interesting and will be exploited later when fixing the 
parameters. The effective Lagrangian gets even simpler if the so called chiral circle 
condition 

(^ + ^ = In 

is imposed. Then one can use the parameterization 

with a SU(2) unitary matrix U and (p as the non-linearly transforming pion field 
(Coleman et al. 1969). The Lagrangian reads then 

£ = NcM'^I2{M) tYr d^Ud^'U^ (2.36) 

which resembles the Weinberg (1967) non-linear Lagrangian. The former discus- 
sion carries along also for the currents. We just quote the final result for the vector 
and axial currents for linear chiral fields in the same approximation respectively 

1^ = aa^yf - Tvd^a 
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Gradient Expansion of the Imaginary Part - Wess-Zumino-Witten-Term 

The gradient expansion of the imaginary part of the fermion determinant 
eq. (2.18) is more involved. The problem is, that the corresponding term can 
not be written as 4-dimensional space time integral over a local Lagrangian. In 
order to obtain a closed analytical form for the gradient expansion of Im5'g f f one 
considers the change SjjlmSgff (U) under the variation 5U of the chiral field U 
and performs a gradient expansion of this expression. The actual calculation (Dhar 
and Wadia 1984, Dhar et al. 1985, Diakonov et al.l988) follows very closely the 
Goldstone-Wilczek expansion of the anomalous baryon current as it is described 
in App. B. The result is: 



iSulmSeff^iU) 



iNc 



c / A 



grad " 487r2 J ^ ^^^^^^ 



d X€r 



^'^ flavor 



u'fdaU) MdsU) (u'fdjU) MdsU)] SU 



which can be shown to be the variation of 



iNr 



ilmS,ff^{U) 1^^^^ - 24Q^2 Jb^ " --MiM2A^3M4M5 ^2.37) 



d xe, 



^'■flavor 



U^d^.u) (U'fd^^u) (U^d^^u) (U^d^jj) Md^.U 



where B^ denotes a 5-dimensional sphere, whose boundary is the 4-dimensional 
space-time {dB^ = R^). The r.h.s. of eq. (2.37) is the famous Wess-Zumino- 
Witten term, which has been discovered in context of the integration of the chiral 
anomaly (Wess and Zumino 1971, Wess 1972). It turned out that it is indispensable 
to implement it in any effective chiral mesonic model in order to get a proper 
description of low energy mesonic phenomena (Witten 1983a). In case of an SU{2) 
chiral field U it can be shown to vanish identically as it does always, if U is time 
independent. Further details will be discussed in the context of the SU{3) NJL in 
chapter 6. 

Goldstone and Wilczek (1981) have shown that in the lowest non-vanishing 
order of the derivative expansion the baryon current is proportional to the topo- 
logical current 

< S^(x) >= i_eA*i^"/3 trMdr^UU^daUU^dgU) + ■■■ (2.38) 



In the next chapter the relevance of such topological current in the NJL model 
will be discussed. 
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2.5. Mesonic spectra from effective actions 

The calculation of on-shell mesonic two-point functions within the NJL model 
has been undertaken by several authors, both in a pure fermionic language (Blin 
et al. 1988; Bernard et al. 1988; Bernard and Meissner UG 1988; Klimt et al. 
1990; Takizawa et al. 1990) and in a bosonized version ( Alkofer and Zahed 1990; 
Jaminon et al. 1992). As we shall see in Sect. 2.7. this is actually a natural step 
because the parameters of the NJL model will be fixed by reproducing the meson 
masses. In most cases a diagrammatic procedure has been considered within a 
Bethe-Salpeter formalism in the ladder approximation. Similar results have been 
also obtained in a time-dependent Hartree-Fock approach (da Providencia 1987 ). 
The virtue of the on-shell definition method is that no approximation is involved 
apart from working in the leading order of the large Nc expansion. As it can be 
shown, the heat kernel and derivative expansions are low energy approximations to 
the complete on-shell two-point function. In this section we sketch the calculation 
of mesonic spectra in the bosonized version of the model as it has been done 
by Jaminon et al. (1992). This is based on the observation, that the relevant 
quadratic part S^rF of the effective action (2.23) for SU(2) scalar fields can be 
always written in the form (Jaminon et al. 1992) 



5,;/'2P(a,7r) 



d q 

(2^ 



{a{q)a{-q)Za{q) [q'^ + ml{q) 
n{q)n{-q)Zn{q) \q'^ + ml{q) 



+ 



(2.39) 



Therefore one can define mesonic masses from 

1 1 52^,;/'2p(a,7r) 



Za{qV S{P1 - P2) S<^{Pl)S(^{P2) 



1 



1 



d^S. 



F, 



eff 



2p( 



p=pi=~P2,Pi^q'^ 



cr, TT 



P=Pi=-P2,p'^=q'^ 



^7r(?2) S{pi - p2) Sn{pi)Sn{p2) 

which effectively corresponds to a rescaling of the fields according to 

^71/2/ 2n / ^1/2/ 2 

Zcr' [q ) (J, Ti = Zt,' [q 



p^^q^ 



(j? + ml] 

[p^ + ml) . i2.40) 



p-=q'- 



a 



n 



(2.41) 



The on-shell definitions of the mesons can be found then by evaluation of- Zfj{—q 
ml) and Z-r^^—cf = ml) to obtain ml = ml{—q = ml) and ml = m,l{—q 



25 



rn^). where the Z-factors are given by the Feynman integrals of the form (Janiinon 
and Ripka 1993) 

Z^{q^) = 4N,g^ f -^G{k + q/2)G{k - q/2) (2.42) 

with the propagator G{k) = l/{k + M ). From these expressions it is clear that 
the corresponding off-shell definitions can be obtained by setting g = in all 
the equations. These expressions then coincide of course with the corresponding 
one, which one would obtain from a gradient expansion of S^jf from the very 
beginning. Then the approximated analogue of (2.39) is already given by (2.27). 
The field redefinitions become trivial in this case and the mesonic spectra could 
be read off directly from the effective potential (Coleman 1985). For definiteness 
we present Z^{q ) explicitly for the proper time regularization eq. (2.31) 

Z^q^) = ^^Jdxr{0, [m2 - x{l - x)q^]/A^) (2.43) 

where we made use of the incomplete r(n,x) -function 

r{n,x):= dtf-^e-"^, a; > (2.44) 

Jx 

2.6. Fixing of the Parameters in the Mesonic Sector 

We have mentioned already that if the coupling constant G is higher than a 
certain critical value, the vacuum is unstable if the mass is increased. This leads 
eventually to the Goldstone phase. Assuming that this is the case we will fix the 
parameters of the model. Two ways for parameter fixing may be distinguished 
depending upon the particular approximation involved. We will refer to them 
for brevity as the on-shell and the off-shell condition. The off-shell fixing of the 
parameters is dicussed in detail by Meissner Th et al. (1990b). 

Mesonic Properties 

Looking at eqs. (2.16,2.19) we see that we have rMo,/U , A, ay and the con- 
stituent quark mass M as the five parameters of the model. Thus one can impose 
e.g. the following four conditions 

1) The stationary phase condition SS^j-f/da in the vacuum state 



6S. 



eff 



5a 



= (i? - SNcg^hiM) - ^^^^ ay = (2.45) 

vac \ gay J 
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immediately determines n in the broken phase with ay 7^ 0. 

2) Fixing of the pion decay constant f-n- is done by considering the expectation 
value of the axial current between a pion state and the vacuum, which in the 
semibosonized form corresponds to consider 

< I A^(x) I tt'^ >= A^ fvqVq [q{x)'y ^'y^r'' q{x)] e' -f '^^''^'njl. (2.46) 

where N is the normalization constant. This is a similiar like eq. (2.21) but 
now it has to be evaluated for a plane- wave pion state and not the vacuum. 
In the one-fermion loop approximation we obtain then 

< I A«(a;) I 7r« >= avZi^\q\^,n^ix). (2.47) 

Comparison with the known matrix element for the pion-decay from current 
algebra (Cheng and Li 1984) 

< I A«(x) I tt'^ >= ip^Mq'^ = -ml)n''{x) (2.48) 

gives immediately the relation 

Mq^)=avZl/\q^) (2.49) 

for Zty and ay. 

3) From the definition of the pion mass in eq. (2.40) one obtains straightfor- 
wardly the explicit form 

mi = ^^ (2.50) 

gay Ztt 

4) Now one can either choose the normalization of the pion wave function 
eq. (2.41) according to 

Z^{q^ = -ml) = 1 (2.51) 

or one can equivalently require the classical PC AC relation (Goldberger and 
Treiman, 1958) d^^ A?^ = /jr ?Ti:^7r", which in the present model gives in 
analogy to eq. (2.8) but now from eq. (2.12): 

gl^^a ^ am^^a ^2.52) 
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so that together with eq. (2.50) and eq. (2.49) again eq. (2.51) follows. 

This fixes all the parameters except one. Usually the constituent quark mass 
M is chosen as the free parameter. Then from eq. (2.51), which determines the 
cutoff A for a given constituent quark mass M, and from eq. (2.49), one can deduce 
c^y = /tt, so that eq. (2.50) inserted in eq. (2.45) gives ^ . Using again eq. (2.50) 
gives the current quark mass rriQ. So M = qJt^ is the only free parameter of the 
model. 

The above description refers to the on-shell definition of the parameters, whereas 
the off-shell formulation corresponds to setting g = in eq. (2.51). Because the 
Goldstone bosons and especially the pions have a small mass compared to the con- 
stituent quark mass, the difference between both methods is rather small. However 
this is no longer the case, when heavier mesons like p, uj and Ai are taken into 
account. 

Furthemore, using again eq. (2.22) and eq. (2.50), it follows the famous Gell- 
Mann, Oakes and Renner (1968) relation 

<qq>mo = -mlfl + C(mo) (2.53) 

Numerical Results 

As we have said already the off-shell definition is only a low momentum ap- 
proximation to the on-shell condition. For on-shell pions the external momenta are 
rather small so that the numerical deviations from the two methods turn out to 
be less than 7%. The only important difference is that in the chiral limit (toq = 0) 
the scalar meson mass acquires a value ma = 2M in the gradient expansion, while 
in the on-shell method the a meson lies in the quark-antiquark continuum with a 
finite decay width into this channel. This is the first manifestation of the lack of 
confinement within the model. 

The dependence of A against M can be seen at Fig. 2.1. In the relevant region 
of M = 350 — 450MeV, the numerical value of the cutoff A depends noticeably on 
the regularization scheme chosen. 

From eq. (2.22) we can now evaluate the quark condensate and compare it 

with the generally accepted value (Shuryak 1986) of 

-I 

< -{uu + dd) >^/3= -(225 ± 24)MeV (2.54) 
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The results for the vacuum sector, i.e. quark condensate, current quark mass and 
the vacuum energy density (de Grand et al. 1975) can be seen at Tab. 2.1. For 
these observables, i.e. for the quark condensate and the current quark mass a 
plateau has been observed for constituent masses above 400MeV. 
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3. Systems with Finite Baryon Numbers - Solitonic Solutions 

In this chapter we consider static (time independent) meson field configura- 
tions. First we show the expression for the total energy (sect. 3.1). A system with 
baryon number S = 1 is obtained by adding Nc = 3 valence quarks, which is for- 
mally done by introducing a thermochemical potential /U. Furthermore we discuss 
the connection between baryon number B and topological winding number n of 
the n field (sect. 3. 2). The mean field equations of motion for the non-linear model 
(meson fields restricted to the chiral circle) are derived and solved for a B = 1 
system (sect. 3. 3). Systems with higher baryon and winding number are briefiy dis- 
cussed (sect. 3. 4). Finally we show that for the linear Nambu-Jona-Lasinio-model 
(off the chiral circle) no solitonic solution exists but the system collapses to a zero 
energy and zero size configuration (sect.3.5). 

3.1. Static Meson Field Configurations - Energy and Static Expec- 
tation Values 

1-Particle Hamiltonian h and Dirac Spectrum 

In order to construct a system with finite baryon number B one has to consider 
meson field configurations different from the vacuum. Because we are interested in 
static properties of the baryon, we will restrict ourselves to time independent meson 
fields a and n. For those the Euclidean Dirac operator iD can be separated into a 
trivial time derivative -^ as well as the time-independent 1-quark hamiltonian h: 



-iD = P{^ + h) (3.1) 



where 



and in the vacuum: 



h = — +gP [a{r) + zf7f(r)75] (3-2) 



hv = ^+l3M (3.3) 



It is our aim to express the operator Spln(— zD) through the eigenvalues of h. One 
should note that h is hermitian and traceless: Tth = 0. Furthermore we make the 
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hedgehog ansatz for a and yf: 

a{r) = a{r) 

(3-4) 
7r(r) = f7r(r) 

For chiral models which use a gradient or heat kernel expansion for the highly 
non local fermion determinant Spln(—iD) like the Skyrme model (cf. sect. 8.1) or 
the Gell-Mann-Levi chiral sigma model (cf. sect. 8. 2) it has been proven (Ruiz- 
Arriola et al. 1989) that the hedgehog shape eq. (3.4) is a necessary condition 
for the meson fields to be a solution of the time independent classical mean field 
equations of motion for a baryonic system. One has to admit that there exists up 
to now no such proof for the case if Spln(— zD) is treated exactly as it is done in 
the present approach. 

Due to the hedgehog ansatz (3.4) the 1-particle hamiltonian h (3.2) commutes with 
the Grand spin G (where G = J + T, J:total spin, T: isospin), its ^-component 
Gz as well as the parity 11 and the 4 observables (/i,G jG^^n) form a complete set 
of commuting operators. Therefore the eigenstates | A) of h: 






(3.5) 



can be characterized by the 4 quantum numbers ex.G.G^ and (— ) , where G^ is 
degenerated. Following Kahana and Ripka (1984) the eigenvalue problem eq. (3.5) 
can be solved numerically by putting the system in a large but finite sphere with 
radius D and appropriate boundary conditions for the radial part of the free eigen- 
functions at D. In doing so one obtains a basis consisting of discretized wave num- 
bers km which belong to a given grand spin G and parity (— ) . For the numerical 
diagonalization of the matrix of h one takes wave numbers smaller than a given 
numerical cutoff frequency Kmax into account. Kmax and D have to be chosen 
large enough so that the value of any calculated observable does not change any 
more by further increasing Kmax and D. It should be emphasized that K^ax is 
a purely numerical cutoff and has nothing to do with the model intrinsic physical 
UV cutoff A introduced and fixed in Chap. 2 in order to reproduce the pion decay 
constant /tt. 

Fig. 3.1 shows the spectrum of h for a meson profile on the chiral circle (non-linear 
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realization of chiral symmetry) as example: 

<j{rf + At? = al + nl = fl 

a{r) = U cos e{r) (3.6) 

7r(r) = /tt sin 6'(r) 
where the chiral angle 9{r) is parameterized by 

0{r) = -rrne'Tl (3.7) 

with the topological winding number n = 1 (Kahana et al.l984). The parameter R 
characterizes the size of the meson profile and in a way describes the magnitude of 
the deviation between the actual and the the vacuum field configuration (i? = 0). 
If R is large enough (R-M > 0.5) one finds bound orbitals both in the positive and 
the negative spectrum. The positive bound state with lowest single particle energy 
and quantum numbers G = 0+, will be called the valence orbit | A) = | val). Its 
single particle energy decreases with increasing R, switches sign and gets finally 
part of the negative spectrum. This level originates from the positive continuum 
and gets bound and localized by interacting with the negative continuum (Dirac 
sea). This simple non-linear mechanism creates the soliton. 

Static Energy 

Because the meson fields are time independent, the trace over the Euclidean time 
coordinate r in the fermion determinant can be performed trivially. Assuming anti- 



periodic boundary conditions in the Euclidean time interval 



r r 



the Euclidean 



time gradient -^ has the spectrum Un = -^-^ — -'re , n = 0, ±1, ±2, . . .. Therefore 
the expression for any functional F of the form F -g-, f{r) simplifies in the zero 
temperature limit (T ^ oo , Aun -^ 0) , which leads to the the ground state of the 
system, to: 

d 



SpF, 



^Jirl 



^T j SxE j '^TrF^[-iu, f{r)] (3.8) 

Especially if we take for F the effective action S^ff we find for the imaginary part: 

Im5'ejj((T, 7?) = Im5'ejj^(cr, 7f) 

/ N v^ /■ '^^T 1 / • , y-^ f du — w (3.9) 

= (-) Z^ / 7^Imln(-zw + ca) = 2^ / 7^ arctan = 
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which is based on the fact that h is hermitian and therefore e;^ real. As we have 
already mentioned this is not true if cu-mesons are included in the model (cf. chap. 
7). We therefore have: 



Seff{(y, 7?) = Re5'e//(cr, Tv) = T ■ E 



(3.10) 



where E is the total static energy of the system. The explicit expressions for the 
various parts of E in terms of the eigenvalues e;^ read: 



J^ ^Sea T~ tjrnpg -\- tj] 



br 



(3.11a) 



with: 



E 



sea 



Nc 
' 2 



E7^l(eA,A)-E^l(^Av.A) 

.A Ay 



(3.116) 



Ernes 
Ebr 



2 



j d^r{a'^ + n^ - fl) (3.11c) 

-mlU j d\{^ - U) (3.11rf) 



In case of the proper-time regularization method, which we will restrict ourselves 
on for explicitness in the following, the regularization function TZi is given by 
(Meissner Th et al.l988): 

A 



7^l(eA,A) = (-)■ 



'An Jl 



dss-''/\--<'^/^y 



(3.12) 



A detailed and complete discussion of the various regularization schemes can be 
found in Blotz et al. (1990) and in Doering et al. (1992). Furthermore it should 
be noted, that because of Tr/i = in the unregularized case (A -^ oo) eq. (3.116) 
with eq. (3.12) can be written as sum over the negative (occupied) states 



E, 



E ^A 
A,eA<0 



E 



fA 



V 



(3.13) 



A\/jeAy<0 



which corresponds to the naive picture of the occupied negative states (Dirac sea) 
(cf.Fig. 1.3). If a regularization is performed tr h = is only approximately true. 
Since its numerical value in realistic cases of baryon number S = 1 is about 50MeV 
it actually does not matter if one runs over the whole spectrum or only over the 
occupied states. 
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Static Qbservables 

Finally we consider the expectation value of the static observable O := J d rq^{r)Oq{r)^ 
where O represents an arbitrary time independent operator in Dirac- and isospin 
space. Because of O = :^ /(i x^(x)/3(9q'(x) we have after Wick rotation into Eu- 
clidean space: 



E 1 JI^Ql^Qld^^E \q{x)PO^q{x) 



J d^XEqix){-iD)q{x) 



I sea 



^ j'-pg'Dge^/^^^^^(^)(^^^)'^*^^) (3 14) 



TduE 
Using the relation 



^ ^ Spln(-^L> - w^/JO^; 



„ — vac. contr. 



f+k du 1 1. 

iim / = -signx (3.15) 



one gets in the unregularized case: 

{0^)sea = ^NcY.{-signex)Of - vac. contr. (3.16a) 

A 

with: 

of := I d^r(j)\{r)0^(j)x{r) = j d\^x{r)f50(j)x{r) (3.166) 

In the regularized case one has to substitute [Spin] -^ [Spln]re.g, which gives after 
rotating back into Minkowski space: 

{0)sea = (+iVc) E ^2(eA, A) ' ^A " vac. contr. (3.17) 

A 

where for the proper time method the explicit expression for the regularization 
function reads: 

1 r- , 1/2 (ex\ -sie,/A)^ 



7^2(e,, A) = (-)-= y^ rf..-V^ (^lA ) e-n^A/Aj (g^^g^ 
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3.2. Baryon- and Topological Winding Number - Chemical Potential 
- Valencequarks 

Baryon Number and Baryon Density 



The baryon current bn is generally defined as: 



1 



b^{x) := —q{x)-f^q{x) 



(3.19) 



Its 4-divergence vanishes 9^6^ = due to the t/(l)-symmetry [q h^ e^'^q) of the 
Lagrangian giving rise to a conserved charge, the baryon number 



B= f d^rboix) 



(3.20) 



Applying eqs. (3.16, 3.17) for (9 = X we obtain for the expectation value of B in 
case of time independent fields without regularization: 



{B) = ^Nc 



^(-signcA) - 5I(-signeA, 



A Ay 

For the baryonic density 6o(r) we find correspondingly: 



(3.21) 



(&0(r1) 



1 



-Nr. 



^(-signeA) 
A 



<P\ir)<P\ir) 



^(-signcA^; 
Ay 



■</'L(^'^Ay(^ 



(3.22) 

Neither {B) nor (6o(r)) are UV divergent. For (B) this is clear from eq. (3.21) 
whereas for (6o(^)) it follows from the fact that the 2nd order gradient expansion, 
which gives rise to the logarithmic divergence, vanishes and the 4th order is already 
UV convergent (cf. app. B). 

Therefore there is no need for a regularization of the baryon number. Further- 
more as we see at the end of the next section, this is consistent with a more refined 
derivation of the baryon number, which in that case originates from the imaginary 
part of the Euclidean effective action. It is therefore consistent to start with a 
proper time regularized real part of the effective action as it is done in Chap. 2 
and treat the baryon number unregularized, as long as one treats the imaginary 
part unregularized. However, one should mention that recently by Schlienz et al. 
(1993) also the imaginary part was regularized though this treatment does not 
preserve the anomaly structure of the theory. 
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Furthermore it should be noticed that due to the general method of calculat- 
ing static expectation values of observables as it has been described in the last 
section, our expressions for {B) and (69 (r)) differ from the usual ones obtained in 
Minkowski space, where the sum J2\^s performed only over the negative (occupied) 
states. For the global quantity (B) it is easy to see that both cases lead to the 
same result, whereas the local densities {bQ{f)) may differ. 

Thermochemical Potential n; Separation of the Valence Part 
In order to constrain the baryon number of the system to a given value (e.g. 
(B) = 1 for nucleons and hyperons) we proceed the way known from statistical 
mechanics (Huang 1987, Negele and Orland 1987), and introduce the thermochem- 
ical potential /U as a Lagrange multiplier into the generating functional Z which 
becomes now the grand canonical sum of states Z{p) (Williams and Cahill 1983, 
Meissner Th et al.l990) 

Z{fi) = J VqVqe~ ^ d^^Eq{~iD~t^P)q ^323) 

After integrating over the fermion fields we arrive at the fermionic part of the 
grand canonical effective action 

S^jj{a, 7f, ^i) = {-)Sp\n{-iD - ^(3) (3.24) 

The mesonic part remains unaffected. After subtracting the vacuum contribution 
the SKr can be splitted in a natural way into 2 parts: 

S^ff{a, TT, Ij) - Spjj{ay, TTy, ^y = 0) = 

^Slfja,n,f^) + S'^%a,n) (3.25) 



where: 



and 



Qval 
^eff 



Seffi^^^T^^IJ')- S^ff{a,n,n = 0)\ (3.26) 



Qsea . 

^eff - 



'^F ( ^ z± .. r\\ oF 



S^ffia, n,ij = 0)- S^ff{ay, ny, jj, = 0)J (3.27) 

The S^f^Aa, yf) is of course nothing but the fermionic contribution (Dirac sea) to 



the effective chiral action S^rr{a,n) = TEsea considered in the last section. In 
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general it needs regularization. The valence contribution S'^'iAa, tt, h) is finite and 
needs not to be regularized. It can be expressed in terms of the eigenvalues e;^ by: 

Slfj{a,TT,^i)=N,T Y. (^A-^) (3.28) 

0<eA</" 

where we have used the relation: 

ln(— i-u -\- 1\ — n) — ln(— m + e;^) = / dx — -. (3.29) 

as well as eq. (3.15). Similarly one can separate the grand canonical expectation 
value of any observable O := J d rq'Oq into a valence and a sea part: 

{0){fi) = {0),,iifi) + {0)sea (3.30) 

with 

(O)„a/(^) = (O)(^)-(O)(^ = 0)=iV, Yl Ox (3.31) 

0<ex<fi 

Especially for the baryon number (B), which is the integral of the time component 
of the baryon current bf^ = {i/Nc)q'yj]^q in Euclidean space, one obtains with 
setting O = 74 in eq. (3.14) and using 6q = —ib^ (cf. app.A): 

d ^ , / .„ ^„^ 



with 



and 



(S)(^) = ^^ ■ ^Splog [-zD - UE^O'' - li 

_ 1 r+oo rfw 1 _ 

^ Wc J-00 2^1-^ u + i{ex -fi)^ ^^'^^^ 

= i-)hYl sign(eA - //) = By^ii^) + Bgea 
A 



(5)mK^)= E 1 (3-33) 

0<eA<M 



(5).ea = (-)^EsigneA (3.34) 

^ A 

Analogous expressions hold for the baryon density. 

In any case one notices that the valence contribution is finite and that there is 
no need for a regularization of this part. It is also important to realize that 
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the thermochemical potential /U is a real number both in the Minkowski space 
and in the Euclidean space, which causes finally simply a shift in the eigenvalues 
e;^ ^ e;^ — fi. It is therefore different from the treatment of the time component 
of a four vector uj^, which has to be Wick rotated from Minkowski into Euclidean 
space like the real time: ujq -^ iuj^, where ujq and 0^4 are both real (cf. App.A). 
From that and from the antihermiticity of 74 it is clear from eq. (3.32) that the 
baryon number B{p) originates from the imaginary part of the effective Euclidean 
action. Therefore the baryon number is finite and need no regularization. This 
philosophy will be prosecuted throughout this work. Actually the imaginary part 
is not zero in the present case due to the additional term: ujj^/30 

Using eqs. (3.32) and (3.28) we can write the valence part of the grand canonical 
effective chiral action as: 



^S^eff = Nc 



(3.35) 



J2 ^A-^S^aK^) 
LO<eA<M 

As it is known from statistical mechanics the total energy of a system at zero 
temperature T ^ 00 (ground state) is given by: 

Etoti^^) = ^S.ffi^i) + N,ii{B){ij) (3.36) 

which leads finally to: 

Etotip) = Esea + <^val{p) + ^mes + -£^5r (3.37a) 

with eq. (6) and 

et-oz(^) =Nc Y. ^A = Nce^al (3.376) 

0<eA<M 

Valence Picture and Bosonized Picture 

From Fig. 3.1 we can distinguish 3 different regions concerning the behaviour 
of the valence level O^ and, associated to it, the behavior of E and {B) {Emes = 
Ebr = 0)7 which have been discussed already by Meissner Th et al. (1990a). 

(A): 

(^val > ^ > 

(5)(^)=0 ^ ^ 

(3.38) 



Etot = E, 



sea 
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(B): 



(C): 



e,ali^^)=Nc E ^X = Nceyal (3-39) 

0<ex<fi 

Etot = Esea + Nc E e;^ = E'gga + ^c^val 
0<eA<M 



(3.40) 



-E'toi = E, 



sea 



Case (A) is uninteresting for our purpose, because the valence particle is part of 
the positive continuum and {B) = 0. The cases (B) and (C) give (B) = 1 as 
desired. Hereby the valence orbit is counted explicitly for EtotilA as well as for all 
other static observables if e^^^ > (B)- If ^val < ^^^ valence particle gets part 
of the negative spectrum (Dirac sea) which now as itself carries baryon number 
(-B)(|u) = {B)sea = 1 due to the fact that it contains now 1 orbital more. For the 
following it is convenient to introduce the notation: 

Vval = „ ., ^ „ (3.41) 

I , if e^al < 

so that we can write for a system with B = 1: 

Eb=1 = NcVval^^val + Esea + Emes + E^^. 

(3.42) 

(0)s=i = iVcr?w(0)w + (0).ea 

Fig. 3.2 shows eyg^i^Egea and Etot {Emes = Efjr = 0) for a meson profile with linear 
shape: 

on the chiral circle with winding number n — 1, M = 465MeV and TOtt = in 
dependence of the size parameter R. One clearly recognizes that Etot{R) has a 
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local minimuni at i? ~ 0.8fm corresponding to a solitonic solution with B = 1. 
In the next section we will construct these solutions selfconsistently by variation 
of Etot[(^{r)] with respect to all degrees of freedom {9{r)}. On a first glance to 
eqs. (3.41,3.42) it appears as if observables change discontinously their value if 
the valence single particle energy changes sign. In fact this is not true. The 
regularization function takes care of this automatically and indeed one can show 
analytically and numerically that no discontinuity and no kink appears (see Fig. 
3.2. as examples). 

Connection between Baryon Number B and Topological Winding Number n; 

Gradient Expansion 

From Fig. 3.1 one can see that for very large values of the size parameter R the 
valence orbital approaches the states emerging from the negative continuum and at 
the end cannot be distinguished any more from them (Kahana et al.l984, Kahana 
and Ripka 1984). On the other hand we know that for those large profile sizes R 
the gradients of the meson fields get small drcr oc 9r7r oc -^ so that the gradient or 
heat kernel expansion, discussed in sect. 2. should be valid. The gradient expansion 
of the baryon current has been performed by Goldstone and Wilczek (1981) and 
is shown in appendix B. Indeed it turns out that the gradient expanded baryon 
number of the Dirac sea {B)sea is identical to the topological winding number n 
of the meson profile 

lim {B)sea = n (3.44) 

which is n = 1 in the case considered above. This feature generally holds in the 
region (C). Therefore in the case for large R the valence particle has got part 
of the negative spectrum and the baryon number coincides with the topological 
winding number. Similar considerations hold for higher winding numbers (Kahana 
et al.l984). For sufficiently large R one therefore gets close to the philosophy of the 
topological soliton models like e.g. the Skyrme model which contains no valence 
quarks but relates the baryon number of the soliton purely to the topological 
winding number of the Goldstone field. Furthermore it turns out, that for large R 
the energy of the Dirac sea Esea{R) approaches the kinetic energy of the mesons, 
which is the leading order in the gradient expansion (Meissner Th et al.l988): 

lim E,ea{R) = Ek,^{R) = \ld^r \{Vaf + (W)^] oc R (3.45) 
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as long as this expansion is convergent or at least an asymptotic series in -n (Zuk 
and Adjali 1992). 

3.3. SoLiTONic Solutions of the Nonlinear Model 

Mean Field Equations 

The equations of motion for a system with baryon number B = 1 are given by the 
stationary points of the grand canonical effective chiral action Sgff{p, B = 1). In 
practice the hedgehog ansatz eq. (3.4) is used and the system is assumed to be 
constrained to the chiral circle (3.6). The restriction to the chiral circle is quite 
essential and its origin will be discussed in sect. 3.5 and 8.4 (Sieber et al. 1992; 
Meissner Th et al. 1993; Weiss et al. 1993). The stationary points of the grand 
canonical effective chiral action 5'gjj(//, B = 1) (see eq. (3.25)) with respect to the 
chiral angle 9{r) are then given by: 

S{eir)}SefMB = l) = (3.46) 

which reduces in case of time independent meson fields to the variation of the 
energy 

5{^(,)}i?(^,S = l) = (3.47) 

In order to perform the variation we start from eqs. (3.42, 3.11) and use the spectral 
representation 

ex = (A \h\ X) = I dW^{f)h(Px{f^ (3.48) 

Because the variation 5(A|A) = 6 J d r(p'^{r)(px{r) vanishes we get for the variation 
of the single particle energy e;^: 

,2 



(3.50) 



^{e{r)} =9 drr^[{- sin e{r)) ■ sx{r) + (cos 6'(r)) ■pxir)]S6{r) (3.49) 

with 

sx{r)=ldn^xir.n)cPx{r,n) 

Pxir) = J dQcpxir, 9)i'^^{Tf)(j)x{r, O) 
which lead to the equations of motion: 

sin 6'(r) [Ncg {S{){r) + Vval^vali^)) - 47r/^mJ = 
cos^(r) [Ncg {Po{r) + VvalPvalir))] 
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(3.51) 



where: 



5o(r)=E^2(eA,A)«A( 



r 



•^ (3 52) 

^0(r) = E^2(eA,A)PAW 
A 

and the regularization function 7^2 (^A 7 A) has been defined in eq. (3.18). 

Numerical Treatment 

Eq.(3.51) has been solved numerically by Reinhardt and Wuensch (1988) as well 
as Meissner Th et al.(1989), who used a standard selfconsistent procedure as it is 
known from Hartree and Hartree-Fock calculations in atomic and nuclear physics. 
One starts with a reasonably chosen profile 9Q{r) given e.g. by eq. (3.43), diagonal- 
izes the h as described in sect. 3.1, calculates sx,px,sq,pq from the eigenfunctions 
4>x{r) and the eigenvalues e;^ from eq. (3.52) and obtains a new profile function 
6*1 (r) from the equation of motion (3.51). The procedure is iterated until a desired 
degree of selfconsistency is reached. An equivalent method for solving eq. (3.51) 
consists in parameterizing 9{r) in terms of a parameter set {aj^ , /c = 0, 1, 2, , . . .} 
covering the whole r-dependence of 9{r) as fully as possible and performing a 
minimization of the total energy in the {a^j-space. First calculations (Diakonov 
et al.l988, Meissner Th et al.l988) used a simple, one dimensional parameteriza- 
tion in terms of the size R of an appropriate chosen meson profile 9 [i] like it 
was done in eqs. (3.7,3.43) and minimized just with respect to R (cf. Fig. 3.2). 
The accuracy of this very rough ansatz was improved by increasing the number 
of parameters to k = 3 (Diakonov. et al.l989), which already gives results close 
to those of the selfconsistent method. In a very recent work 9{r) is calculated 
at a given number /c = 10 — 30 meshpoints in a sufficiently large interval [0, D] 
and a systematic search for the minimum of Efot in the /c— dimensional space 
spanned by these meshpoints is performed using an elaborate minimization al- 
gorithm (Sieber.et al.l992). In other approaches (Diakonov. et al.l988, Adjali et 
al. 1991, 1992) the energy of the Dirac sea Esea[9] (eq. (3.116)) is approximated 
by expanding SplnD'D (iD = i^ — g{a + Itti"^^)) in terms of the 'perturbation' 
g^{a + iTn'j^) like it was done in sect. 2. 3. The corresponding expression is exact in 
the case of a very large meson profile size as well as in the case of small deviations 
from the vacuum configuration (a = ay = /jr, tt = 0). In doing so the numerical 
diagonalization of h is avoided and the variational problem reduces to the solution 
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of an integro-differential equation in 9. 

Generally the variation is performed for a given value of the constituent quark mass 
M or equivalently the ratio A = -^, which is uniquely related to M by (2.51) and 
(2.45). Furthermore finite pion masses (tott = 139MeV) as well as the chiral limit 
(tott = 0) have been considered. It turns out, that solitonic solutions of eq. (3.51) 
exist if and only if M exceeds an critical value M > M^r or equivalently A < Acr- 
The corresponding values are shown in Tab. 3.1. This cusp behavior is typical for 
localized (solitonic solutions) of a system of coupled non linear equations (see e.g. 
Lee 1981 and ref. therein) and has been observed also in other chiral quark meson 
models (Birse 1990 and ref. therein). 

Selfconsistent Fields; Mean Field Quantities 

Fig. 3.3 and Fig. 3.4 show the total {mean field) energy E^^^p = Etot as well as 
the isoscalar electric quadratic radius {R ) = An J drr bQ{r) which is nothing but 
the sum of corresponding proton and neutron contributions, in dependence of the 
constituent quark mass M for rnvr = (Meissner Th et al. 1989, Meissner Th and 
Goeke 1990). Both quantities are divided into their valence and sea contributions, 
respectively. One should notice that the mesonic self energy Emes (eq. (3.11c)) 
vanishes due to the chiral circle constraint eq. (3.6). Using m-jY = 139MeV one 
finds that the chiral breaking term E^,^ (eq. (3. lid)) is small and its influence on 
the solitonic solutions is in general rather small of less than 10 — 15% (Meissner 
Th et al.l989, Meissner Th and Goeke 1991). Apparently it affects the form of 
the asymptotic behavior of of the pion tail of the selfconsistent solution 

lim 7r(r) = Ae^^-^^^^^ (3.53) 

Thus it has some influence on quantities which are sensitive to that form, especially 
if they diverge in the chiral limit m-ji -^ 0, like e.g. the magnetic polarizability and 
the E commutator. Actually eq. (3.53) can be obtained by substituting the sea 
energy Egea (eq. (3.116)) in the variational principle (eq. (3.47)) through the gradi- 
ent expanded expression. One has then to solve the resulting differential equation 
in the coordinate space. It contains formally an infinite number of derivatives for 
which in the asymptotic region (r -^ 00), only the contribution from the kinetic 
term (^) survives. 
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In Fig. 3.5 and Fig. 3.6 the meson profiles ^{r) and n{r) of the selfconsistent 
solutions of eq. (3.51) for three representative values of the constituent quark 
mass M with mjr = are presented. The corresponding baryon densities are 
shown in Fig. 3.7. As one can see the actual form of a{r) and n{r) as function 
of the physical distance r is quite independent of M. If we denote with R some 
characteristic size of the profile, we find that the condition 

l^ral ^ \dr7i\ < M ^ 1<^M-R (3.54) 

for the validity of the gradient expansion is fulfilled if M > lOOOMeV and therefore 
in this case the baryon number {B)sea of the Dirac sea and the topological winding 
number n of the pion field coincide (cf.app. B). Indeed as one can see from 
Fig. 3.5 the valence energy becomes negative e^^; < if M > 780MeV. Due 
to the discussion in sect. 3. 2 one therefore comes close to the philosophy of the 
topological soliton models if M gets large (cf.Fig. 1.4), whereas the valence quark 
picture holds for small M (cf.Fig. 1.3). The important point is the fact, that 
for a given M the model can decide between the two pictures dynamically by the 
equations of motion eq. (3.51) and therefore switch continuously between them 
just by the variation of M. In order to make a definite decision which picture 
is favored one has to calculate baryonic observables and see for which values of 
M they are described properly. If one looks e.g. at the isoscalar charge radius 
{R ) = {R )p + {R )n (Fig. 3.6) one recognizes at once that one is restricted to 
rather small constituent quark masses (M ^ 360MeV) just above the critical cusp 
Mcr in order to reproduce the experimental value of 0.63fm . The corresponding 
mean field energy lies at about E]\jp ^ 1200MeV. High values of M corresponding 
to the bosonized picture lead to a baryonic system whose extension is by far too 
small. The region M ^ 350 . . . 450MeV is also consistent with that of various 
non relativistic constituent quark models (see.e.g.Bhaduri 1988 and ref. therein). 
Furthermore one can see from Fig. 3.6 that the contribution of the valence quarks 
is clearly the dominating one. Counterexamples as e.g. the neutron squared charge 
radius or the nucleon polarizability are found as well. We will see in chapter 5 
that this will be true for most of the nucleon observables. If and how this feature 
changes if vector mesons {uj,p,Ai) are explicitly coupled will be discussed in chapter 
7. 
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3.4. Higher Winding and Baryon Numbers 

For the construction of the solutions of the mean-field equation of motion for 
9{r), one has assumed that 9(0) = —mv, with the topological winding number 
n = 1. There is, however a priori no reason in this model, why n should be 
fixed to 1, because the baryon number B is carried by the valence quarks, if one 
chooses an appropriate thermochemical potential fu.. This is different from the 
topological soliton models (e.g. the Skyrme model), where B = n from the very 
beginning. There is even no a priori reason why n should be an integer number, 
if one leaves mathematical arguments like continuity of the tt field in the origin or 
the convergence of the gradient expansion for B (cf. app. B) aside. In contrast to 
to the Skyrme model, where n & Z is required in order to obtain a finite energy, 
the total energy Etot is finite for any value of n in the present approach. Diakonov 
et al.( 1988, 1989) and Berg et al.(1992) have investigated the behavior of Etotip) 
as a function of n for a fixed profile form. It turned out, that the minimum 
of Etot definitely lies at the point n = 1, which means in fact that this feature 
is also a dynamical consequence of the equations of motion. Unfortunately the 
whole analysis is based purely on numerical arguments and there is up to now no 
analytic proof of this fact. 

In addition Berg et al.(1992) studied systems with higher baryon numbers B > 1. 
They found that local minima corresponding to solitonic solutions only exist if 
n < B and in any case these local minima appear for integer winding number n. 
For example for the B = 2 system they found a soliton in the n = 1 sector with 2 
valence orbitals (0^",0^) occupied if M > 370MeV. In the n = 2 sector solitonic 
solutions exist if M > 560MeV. The main difference between the two cases is the 
fact, that for n = 2 the 0^ as well as the 0^ orbital cross the zero line and get 
part of the negative spectrum whereas in the case of n = 1 only the 0^ orbital 
comes down and the 0^ orbital remains in the positive spectrum. (cf.Fig. 3.1). 
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3.5. Collapse in the Linear Model 

Up to now in any case the meson fields (T{r) and 7r(r) have been restricted 
to the chiral circle a{r) + 7r(r) = /^, known as nonlinear realization of chiral 
symmetry (Weinberg 1968) . Due to this constraint the variation SEtot = is 
performed only with the chiral angle 9{r) as degree of freedom. The motivation to 
consider the action as an effective chiral theory of the nucleon from an instanton 
liquid model of QCD (Diakonov et al.l986) even suggests that only these degrees 
of freedom, namely the Goldstone ones, should be taken into account. On the 
other side in spirit of the Nambu-Jona-Lasinio-model there is no reason at all why 
the meson fields should be restricted to the chiral circle and the constraint eq. (3.6) 
seems to be an artificial one. Sieber et al.(1992) performed a detailed analysis of 
the full linear model with both a and yr-degrees of freedom and showed that at 
least in case of the proper time regularization scheme no localized solitonic solution 
with finite energy exists in this case but the system collapses into a configuration 
with both size and energy being zero while S = 1 is maintained. To this end they 
choose a parameterization for a{r) and n{r) in terms of their size R and their 
depths U by: 






(3.55) 



where /(r) is a strictly monotonously decreasing function with /(O) > 0, /(oo) = 
and 9{0) = —mi. Taking the limit w ^ oo, i? ^ oo but UR^ = const., the collapse 
described above occurs if 1 < a < |. This is due to the fact that in this limit for 
a > 1 altogether n valence orbitals cross the zero line and get part of the negative 
spectrum which means that they are not taken into account in the expression for 
the total energy any more (r]^;^; = in eqs. (3.41,3.42)) and the baryon number 
of the system gets (B) = {B)sea = n. Furthermore it turns out that the Dirac sea 
approaches its vacuum configuration, because any deviation from this is suppressed 
by the UV cutoff A. The reason for the occurrence of this behavior is the fact, 
that in solving the quark loop of NJL without vector mesons the sea energy Egea 
gets regularized by A whereas the finite baryon number {B) is not. On the other 
hand Em^s vanishes for i? — > oo as long as a < |, which proves the statement 
above. 
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Earlier works (Meissner Th et al.l989, Reinhardt and Wuensch 1989, Meissner 
Th and Goeke 1991), which claim the existence of solitonic solutions also in the 
linear model turned out to be premature in this respect since too small values of 
Kmax (cf- sect. 3.1.), where chosen in the basis for the diagonalization of the 
single particle hamiltonian h. This fact has been emphasized in the mean time by 
various authors (Watabe and Toki 1992, Sieber et al.l992, Kato et al.l993). 

In sect. 8. 4. it will be shown that a modified version (Meissner Th et al. 1993) of 
the Nambu-Jona-Lasinio-model , which in addition to the spontaneously broken 
chiral symmetry also simulates the anomalous breaking of scale invariance in QCD 
shows perfectly stable solitonic solutions also without the non linear constraint, 
while the properties of the soliton as well as the relevant nucleon observables stay 
nearly unchanged. Thus the restriction of a and n to the chiral circle is justified 
by a model which implements the concept of trace anomaly in QCD. Stability can 
also be achieved by adding the t/^(l) breaking t'Hooft term (Kato et al. 1993), 
at least for a certain range for t'Hooft coupling strength. 
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4. Restoration of Broken Symmetries - Collective Quantization 

As in nonrelativistic many particle physics we face the problem that the hedge- 
hog mean field solution of the classical equations of motion, which we have con- 
structed in the last chapter, breaks the rotational, isorotational and the transla- 
tional symmetry of the full theory. This means that the eigenstates | A) of the 
1-particle hamiltonian h do not carry good spin, isospin or momentum quan- 
tum numbers. Because we want to describe nucleonic systems, which have these 
quantum numbers, those symmetries must be restored. This is done by coupling 
the corresponding expectation value ((J), (T), (P)) through a Lagrange multi- 
plier to the effective chiral action S^ff like it was done in sect. 3. 2. in case of the 
baryon number B with the thermochemical potential /U. This turns out to be 
equivalent to consider the soliton in an (iso-)rotating or moving system, called 
cranking (sect. 4.1) or pushing method (sect. 4. 3), respectively. Because the grand 
spin G = J + T is a good quantum number of the hedgehog single particle state, 
spin and isospin degrees of freedom are coupled and it is enough to consider one 
of them, e.g.. the isospin. Whereas the spectrum of the momentum operator P is 
continuous, the isospin degrees of freedom have to be quantized. For this we will 
use the semiclassical collective quantization method (sect. 4. 2). Zero point energies 
appear, because the expectation values of the 2-particle operators T and P do 
not vanish (sect. 4. 4). Finally we will be able to write down the expressions for the 
masses of nucleon A^ and delta A in a system at rest and give numerical values 
(sect.4.4). 

In this chapter we will throughout consider the non-linear version of the 
Nambu-Jona-Lasinio-model , i.e. restrict a and tt to the chiral circle a + if = 
ay = f^ and use the notation: 

U = a + irif 

.. (4-1) 
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4.1. ISO-ROTATIONAL MOTION: CRANKING 

In order to restore the isorotational symmetry, i.e. to construct a system 
with good isospin quantum numbers (e.g. A^,A) we follow a procedure which was 
established in nonrelativistic many particle physics (see e.g. Ring and Schuck 
1980) and is known as cranking approach. It has been used also quite successfully 
in case of the Skyrme model (Adkins et al.l983) as well as in the chiral sigma 
model with valence quarks (Cohen and Broniowski 1986). Furthermore one can 
generalize the method to 5't/(3)-fiavor and calculate hyperon properties, which 
will be done in chapter 6. 

Adiabatic Isorotation 

The main idea is to perform an adiabatic isorotation of the hedgehog meson 
fields with the angular frequency O: 



a{x) -^ a{t, x) = a{x) 

U5{x) ^ t/5(t, x) = 7^t(t)t/5(x)7^(t) 
with the 5't/(2)-rotation matrix 



(4.2) 



nit) 



e* 2 ^ 



and its 5'0(3)-representation D, defined by: 



7^(t)r''7^t(t)=r^D^«(t) 



(4.3) 



(4.4) 



For simplicity we will first stay in Minkowski space without regularization and 
valence quarks, i.e. we set the thermochemical potential fu, to zero, and comment 
on the general case later on. 

The effective action in the rotating system is given by: 



JSeffitJ^] 






(4.5) 



with: q = W{t)q, and therefore: 

Seffimm = Spin 



iD + (3n- 



Spln 



idt 



h-n 



T 



(4.6) 
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From this we see, that O acts as a Lagrange multipher constraining the isospin 

T" = I d^xQ(x)0—Q(x) (4.7) 



d xq{x)/3—q{x) 



similar to the thermochemical potential n in case of the baryon number (sect. 3. 2). 
The 1-particle hamiltonian h{Q) in the rotating system can be read off from 
eq. (4.6) to: 



/i(0) = h-n- 



(4.8) 



Perturbative Cranking, Moment of Inertia 



We assume O to be small, so that the problem can be treated perturbatively. 
We will comment on the validity of this assumption later on. 

The 1st order in O vanishes 



because of hedgehog symmetry: 



(T«) = 



(4.9) 



E {X{G,G,)\UX{G,G,)) = 

A(G,G^) 

In second order we get: 

s^ffin) = TE{9) = T [emf + ^n'^e^^n^ 

where the moment of inertia tensor G" is defined by: 



(4.10) 



(4.11) 



e 



ab _ 1 s's,ff{n) 



1 



T dfi'^dn^ 1^=0 T 



Sp 






(4.12) 



In terms of the 1-particle eigenstates | A) (cf.eq. (3.5)) G'* is given as a sum over 
all particle {€iy > 0) - hole (e;^ < 0) matrixelements of the 'perturbation' r: 



Qab^Nc^ (A|r«|z.)(Hr''|A) 



(4.13) 



which is known as Inglis formula (Ring and Schuck 1980) 
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Of course it is necessary to regularize 5'g/-/-(0) and G. In case of the proper 
time regularization, which works in Euchdean space, the expression for G" was 
derived by Reinhardt (1989): 



^0 - -T 



1 °° 



7^e(eA,e^„s) 



^/^ (4.14) 



e se^ - e '^'^A exe^^^^ + e^e *'^a 

ej - 4 e^^ + ^A 



Hereby it is essential to perform the continuation of the rotation frequency O into 
the Euclidean space (O -^ Q.j^ = iQ) like it is done for the time component of a 
physical 4-vector, e.g. the time as itself (t ^ r = it) or the time component of 
a p- meson. This means that both O and O^; are real numbers and the fermion 
determinant in Euclidean space is: Spln(— zD — z/30^^). Otherwise it turns out 
that the regularized moment of inertia for the vacuum configuration t/5 = 1 does 
not vanish, which has to be considered as unphysical (Reinhardt 1989). 

The separation of the valence quark contribution by using the grand canon- 
ical effective action in the rotating system Sgff{p,Q) with the thermochemical 
potential p (cf. sect. 3.2.) is straightforward yielding: 



Qab 



^Nc ^ {val\r''\iy)\iy)T^\val) 



which is nothing but the standard 2nd order perturbation theory result for the 
orbital | val) with the perturbed hamiltonian eq. (4.8). 

On simple symmetry reasons G" is diagonal and can be written as G'* = 
5 0. By using the 1-particle eigenstates | A) and | z/) of the selfconsistent solitonic 
solution as they were obtained in sect. 3. 3, Goeke et al.(1991) and Wakamatsu and 
Yoshiki (1991) calculated numerically the moment of inertia, which will serve as 
the essential quantity for obtaining the masses of A^ and A from the soliton mean 
field energy E^p (see Fig. 4.1). 

Observables in the Rotating System 

Generally the expectation value of an observable O = J d xq^ {x)Oq{x) in the 
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rotating system reads (in Minkowski space): 

1 / VqVq J d^x[qf50q\e' / d^xq[iD{Ur>)]c 



(0)(0) 



T 



/ VqVqe^ / d^xq[~iD{U5)]q 



1 '^Q 1 






(4.16) 



}uj=0 



One should be aware that this expression differs from the simpler one 

1 6 



(0)(0) 



if O is an isovector by: 



Tdu 



Spin 






ba I r\b\ 



w^O 



(0°)(0)=D''°(0'')(0) 



(4.17) 



(4.18) 



In the perturbative cranking approach any (0)(0) is expanded up to 1st order 
in O. If this order vanishes generally only the 0th order is taken into account. For 
the case of (7^ and the isovector magnetic moment we will extend this method to 
higher orders and discuss the results shortly in sect. 5.1. and 5.2.. 

In the special case of the isospin operator eq. (4.7) we have: 



(T'^)(0) = L>^«(T^)(0) = D^^'&^VL^ = GD^^O^ 



(4.19) 



On the other side for the total spin: J'^ = J d^xq{x)(3 i^ + xx y) q{^) = G^—T°' 
one finds: 



(j«)(o) = (j«)(o) = (G"^)(o) - (T")(o) = -(T«)(o) = -eo« 



(4.20) 



which is true because [/i, G"] = and therefore {G"") vanishes up to 1st order in 
O as it is known from standard quantum mechanics perturbation theory. From 
eq. (4.19) and eq. (4.20) we read off the relation between spin and isospin expec- 
tation values: 

jjba^jh^ + (f ") = 



and furthermore 



{jy = (fy 



(4.21) 



(4.22) 
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4.2. Collective Quantization of the Iso- Rotational Degrees of Free- 
dom 

Up to now the expectation values (T)(0) and (J)(0) are continuous quantities. 
The semidassical collective quantization method, which was estabhshed by Adkins 
et al.(1983) in the Skyrme model, can be used for quantizing the rotational and 
isorotational degrees of freedom. The main idea is to express (T)(0), ( J)(0) as well 
as the rotational contribution of Sgff{Q), which can be written as J dtLroti^) = 
|T0 G, in terms of the components of the SU{2) matrix TZ{t): 

4 
n = b^ + ib^r" with ^ 6"^ _ 1 (4.23) 

and in terms of its time derivatives bQ,bf 



n^n 



t" 



(T«)(o) = eD^«(o)o^ = {-i)eTTr 

(j«)(o) = -0(0)0'^ = (-z)eTr^ [7^7^V"] (4.24) 

LrotiB) = eTr^[7t7tt] = 2e^[6"]2 

a 

and to consider 6" as space components on the 3-dimensional sphere S . The 
corresponding canonical conjugate momenta are given by: 

n'^(b^) = ^I^r^ = 4eb^ (4.25) 

The quantization is performed in the standard way by substituting the coordinates 
b'^ and the conjugate momenta tt*^ through the operators 

6"->6" 

on which the following commutator rules are imposed 

[6",7r^] = (+i)5"^ (4.27) 

The 6" and tt*^ are considered to act on a collective wave function | \I/(6)) which 
is normalized on the 3-dimensional unit sphere 5" with respect to the surface 
integration measure dn{b): 

i^m^ib)) = f d^{b){^{b)\b)m{b)) = 1 (4.28) 
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The dn{b) can be expressed by the angular representation of S^ by: 

dfx{b) = di/jdcj) sin (j)d9 sin 9 



(4.29) 



with 



60 = cos 6 

hi = sin 9 cos (j) 

62 = sin 9 sin (p cos '0 

63 = sin ^ sin (p sin •0 



(4.30) 



and 9 E [0, 27r], (/) G [0, tt], ^ G [0, tt]. The (T)(0) and ( J)(0) now go over into the 
isospin and spin operators t"' and j", respectively, which are given by: 



i"" 



'}k 



i.a_^ ,0_^ nbci^b 9 



^0_^ _ ^Z " Mmr.1 



db^ db^ 



9 ^klm^l 9 



d¥ 



(4.31) 



and fulfill the desired commutation relations: 






\j\f 



(4.32) 



The collective wave functions for proton (p) and neutron (n) with spin | and 
read: 

{b\ p t) = l(6l + z62) (5| ^ I) = 1(^0 _ ^^3) 



TV 



n 



(6|pi) = l(60 + z63) (6|ni)- ' '^^ ^^^' 



TT 



--{b^-lb^ 

Tl 



(4.33) 



whereas for the A we have e.g.: 

(6|A++ 



Z 11 



(4.34) 



from which the other A-states can be constructed by applying the spin- and isospin 
lowering operators j~ and t~ respectively. 

Finally by knowing those collective wave functions it is possible to calculate 
any matrix element of an observable acting in the collective space (6-space) be- 
tween two collective wave functions. The general recipe is the following: One uses 
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eq. (4.16) and expands the expressions up to 1st non vanishing order in the rota- 
tional frequency fl^. Then one replaces O" by the operator |v. This expression 
can now be sandwiched between the collective wave functions. Actually if terms 
higher than of first non-vanishing order in O are taken into account this simple 
recipe must be modified due to some time-ordering problems (see for details sect. 
5.1. and 5.2.). 

Important matrix elements, which we will need later on are e.g.: 

(pt |D^3|pT) = (-^)5^^ (4.35) 

and of course: 

(pT f\p^) = l (4.36) 

whereas for any collective state with quantum numbers T,J we have: 

{^Tjf\'^Tj)=T-{T + l) 

.o 4.37 

(*rjl/l*Tj) = ^-(^ + l) 

Therefore we see with eq. (4.22) that for the collective quantized eigenstates the 

general constraint 

J^T (4.38) 

holds. 

4.3. Translational Motion: Pushing 

As it does for the (iso-) rotational symmetry the soliton breaks also the trans- 
lational invariance of the effective action. Because we want to consider nucleonic 
systems with good linear momentum, especially for momentum zero, i.e. at rest, 
this symmetry has to be restored as well. The basic idea is the same as in the 
case of the isorotational motion: In order to avoid approximately the complicated 
boosting procedure one couples the corresponding operator, the linear momentum 
P = Y through a Lagrange multiplier v to the effective action SQff[U^{x)]. By 
applying the generators of the translation e*^ it can be seen to be equivalent to 
calculate the action for a field U^iv) = U^{x — vt) which moves with velocity v 
{pushing model): 

SeffP^m = Spln(iD[U5iv)]) = Spln|e--^*(ii?[C75])e*''^^j 

•'•' ^ ^ I J (4.39) 

= Spin (iD[U5])-vP) 
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Assuming that v is small, which we will also comment later on, we can treat the 
problem perturbatively: 

S,ff{v) = TE{v) = T (Emf + ^^'m* + 0{v^)) (4.40) 

Hereby we made use of the fact, that the 1st order in v, which is nothing but the 
expectation value of the linear momentum (P) in the static configuration Ur^{x), 
vanishes due to Ehrenfest's theorem already at the 1-particle level: 

{\\P\\) = \{\\[P\x]\>^)='-{\\[h\x]\\) = Q (4.41) 

and because of eqs. (3.16,3.17) also in general. The inertial mass M* defined by: 

1 ^'^S.fjpr,] 



M' 



n (4-42) 



is generally given by the particle-hole matrix elements of P similar to the Inglis 
formula eq. (4.13). Other than for the isorotational motion in case of the trans- 
lational motion Lorentz invariance of the NJL Lagrangian guarantees that the 
inertial mass M* is equal to the soliton mean field energy E]\^p: For small v, 
Sgff[U5{v)] in eq. (4.39) is the nonrelativistic Galilei limit of the effective action, 
which arises by boosting the soliton field in {i = ~ direction with the boost-velocity 
uj = arctanh(v) (cf. app. C). Indeed the soliton energy E transforms under this 
boost like the time component of a Lorentz 4- vector (Betz and Goldfiam 1983), 
i.e. the energy E{uj) in the boosted system is given by: 

E{uj) = cosh{uj)EMF = Emf + ^v'^Emf + 0{v^) (4.43) 

which after comparing with eq. (4.40) shows the desired identity: 

M* = Emf (4-44) 

An explicit proof of eq. (4.43) for the unregularized action Minkowski space using 
Poincare algebra as well as the selfconsistent mean field equation of motion (3.47) 
is given in appendix C. Pobylitsa et al.(1992) have shown, that this conclusion 
also holds, if one considers a regularized theory with a finite cutoff in Euclidean 
space as long as the regularization scheme is gauge-invariant. This means, in other 
words, that the central transformation identity 

Spin {i3t {u)AB{uj) } = Indet {i^t {cj)AB{uj) } = Indet^ = Spln^ (4.45) 

where A denotes an arbitrary operator and B{uj) the generator of the boost, has 
still to be valid, if Spin = Indet gets regularized. 
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4.4. Spurious Zero Point Energies - Masses of Nucleon A^ and A 

Spurious zero point energies for the translational as well as the (iso-) rotational 
motion, which have to be subtracted from the total energy, if a semiclassical quan- 
tization is performed, arise due to the fact, that even at mean field level the 
expectation values for the operators P and T are finite. 

Expectation Values of 1- and 2-Particle Operators 

Let us first have a look at the expectation values of 1- and 2-particle operators 
in general. For this we consider the pure 1-particle operators [O]^^)'!^ ](l) ^^^ 
the pure 2-particle operator [O ]('2-) defined by (Negele and Orland 1987): 

[0]{i) = J d^xql{x)Oaf3q(3ix) 

[0\^ = I d^xqUx)iOXpqp{x) (4.46) 

[0^]{2) = J d^xid^X2qi^{xi)ql^{x2)Oa^p^Oa^f3^qfi^{x2)qf3^{xi) 

where the a's and /9's denote some Dirac or isospin indices. Using the anticom- 
mutator relations for the fermion operators q and q' 

{qaix),qliy)} = S^{x-y)6c,f3 

{qa{S),qpm = (4.47) 

{qUx),qlm = 

we can decompose the squared 1- particle operator {[O]^-)} into: 

{[0]ii)}^ = [0\) + [0%) (4.48) 

Using the notation: 



the expectation value of the last summand in eq. (4.48) can be written as 
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{[O 



(zTy 



zR^H..o^^^^^«^'^^^'-" ^'-''^ 



a=0 



For static field configurations S{a) is proportional to T so that the last term in 
eq. (4.50) vanishes in the ground state (zero temperature limit) T -^ 00. Therefore 
we finally obtain: 

^lirn^ ({[0](i)}2) = {[0\)) + ([0](i))2 (4.51) 
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Zero Modes of the (Iso-)Rotational Motion, Band-Head-Energy 

From eq. (4.51) we obtain for the expectation value of the operator T in the 
rotating system (cf.eqs. (4.16,4.19)): 

{{[\l)?)m = {[T\i))m + (([f](i))(0))' = 0^02 + ([f'](,))(0) (4.52) 

The expectation value of the 1-particle operator [T ]n)(^) is independent of O, 
because [? ](i) commutes both with h and h{fl). Then the expectation value can 
be shown to get no contribution from the Dirac sea and hence simply reads: 

{[T]^l))m = {[T\i)) = Nc{B)l{l + 1) = iVc(S)„«4(^ + 1) = i (4-53) 

If we consider O as a Lagrange multiplier for adjusting the quantum number T, 
we have the condition: 

{{[f]ii)fm) = nT+i) (4.54) 

and therefore find: 

e2o2 = T(T + l)-([f2](l)) (4.55) 

Inserting this equation into the expression for the energy in the rotating system 
(cf.eq. (4.11)): 

E{n) = Emf + l'S>n'^ (4-56) 

one obtains for the energy of a system with isospin T and spin J: 

E -E ^I^l±A <[^'](1)^ r457^ 

^JT - ^MF H ^Q ^Q — (4.57) 

which means that the hand head term — 20 has to be subtracted from the 
cranked energy E^p -\ — 29 • The expression eq. (4.57) is familiar from nuclear 
many body physics, see e.g. Ring and Schuck (1980) and Blaizot and Ripka (1988). 

Zero Modes of the Translational Motion, Center-of-Mass Energy 

Instead of a proper boosting for the translational degrees of freedom we adopt 
a purely non-relativistic treatment, which consists in separating the whole motion 
into an collective as well as an intrinsic one (Ring and Schuck 1980). Thereby we 
hope that at least the magnitude of the effect is described reasonably. 
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For this let us look at a A^-particle system interacting by a purely local time- 
and velocity independent one-body force V{x). We consider the static and the 
pushed total Hamiltonians: 

^^^ J ^ (4.58) 

N P 2 



«-£2^+^(^>.-"^" 



respectively. The corresponding solutions of the time-dependent Schrodinger equa- 
tions 

at|vi>(t)) = if|vi>(t)) 

(4.59) 
dt\^v{t))=Hy\^y{t)) 

are connected by the unitary generators of the Galilei transformation (see e.g. Fonda 
and Ghirardi 1972): 

I ^„ (t) ) = ^^hMvh^^^Pvt^^MRv^ ^ ^^^ ^ (4_gQ^ 

with: 

N 

^ = E m 

N 
M = E ^(k) (4-61) 

^ 1 ^ 

^ = ]^ E ^(fc)%) 

Applying this transformation one can easily show that the expectation values of 
the operator ({[-P](i)} ) = J2ij'P(i)'P{j): in the moving and in the rest system are 
related by: 

(T.~P{^)~P{J)) (^) = (Ep(z)P(j)) + ^'^' (4-62) 

\ ij I \ ij I 

J2ij P{i)P{j) describes the center of mass motion and is also called the square of the 
collective momentum Il^^ii{v) = J2ij'P(i)P(j)- Using eq. (4.52) as well as the fact, 
that ([-P](i)) = in a static system (Ehrenfest's theorem eq. (4.41)) we obtain a 
relation, which fixes v : 

fililiv) = M^v^ + {[P\)) (4.63) 
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Substituting eq. (4.63) into the expression for the energy of the sohton in a moving 
system eq. (4.43) one recognizes that the total energy can be separated into a 
collective and an intrinsic part: 

([^'](l))' 



ZtjMF 



E 



MF 



2E 



MF 



(4.64) 



where we have set E^p for the total mass M. The collective term 2e°^^ i^ ^^^ 
kinetic center of mass energy and vanishes, if we transform to a system with resting 

/rp2i \ 

c.o.m. The relevant part is therefore the intrinsic energy Ej^p OE ^ ^^ which 

analogous to eq. (4.57) the spurious zero point energy ^p^ *■ ' has got subtracted 
from the classical mean field energy E^p. 

Summarizing we find finally for the total intrinsic energy for a system with 
spin J and isospin T whose structure is known from e.g. Ring and Schuck (1980) 
or Blaizot and Ripka (1988): 

E -E .I^I±A <[^'](l)) ([^'](l)) .4,,. 

General Aspects of the Zero Mode Treatment 

Aside from the fact, that we have handled the center of mass motion purely 
nonrelativistically, our treatment of the spurious zero modes goes clearly beyond 
the semiclassical approach, which can be seen e.g. from the condition eq. (4.55) for 
the Lagrange multiplier O , which is incompatible with the one for O" in eq. (4.19) 
on a classical level since eq. (4.55) allows imaginary omega. One should also note 
that the 'correction' terms to E^p, which as itself is of 0{Nc) are of different 
order in Nc^ namely 0{j^) for the centrifugal term 20 ^^^ ^(-^c) ^^ ^^'^ 

band-head term <!^ aid the center of mas. term <!^. respectively. 

In nonrelativistic many particle physics the form of eq. (4.65) can be obtained 
within certain approximation by using Peierls-Yoccoz projection techniques (Ring 
and Schuck 1980) for the angular and linear momentum. Those techniques have 
been successfully applied in soliton models with valence quarks (Birse and Banerjee 
1984, Birse 1885, Liibeck et al.l986, Fiolhais et al.l988, Fiolhais et al. 1991, 
Neuber and Goeke 1992, for a review cf. Birse 1990). Unfortunately it is up to now 
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not clear how to establish them in case of the NJL soliton, because the regularized 
Dirac sea does not allow the definition of a Fock state, which is necessary for 
applying projection methods. 

In case of the Skyrme model it turned out that considering RPA fluctuations 
around the mean field solution, which breaks the (iso-) rotational and translational 
symmetry of the full theory, allows a treatment of the corresponding eigen modes 
(Moussallam and Kalafatis 1991, Holzwarth 1992), analogous to the case of two 
dimensional soliton models (Rajaraman 1982). Altogether these modes lower the 
classical soliton energy and about 80% of the lowering originates from the rota- 
tional and translational zero modes. Hence it is probably a good approximation 
in (4.65) to concentrate on the zero-modes of the baryon. 

A full treatment of the RPA-modes in the present NJL model has not been 
done by now since the Dirac sea complicates the formalism tremendously. On the 
other hand, as long as there exists no consistent treatment of the zero modes in our 
model, the method described above is as good as the semiclassical quantization 
by itself and one can hope that eq. (4.65) describes at least roughly the mass of a 
particle at rest ((11^^;/) = 0) with spin J and isospin T. 

Masses of A^ and A, Numerical Results 

From eq. (4.65) we can read off the expressions for the total masses of a nucleon 
A (J = T = ^) as well as a A (J = T = ^), respectively (Pobylitsa et al.l992): 

89 86 2Emf (4gg) 

Af A - -Bmf + ge - ge - ^^^^ 

For the numerical calculation we take the mean field energy of the selfconsistent 
solutions (cf. sect. 3.3). The G as well as the expectation value ([-P ](i)) can be 
obtained from the 1-particle eigenstates | A) of this solution using eqs. (4.15, 3.17). 
The numerical results in dependence of the constituent quark mass M are shown 
in Fig. 4.1, 4.2 (Pobylitsa et al.l992). In the relevant region M ^ 400MeV it turns 
out that the rotational zero point energy -^ lies around lOOMeV whereas the 
translational zero point energy amounts to about 300MeV. This is the order of 
magnitude obtained also in nonrelativistic quark models (see e.g.Bhaduri 1988 and 
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ref.therein) as well as in relativistic soliton models using Peierls-Yoccoz projection 
techniques (Fiolhais et al.l991, Neuber and Goeke 1992). On the other hand 
especially the value for the translational zero point is quite high in comparison 
with the soliton mean field energy Ej^p ^ 1200MeV, so that the perturbative 
treatment in O and v appears more than questionable. It is also interesting to 
note that the 'cranking' term 29 can be numerically of the same order of 

magnitude as the band-head energy — 29 ■, although they are of different order 
in iVc. 

Finally M^ and M^ (Fig. 4.2) are oX M ^ 400MeV with M^ ^ 900MeV and 
M/^ ^ llOOMeV, respectively, somewhat close to their experimental values (M^v = 
938MeV, M/^ = 1230MeV). Especially one notes that the nucleon becomes stable 
against decay into A^c = 3 free quarks as consequence of the subtraction of the 
spurious zero point energies. 



62 



5. Nucleon Form Factors and Observables 

We are now supplied with the basic techniques to calculate various nucleon 
observables from the mean field solutions obtained in chapter 3. We will consider 
especially the electromagnetic form factors of proton and neutron including radii 
and magnetic moments (sect. 5.1) as well as the axial vector coupling constant (7^ 
and the axial form factor g^iQ ) (sect. 5. 2). Furthermore the pion nucleon coupling 
QttNN 3.nd the form factor g-n-NNiQ ) ^.re discussed and the Goldberger-Treiman 
relation between (7^ and Qt^jsin is exhibited (sect. 5. 3). Finally we deal with the 
nucleon E-term (sect. 5. 4). The spin of the proton that is carried by the quarks 
will be discussed in the context of the SU{?>) soliton in chapter 6. 

5.1. Electromagnetic Form Factors, Moments and Radii 

Electromagnetic Current 

On a classical level the electromagnetic current jfl^ arises as Noether current 
of the electromagnetic U{1) phase transformation: 

q -^ e'^'^q (5.1) 



with the charge matrix 



and is therefore given by: 



«^(' 4)4^ + T (^•^' 



fr = 9>Q9 (6.3) 



It is the sum of an isoscalar and an isovector part. Applying the general method 
described in sect. 4. 2. we find for the sea part of the expectation value of if[^ in 
the rotating (cranked) soliton in Minkowski space (cf.eq. (4.16)): 



^■""(^» - ^iii^^p'" 



iD + (3— + A^(x)7^(/37^g)7^t 



(5.4) 



If one performs a Wick rotation to Euclidean space the time component Aq has 
to be handled in the same way as the cranking frequency O (cf. sect. 4. 2), i.e. the 
time component of a physical 4- vector: Aq -^ A/^ = —iAq, where A4 and Aq are 



63 



both real. The valence part arises in the standard way by introducing a chemical 
potential (cf. sect.3.2). 

One should mention that d^^ (jf[^{x)) does not necessarily vanish on the 
fermion 1-loop level, because of the regularization involved. The same is true also 
for the divergence of the axial current (cf.eq. (2.7)). Apparently the result of the 
divergence depends on the regularization scheme used and indeed, for momentum 
cutoff schemes the Noether theorems are violated in the solitonic sector (Doring et 
al.l992). On the other hand the proper-time method preserves those invariances. 
This has been checked for t/(l), SU{Nj)y and SU{Nf)j^ transformations. 

Electric Form Factors of Proton and Neutron 

The electric form factor Gj^{q ) is defined by the matrixelement of the time 
component Jq^ between nucleon states | Nj{pi, Si.ti)) with 4- momentum pi, spin 
z-component si and isospin z-component tj (i=l,2) (Bernstein 1968) 

(iV2(p2, S2, t2) IjT(0)I iVl(Pl, SiM)) = u{p2. S2, t2)(3Qu{pi, Si, h)GE{q^) (5.5) 
In eq. (5.5) q = P2 ~ Pi denotes the 4-momentum transfer and 



the spinor of a free pointlike nucleon with mass Mj^. It is convenient to use the 
Breit system defined by: 

pi = {E,-q/2) p2 = {E,+q/2) q={0,q) 

. (5.7) 

E = ^Ml + cf/A Q^ = -q^ = i Q = \q\ 

There eq. (5.5) simplifies to: 

^.N2{P2.S2M)\3TmNl{VX.SiM))=GE{q^)x\^x\^QXs^Xt^ (5.8) 

The nucleon state | Ni{j>i, ^ijii)) is treated in the static approximation and center 
of mass corrections are neglected. This means: 

— * — * 

I N{^, s„ t,)) = \ iV(|, s„ t,)) = v/(27r3)5(0)| s,t,) (5.9) 

In our approach we take for | Sjtj) the collective quantized nucleon states \p ^), 

\p I), \n^) or \n |). 
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Like the electromagnetic current jfl"' itself G]^{q ) separates into an isoscalar 
and an isovector part 

Gsiq^) = Ige^-^Q^) + kGE^-^Hq^) (5.10) 

It is easy to see that for the isoscalar part the first nonvanishing term is of zeroth 
order in O, whereas for the isovector part is the term of first order in O. These 
two quantities are then related to proton and neutron experimental currents by 
GE^=^{q^) = GE{q\ + GE{q% and GE^^Hq^) = GE{q\ - GE{q%. 

From eq. (5.5) one can extract Ge {q ) and Ge {q ) after some algebra 
to (Gorskiet al.l992): 

Ge^=%^) =ld^xe-''^''^l\<P,alix)\^ - ^EsigneAl<^A(^)l 



GE'-Hq')= d'xe-^^^^^-^ d^l-Y. 



exe """x + ei.e~*^^ 1 e '"'^a - e"'*'^^ 



ds .,-^ 

■[<p{{x)TMm<PUy)r<Pxm] 

(5.11) 
where denotes the moment of inertia (cf. sect. 4.1) and (pxi^) = (^1-^) the 1- 
particle eigenfunctions in ^-representation. 

The results for the proton and neutron form factor as well as the corresponding 
charge densities are shown in Fig. 5.2-5.5 (Wakamatsu 1991, Gorski et al.l992). 
As one can see, the proton form factor G^{q ) is described very well in contrast 
to the neutron form factor G^{q ). Furthermore we recognize that the neutron 
charge density is dominated by the sea quarks at large distances r, which confirms 
the popular picture that the long negative tail of the neutron charge distribution is 
made by the pion cloud (Thomas 1983), which is connected in our model with the 
polarization of the Dirac sea by gradient or heat kernel expansion of the fermion 
determinant. 

The corresponding quadratic radii 



?2\ 



JG_E{f_ 



{R1e = -^ T, ' \a2=0 (5.12) 
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are shown in Fig. 5.1. The fact that the < B? >^ is negative is due to the long 
negative tail in the corresponding charge distribution. It originates in the present 
model from the polarization of the Dirac sea, which in other models corresponds 
to the pion cloud. This negative tail is obviously overstressed in the present model 
< R >^= — 0.21fm compared to the experimental value of < R >^= — 0.12fm . 
One should compare these values with the Skyrme model, which gives — 0.36fm in 
the scalar version (Braaten et al. 1986a, 1986b) and — 0.24fm with vector mesons 
(Kaiser et al.l987). 
Magnetic Moments and Form Factors 

The same considerations can be applied in case of the magnetic form factor 
which is related to the space components of the electromagnetic current jf^ in 
the Breit frame through 

(iV2(p2,«2,t2)|?"^(0)|iVi(pi,.i,ti))=^^^xt^xI^(aXQ-)X,s,X,, (5.13) 

The main difference compared to G]T;{q ) is the fact, that in case of GM{q ) the 
leading order of the isovector part is (9(0 ) whereas the leading order of the 
isoscalar part is 0(0 ). 

At g = one obtains the magnetic moments: 

which have been calculated in the present model by Wakamatsu and Yoshiki 
(1990). The results are shown in Tab. 5.1 and compared with the experimental val- 
ues {^P = 2.79, ;u" = -1.91, fi^^'^ = 4.70 and ^'^^^ = 0.88). One recognizes that 
for the relevant values of the constituent quark mass M ^ 400MeV the isovector 
part /U has a relatively large sea quark contribution and its total value comes 
out too small. The isoscalar part fi is clearly dominated by the valence quarks 
and can be reproduced quite reasonable. 

A calculation of the full q dependence of G^[{q ) has been performed recently 
(Gorski et al.l993) and shows a good agreement with the experimental results as 
far as the g dependence is concerned (Fig. 5.6,5.7). 

Similiar to the axial vector coupling constant, to be discussed in sect. 5.2, 
actually the first non- vanishing term for the isovector magnetic moment is of zeroth 
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order in the rotational frequency O. However there are important corrections in 
hnear order of O, which are entirely due to the time-ordering of operators in 
non-local theories (Christov et al. 1993b). Without going into mathematical 
details the final results ~ (O , O ) are presented in Fig. 5.8 and compared with 
the experimental data (Christov et al. 1993b). Apparantly the corrections ~ O 
are quite important (Compare with the axial vector coupling constant in Fig. 5.9). 

Furthermore one has to keep in mind that all calculations have been per- 
formed in zeroth order of the pushing velocity v, which means, that any center of 
mass corrections to nucleon observables other than the nucleon mass (cf. sect. 4. 4) 
have been completely neglected. On the other hand one knows from nonrelativis- 
tic quark models (Bhaduri 1988 and ref.therein) as well as other effective quark 
models (Betz and Goldfiam 1983, Dethier et al.l983, Fiebig and Hadjimichael 
1984a, 1984b, Braaten et al. 1986a, 1986b, Luebeck et al.l986. Leech and Birse 1989, 
Stern and Clement 1989, Fiolhais et al.l991, Neuber and Goeke 1992) that those 
corrections might be noticeable. 

5.2. Axial Vector Coupling Constant and Axial Form Factor 

The general form of the matrix element of the axial current (cf.eq. (2.7)) 

A«(x) = ^(a:)7^75r^(a:) (5.15) 

between nucleon states (eq. (5.9)) can be written as: 

(iV2(p2,S2,t2)|A«(0)|Ari(pi,Si,ti)) 

^a (5.16) 

where u{pi,Si,ti) , z = 1,2 denotes the free nucleon spinor (cf.eq. (5.6)). gAio. ) 
is the axial and /i^(q' ) the polar form factor of the nucleon. At zero momentum 
transfer QAiQ ) reduces to the axial vector coupling constant qa = 9AiQ = 0) 
given by: 

9A-2{p^ ||rf3x(ip|(x))(0)|pT) (5.17) 

The Axial Vector Current in Zeroth Rotational Order 

Using the method of collective quantization described in sects. 4.1 and 4.2 we 
have in zeroth order O (cf.eqs. (4.16,4.35)): 

g^ = '^Ijd^'x (q\x)^q{x)) (5.18) 
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Furthermore the sea contribution of (7^ gets regularized due to eqs. (3.17,3.18) 
with O = ^. 

Using the mean field equations of motion eqs. (3.47, 3.51) it can be shown that 
the PCAC relation eq. (2.8) is maintained for the expectation value of the axial 
current, i.e. it is valid also at the 1-quark loop level to order ~ Vt^ > as long as the 
regularization scheme applied respects chiral invariance (Meissner Th and Goeke 
1991). 

^^'{Al{x)) = -mlUTvyx) (5.19) 

As we have mentioned in the last section this is true, if one uses e.g. the proper 
time regularization, but violated in case of momentum cutoffs. 

Fig. 5.9 shows the numerical results for various constituent masses M cal- 
culated from eq. (5.18) with the selfconsistent profiles. In the relevant region 
M ^ 400MeV one obtains in zeroth order of O an axial vector coupling constant 
9 A ~ O-S? which is too small compared to the experimental value of (7^ = 1.23 
(Wakamatsu and Yoshiki 1990, Meissner Th and Goeke 1991). The valence quark 
contribution is clearly dominating. 

Applying eq. (5.19) to eq. (5.17) (7^ including the valence part can be written 
as radial integral over the pion field: 

9 A = -Y'^lfiT I drr^Tvir) (5.20) 

This sum rule turns out to be fulfilled numerically as well (Meissner Th and Goeke 
1991a). 

The momentum dependence of the form factor QAiQ ) is usually parameterized 
by a dipole form: 

SAi^')-fl.I_)-' (5.21) 



9a{0) V Ml. 

where the most recent experiments determine the dipole mass to M^ = 1.09 ± 
O.OSGeV (Ahrens et al.l988). Meissner Th and Goeke (1991a) obtained in the 
present model a value of M^ = 0.95GeV, which is almost compatible with the ex- 
perimental result. Hence the q dependence of the axial form factor is reproduced. 

The axial vector current with first-order rotational corrections 
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Recently it has been shown, that due to the non-commutativity of some collec- 
tive operators (Wakamatsu and Watabe 1993) and due to the explicit time-ordering 
of time-dependent operators (Christov et al. 1993b, Blotz et al. 1993c), that there 
are important contributions from the rotational corrections in the first order of 
the frequency O ~ l/Nc- Their origin lies in the fact, that the generators of the 
SU (2) gpij^-group J^ do not commute with the D- functions D . Furthermore one 
has to take care of the time-dependence of the operators, because one has to refer 
to the operator formalism (and not to the functional integral) when dealing with 
the rotation matrices 'R-{t) (Christov et al. 1993b). Therefore in the expansion of 
(4.16), there emerges a term of the form ~ [J^^, -D ], which is non- vanishing. 

Actual calculations (For details see Christov et al. 1993b) including these cor- 
rections improve the value for (7^, such that the total value of (7^ for M ~ 400Mey 
in NJL comes out to be (7^ = 1.15 ^ 1.30, which is quite close to experiment. The 
formulas for these additional terms are presented together with the corresponding 
terms in SU(3) in chap. 6. Detailed values for g/^ in SU(2) including these cor- 
rections are presented in Fig. 5.9 and the importance of the O corrections to the 
isovector magnetic moment has been shown in Fig. 5.8. 

There are two philosophies (Blotz et al. 1994). First: One evaluates the pion 
and sigma field from the equation of motion eq. (3.51), which are given in (9(0 ). 
Then one evaluates with this the axial vector current and a new pion field, both 
in order (9(0 + O ). This treatment is strictly consistent with an expansion 
in the 1/Nc spirit as far as rotations are concerned and it is consistent with the 
perturbative cranking procedure. On the other hand PCAC is violated. If one uses 
the numbers of Alkofer and Weigel (1993), one may obtain a rough feeling how far 
PCAC is violated and for reasonable values of the constituent mass this violation 
seems to be small. Second: One derives equations of motion for the pion field and 
also for the axial vector current, both up to first order in O. In this treatment 
PCAC is strictly conserved. However the treatment does not comply with a strict 
1/Nc expansion and the perturbative expansion of the cranking procedure. On 
the present level of investigation one cannot prefer one procedure to the other. 
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5.3. PioN NucLEON Coupling and Form Factor 

The pion nucleon form factor QT^NNiQ ) is defined through the nucleon matrix 
element of the pionic current: J^{x) 

(D + m2)7r"(x) = J^{x) (5.22) 



(5.23) 



as: 

{N2{p2.S2.t2)\4iO)\Ni{pi,Si,ti)) 

= {-i)u{p2, S2, t2)r''-f5u{pi, si, ti)g^^^{q^) 

For static pion fields one has: J^{r) = (—V +m^)7r"(f), which remains also valid 
in the cranking approximation because the 2nd time derivative in eq. (5.22) is of 
2nd order in the cranking frequency O and therefore neglected. From eq. (5.23) 
by using eq. (5.22) the g-j^NN^Q) (Q = ^Q ) can be finally be extracted to: 

Especially at Q = —q = one has: 



2 47r /■ o 
— mj^— / drr Ti[r) (5.25) 

9 J 



2Mn 

Comparing eq. (5.25) with eq. (5.20) we recognize the famous Goldberg er-Treiman 
relation (Goldberger and Treiman 1958, Cheng and Li 1984): 

9nNNiO) = ^9A (5.26) 

Jtt 

which arises as a consequence of the PCAC relation eq. (5.19) and holds as long 
as the regularization scheme respects chiral symmetry. 

The above formula deals with spacelike q = —Q < 0. However, direct 
experimental data for the ir — N interaction are only available for timelike g , even 
on shell q = —Q > m^. The eq. (5.24) can be analytically continued to the full 
timelike region q = —Q > (Cohen 1986): 

9ttNn{Q^ < 0) r / 2^^2x 47r ri? 3( h{Qr) \ , V 

— 2M-^ — =H--+^)-Yyo ^^^ i^^j"^'^ 

+ [-{ml + Q') ■ '-^Ae-^^^R (^)] (5.27) 



+ 



.^A^-m^R_L (^^ smh{\Q\R) + \Q\ coshdQIi?))' 
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The A is given by the asymptotic form of the pion field in eq. (3.53). Especially 
for on shell pions q = —Q = m^ one finds 

9ttNn{q'^ = ml) ^ Z^A (5 28) 

with the experimental value g-n-NNiQ = "^n-) — 13.6. Tab. 5.2 shows the values for 
9t:Nn{h = ''Thr) ^^^ dirNNi^) ^^ ^^^ model, where the nucleon mass from eq. (4.66) 
and Fig. 4.2 have been used. Because of the Goldberger-Treiman relation eq. (5.26) 
it is clear that the discrepancy between the theoretical and the experimental value 
of QtyNN = 13.6 is the same as in case of gy\ evaluated to fl^ > considered in the 
last section. 

In the spacelike region the q dependence of the pion nucleon form factor 
9-kNn{h < 0) can be fitted to A^A^ scattering data by using one boson exchange 
potentials (OBEP). In a monopole parameterization one obtains in such an ap- 
proach 

9ttNn{q'^ < 0) ^ ^ttNN - ml 
9ttNn{q'^ = ml) At,nn - q'^ 

with a monopole mass of A^^ya^ = 1530MeV (Machleidt et al.l987). In our model 
we obtain a monopole cutoff of A^^jYiV — 790MeV. Also all other chiral models 
predict a much lower value than the OBEP, e.g. the Skyrnie model without vector 
mesons 580MeV (Cohen 1986), the Skyrme model with vector mesons 850MeV 
(Kaiser et al.l987), the projected linear chiral sigma model 690MeV (Alberto et 
al.l988). The same order of magnitude for At^tva^ is needed in charge exchange 
reactions (Esbensen and Lee 1985), in more recent estimates within meson ex- 
change models (Deister et al.l990, Janssen et al. 1993). A similiar value of A is 
also obtained by substituting the experimental values for g-j^^^iq = ml) = 13.6 
and g-n-NNi^) = "T'gA = 12-4 from the Goldberger-Treiman relation (5.26) into 
eq. (5.29) which then yields A-j^nn = 468MeV. The discrepancy between the 
OBEP calculation and all the other approaches might lie in the fact that for a cor- 
rect description of the NN interaction other, more complicated processes than the 
TT AA-vertex have to be taken into account. A refitting of the AA-phase shifts by 
means of OBEP models with different exchange processes yields also lower values 
for A^TViV (Holinde 1992). 
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5.4. NucLEON Sigma Term and Form Factor 

The nucleon sigma term E^y is defined as analogon to the quark condensate 
{qq)v (2.22) in the nucleon sector 

T.N = mQ [ d^x{N \q{x)q{x) \ N) (5.30) 



As an isoscalar quantity the E^y stays uninfiuenced by the cranking procedure in 
SU(2). Because the current mass niQ couples to the Dirac operator like a Lagrange 
multiplier for the condition niQqq (cf.eq. (3.14)) one can write E^y in the convenient 
form: 

^^ = ^° dm rh^O (^-31) 

The corresponding form factor is defined by: 

(^{q )u{P2,S2,t2)u{pi,si,ti) = mo{N2{p2, S2,t2) \q{0)q{0)\ Ni{pi, si,ti)) (5.32) 

which gives with eqs. (5.6,5.9): 

^(Q^ = V) = j d\jQ{Qr){q{r)q{T)) (5.33) 

At the Cheng-Dashen point q = 2?ti^, the cr{q ) can be determined from tt — A^- 
scattering data to a(2m^) ~ 60MeV (Gasser and Leutwyler 1982). From dis- 
persion relations an extrapolation to q = has been performed, which gives 
A(j = cr(2m^) — (t(0) = 15MeV (Gasser and Leutwyler 1991). Hence the present 
NJL model should evaluate the sigma term E^y of eq. (5.30) to 45 ± 5MeV. 
One should note however that from chiral perturbation theory a much smaller 
value A(j ~ 5MeV is obtained (Gasser et al.l988a,1988b). In contrast to all 
other observables in our model E^ turns out to be somewhat sensitive to the 
regularization scheme applied. For the proper time scheme one gets a value of 
E^ = (t(0) ^ 35...40MeV (Meissner Th and Goeke 1991, Wakamatsu 1993), 
which is also obtained in the Pauli-Villars scheme (Schueren et al.l992), whereas 
with an extension of the proper time method 



^^^^drF{T)^ J^ <P{T)FiT) (5.34) 

also higher values for E ^y can be obtained for an appropriate chosen cutoff function 
4>{t) (Blotz et al. 1993b). Details for SU(3) will be given in sect. 6.1.. 
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The eq. (5.33) can be analytically continued to the time like region cp' = —Q > 
0. The scalar form factor has been calculated by Schueren et al. (1992). One finds 
from the deviation from the Cheng-Dashen point (Schueren et al.l992) 

A^ = a{2ml) - a{0) = 7MeV (5.35) 

similar to the value obtained from chiral perturbation theory (Ajj ^ 5MeV, Gasser 
et al. 1988a, 1988b) but only half as large as the one by means of dispersion 
relations (A^ ^ 15MeV, Gasser et al. 1991) and by a Bethe-Salpeter approach in 
the meson exchange picture (Pearce et al. 1992). 

Wakamatsu (1992b) has also calculated the isospin violation of the qq content 
in the nucleon (Gottfried sum rule) and found a reasonable agreement with the 
experimental result from NMC data. 
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6. The SU(3)-flavour NJL-model 

In the foUowing chapter we wiU investigate the extension of the SU(2)-NJL 
model to the larger symmetry group of SU(3) (Blotz et al. 1992, 1993b, 1993c, 1993d; 
Weigel et al. 1992a, 1992b). In the vacuum sector, we extend in this way the 
number of Goldstone bosons to the kaons and the eta. Especially these lightest 
mesons are often considered as the dominating degrees of freedom for the low 
energy regime of the strong interactions. Because the masses of these particles 
are still much lower than the nucleon mass, the extension to this larger symmetry 
group is a priori sensible. 

Historically the success of hadronic isospin symmetry SU(2) (Heisenberg 1932), 
which is based on the almost mass-equality between the up and down quarks and 
the flavour independence of QCD, leads to the discovery of flavour SU(3) (Gell- 
Mann and Ne'eman 1964). It was found that, assuming that u,d and s quarks 
transform as the fundamental representation of SU(3), the spin 1/2 nucleons be- 
long to the 8-dimensional representation and the spin 3/2 particles to the 10- 
dimensional representation of the group. From the experimental observation of 
these symmetrically arranged particles in the multiplets one is usually inclined 
to consider SU(3) as the more 'physical' symmetry. That this can be supported 
within a selfconsistent chiral model for mesons and baryons, is shown below. After 
considering the modiflcations of the vacuum sector due to SU(3), which include 
the mixing of the ry-ry'-system, the t/^(l)-anomaly and the Gell-Mann Okubo mass 
relation, the baryon number one sector is described. 

To this aim we will concentrate here on the quantization of rotational zero 
modes (Adkins et al. 1983), which rests on the assumption that SU(3) is indeed a 
good symmetry of the strong interaction. An alternative treatment considers the 
baryon as a bound state of a heavy meson and the background fleld of the SU(2) 
soliton (Callan and Klebanov 1985,1988; Callan et al. 1988; Weigel et al. 1993). 
It is believed that this treatment gives an exact answer for extreme heavy mesons, 
such as those from a SU(4) representation, but it is still not clear whether it is 
more appropriate for the intermediate energy scale of the kaon system of SU(3). 

Therefore special attention is drawn here to the quantization of the generators 
of the full SU(3) group and to second order corrections in the strange current quark 
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mass for the splitting of the multiplets. Finally the axial vector coupling constants 
9% 9 A ~ 3 A ^^^ 9a- currently under refined experimental investigation, are pre- 
sented with some recently found corrections arising from some non-commutativity 
of generators (Wakamatsu and Watabe 1993) and explicit time-reordering (Chris- 
tov et al. 1993, Blotz et al. 1993d). 

6.1. The NJL with SU(3)-flavours: The vacuum sector 

The vacuum or mesonic sector of the SU(3)-NJL model is intensively discussed 
in the review articles of Klevansky (1992) and Vogl and Weise (1991). However 
they mainly used an operator formalism based on the pure quark NJL Lagrangian. 
So we will summarize here the most important parts within the functional formal- 
ism. This has advantages in the baryonic sector, which is our main concern. 

SU(3) Symmetry Breaking and the Redundant t/^(l) 

The Nambu-Jona-Lasinio Lagrangian with scalar and pseudoscalar couplings 
on the level of the four-fermion interaction and with a SU(3) symmetry is conve- 
niently written as 

Cnjl = qix){ip - m)q{x) - ^ [{q{x)X\{x)f + {q{x)i^5>^\{x)f] (6.1) 

where m = diag(mu, m,;, m^) = mil + m2A3 -|- m^X^ is the current quark mass 
matrix. The Aa, a = 0, . . . , 8 are the usual Gell-Mann matrices with A = a/(^)1. 
Under infinitesimal chiral SU(3) transformations the quark fields transform as 

g-(l-z^AV-z^A«/3^5)? (6.2) 

Q--g-(l + ziAV-zlA«r75) (6.3) 

From the eqs. (6.2, 6.3) it is clear that the singlets {qX q) and {qi'j^X q) cor- 
responding to U(3) chiral meson fields have to be included in the Lagrangian in 
order to be invariant under the chiral SU{3)ji (E) SU{3)]^ transformation. 

This is in contrast to the SU(2) case, where we had the freedom to choose the 
chiral fields from a SU(2) or U(2) representation. In addition, this classical La- 
grangian (6.1) is invariant under t/^(l) 0^/^(1) transformations, where the singlet 
and octet parts of the quark bilinears qX'^q and qi'^^X^'q transform independently. 
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Currents and Divergences - PCAC 

The Noether currents of the classical Lagrangian from the infinitesimal trans- 
formations (6.2, 6.3) are given by 



dC 1 

"^^ = "m:^)^''' = -^-7m752 A'^g (6.5) 

and their divergences from the Gell-Mann Levi equations are given by 

^^V^ = ^=q^l[X-,m]q (6.6) 



d£ 1 

9M^^ = ^ = ?-^752{A^mh (6.7) 

Assuming isospin symmetry, i.e. niu = rrif^ = rh, one observes that 

a^lAi = 0, for a = 1,2,3,0,8 (6.8) 

from which follows consistently the conservation of isospin and hypercharge. Fur- 
thermore the divergences of the axial currents are proportional to the correspond- 
ing pseudoscalar quark bilinears, which will be identified later with the correspond- 
ing composite meson fields. We have 

a^A«=l'2'3 = -qt^^X^qm (6.9) 



a^A«-4,5,6,7 ^ _qi^^x\^{m + m.s) (6.10) 

In matrix form the remaining two currents are given by 






9,{''.]=^{\] (6-11) 



A*- 



with 
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and v^ = —qi'j^X'^q- So we observe that the singlet and octet currents are mixed. 
However, we can disentangle these terms by means of an orthogonal transformation 
according to 



d, xn=^LJ (6-13) 



with 



Here V — IZVIZ is a diagonal matrix and the rotation matrix is given by 

/ cos 9 — sin 6 \ 
nO) = ( . „ „ . (6.15) 

\ smt' cose/ / 

The result of this diagonalization is an angle ^, tan2^ = —2\pl i.e. d ~ —35.26, 
which is called the ideal mixing angle. We wilU come back to this point, when we 
consider the mesonic mass spectrum. 

The eigenvalues of T> are simply m and m^. This can be understood in the 
following way. Whereas the former states, which we denote by ?7o and ryg are 
group theoretically given by ?7o ~ ww + rfrf + ss and ryg ^^ uu + dd — 2ss, the rotated 
states are given hj fJQ = rj ^^ uu + dd or t)^ = r]' ^ ss. This is reflected by the 
two eigenvalues of "D, which are either strange or non-strange masses. This will 
become clearer when we consider the mesonic two point functions. 

Ujjyl) Symmetry Breaking 

As a consequence of the 5't/(3)^ Cg> SU{?>)i invariance in the chiral limit and 
assuming a spontaneously broken vacuum by a non-vanishing expectation value 
of some scalar composite fleld, we immediatedly get from Goldstone's theorem 
that according to the 9 broken generators of the symmetry group there will be 
9 massless Goldstone bosons. In nature however, the tj' with rrirji = 957MeV 
cannot be regarded as a Goldstone boson, because its mass is similiar to two twice 
the constituent quark mass, 't Hooft (1976b) proposed a mechanism supported 
in QCD, which breaks the redundant Uj\{l) symmetry, but conserves the SU(3) 
symmetry. Such a term is given by 



^det = ^ 



det(qiPRqj)+det(qiPLqj)l (6.16) 
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where i = u.d^s and -Pr/x = ^(1 ± 75) are the right and left heUcity projection 
operators ^^. Under an infinitesimal chiral transformation it transforms into 

Cdet = «{det(l - iA^*/?^)det(qiPRqj) + det(l + iA^*r)det(qiPLqj)} (6.17) 

This expression is obviously invariant if and only if the matrices \a are traceless 
and therefore det(l + iA^*/?^) = 1. This is the case if the transformation corre- 
sponding to /3o is excluded from the full axial group t/4(3). So one is left with a 
SU{2))y ® SU{?>)j^ ® U{l)y symmetric expression. 

Bosonization 

In order to bosonize expressions like C^^^i, which are not only quadratic in 
the fields but also of higher order, one has to modify the Gaussian bosonization 
procedure of sect. 2.1. Therefore we introduce according to Zahed and Brown 
(1986), Reinhardt and Alkofer (1988) 

1 = JvSaVPad{Sa - qh\)d{Pa - ^^^sA"?) 



= I VSaVPa I VaaViVae^pi J d'^xaaiSa - q^X^'q) + Tia{Pa " ^^HsA^) (6.18) 

In this way, any quark bilinear in £^g^ can be replaced by the fields Sa.Pa- A 
stationary phase approximation for Sa^ Pa and a^, tIq immediatedly relate Sa^ Pa 
with a a and tTq, according to 

^^ AGSa + aa + ^ = Q (6.19) 



dSa dSa 



^^ AGPa + 7r« + ^ = (6.20) 



dPa dPa 

Flavour mixing in this context means that in general 5*0 = 5*0 [c^]. In the case 
of U(2) chiral fields and the presence of the 't Hooft term it holds that Sa ~ ctq, 
whereas for U(N), A^ > 3, chiral fields the 't Hooft term leads to a mixing, such 
that the Sa are in general non-linear functions of the aa- 



PlSee Witten (1979a), Veneziano (1979) and Alkofer and Zahed (1990) for a dis- 
cussion of a minimal t/^(l) breaking expression that is free of fiavour mixing. 
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Because we consider here the case of SU(3) symmetry and no £^g^-term we can 
use immediately the stationary phase equation (6.19,6.20) in order to ehminate 
the Sa, Pa- We obtain for the Lagrangian after a rescahng of the meson fields 
4>a -^ 94>a the form 

"^NJL = q[-ip + m + g{aaXa + il^'^aK)q + -l^'^{(^a(^a + T^aT^a)) (6.21) 

as one would obtain also from a Gaussian integral multiplicator (Eguchi 1976). 
After Grassmann integration over the quark fields one obtains finally the effec- 
tive action in terms of classical meson fields, corresponding to a one-fermion loop 
approximation : 

Seff = -Splog(-i^ + m + g{aaXa + ^IhT^aK) + ^/"^ / d'^^icfaCTa + T^aT^a) (6.22) 

This is expression is in a sense formal, because we have to apply a specific regu- 
larization scheme. We will choose here the double step proper time scheme, given 
by eq. (2.30) where (/)(r) = c6'(l - 1/A?) + (1 - c)6'(l - 1/A^) contains now two 
more parameters. This additional freedom will be used later in the mesonic sector 
to fix the non-strange current quark mass to the preferred value. 

6.2. Fixing of the parameters 

As we have said already the expression (6.22) has to be regularized and as we 
know from sect. 2.2. it should be done from the very beginning in order to avoid 
discrepancies. For pedagogic reasons the regularized form will not be written out 
explicitly. From the action eq. (6.22), which always can be cast into the form 

SefM] = J d^xV,ff[4>] + ^Z[4>]d^<pd^<P + . . . , (6.23) 

one obtains for the non-derivative part of eq. (6.23) the effective potential V'e// 

f rJ^k 1 

Veffi^a, TTfl) = -Tr / -777-rj(^ + 171 + g{aaXa + h^'^aXa)) + -j^'^ic^aC^a + T^aT^a) 
Jreg [^T^j ^ ' ^ 

(6.24) 
which differs for homogeneous fields from the effective action only by the four 
dimensional space-time volume ^^. Actually, spontaneous breaking of chiral sym- 
metry takes place and the ag and ag are the only candidates of fields, which can 



^''JCompare with sect. 2.6., where the parameter fixing is done with using the 
effective action only. 
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acquire a non- vanishing vacuum expectation value. This follows from conservation 
of parity, strangeness and isospin. 

However it has been shown (Pagels 1975) that the vacuum is SU{2)ji^ SU{2)i 
invariant in the chiral limit, if ctq and ag have non-zero vacuum expectation values. 
From Goldstone's theorem (Goldstone 1961) we know that the number of broken 
symmetries of the Lagrangian is related to the number of massless Goldstone 
bosons. Actually the physical world knows 8 pseudoscalar light particles, namely 
TT, -ftT and i]. If we neglect the t/^(l) breaking by instantons, even the ry' can 
be regarded as Goldstone boson. Since other would-be Goldstone bosons are not 
known experimentally, the true vacuum state of the NJL model, if it is physical, has 
to be SU(3) invariant in the chiral limit in accordance with Goldstone's theorem. 
Fortunately the NJL vacuum is indeed characterized by a large vacuum expectation 
value of the ag and a vanishing one for ag. This situation is changed for explicit 
symmetry breaking, like the presence of the 't Hooft determinant, where the t/^(l) 
symmetry is broken, or current quark masses, where the SU{3)y symmetry is 
broken. 

Now from the actual form of the potential (6.24) the non-trivial stationary 
phase conditions are given by 



dV, 



eff 



daQ 



dV, 



eff 



dag 



0. (6.25) 



After some algebraic manipulations, one is left with 



^iH'^-^) = 8Ncg^h{Mu) (6.26) 



//2(i - ^) = 8iVc//l(M,) (6.27) 



M. 



s 



where we have set Mu = y^go-Qy + mi H — y^igcrgy + m^) and Mg = y^go-Qy + 
mi - -/^{9(^8v + ms) and 

h^^^)=l ,f ^4 .2 1 ,,2 (6-28) 

Jreg (27r)* k^ + M/ 

Mesonic Spectra 
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The masses of the pseudoscalar and scalar meson nonets follow directly from 
the second functional derivative of the effective action. This gives an expression, 
which is up to a normalization factor Z^i, proportional to the inverse of a bosonic 
propagator. So it follows from general arguments that the analogue of eq. (2.40) 
is 

Zab{Q')SHPl-P2)Snaipi)5na{p2) " . = (p' + ^Ib) .2 .2 (6-29) 



Pl = -P2,P^=-'?^ 



p^^—q^ 



As described in sect. 2.5. one can now either use the normalization point q = 0, 
which correspond to a derivation of the masses from a gradient expansion, or at 
the on-shell point q = fn^^, so that the free propagator is defined with its physical 
mass. 

Defining 

/2(M„M„q2)=/ "^ ^ . ^ -^ — ^ \ (6.30) 

^ ' ^^ Jreg {2nf k^ + Mf {k + qf + Mf ^ ' 

we obtain for the normalization factor 

Z^,, = 4Ncgh2{Ma, Mf,, q^ = -m^^) (6.31) 

depending on the fiavour content of the given meson. Therefore the corresponding 
fields have to be rescaled according to (/>"' = Zaa (p'^- This gives 

For the ryg and ryg, the situation is more involved, since we have again mixing 
between them. However, similiar to the PCAC relations (6.9,6.10) above, we can 
remove the mixing by the same ideal mixing angle 6 = —35.26, if q = 0. Then 
the masses of the new ry and r]' are given by 



ml = ^^^^ 2ml -ml (6.34) 
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One can formally assign masses to the mixed rjQ and rj^ states and gets 

^Vo - o"^^ + o^K (6-35) 



m 



2 _1 2 I f _2 

V8 - 3"^7r + 2 



-i^rn^ + n^j^ (6.36) 



Within the NJL model it is clear that it is the r^g particle, which fulfills the Gell- 
Mann (1962) and Okubo (1962) (GMO) relation 

ml + 3m^g - 4m^ = (6.37) 

In nature however it is the ry, which should fulfill the GMO relation. This problem, 
however, can indeed be solved if one includes the 't Hooft interaction, which has 
the main effect of pushing the r]Q mass, such that the r] and r]' are driven back to 
their group theoretical states ?7o ^^^ V8 (Hatsuda and Kunihiro 1991). 
Decay Constants 

From the general form of axial vector matrix elements between pseudoscalar 
mesons and the vacuum < | A" | 7r"(p) >= —ipfifaT^^'ip), we immediately deduce 
the pion and kaon decay constants as 



Mu yl/2 , _ Ms + Mu ^1/2 

2(7 



f^ = -^Zi^\ fj^ = ^^-Jiz\l' (6.38) 



Chiral Perturbation for m 



s 



Because we will treat the baryonic sector perturbatively in the current quark 
masses, we will give here also the mass ratio for kaons and pions perturbatively in 
first order in rus 

!l| = !I^£±^ + or^)+of^V (6.39) 

ml 2mo \Mu) \M.J ^ ' 

It was already noted by Hatsuda (1990) and Hatsuda and Kunihiro (1991), that 
this relation gets large corrections order by order in perturbation theory. How- 
ever there are cancellation effects of the non-linear terms in the denumerator and 
numerator of the exact expression for rn-^/mq-, so that the exact result almost 
coincides with the approximate relation (6.39) (Schneider 1994). That means that 
although perturbation theory for the whole vacuum does not work, we are not a 
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priori going into severe troubles if we make use of eq. (6.39), which coincides with 
the exact answer. Later in the baryon sector we wiU examine the vahdity of per- 
turbation theory for mg again. So we summarize that by this equation the current 
quark mass ratio is fixed to nis/mQ ~ 24.5 for given experimental mesonic data, 
ttitt = 139MeV and m^ = 496MeV. Chosing our regularization function (/}{t) 
such as to reproduce the most reasonable value of ruQ ~ 6MeV, this corresponds 
immediately to nis ~ IbOMeV, which is also a reasonable value infered from ex- 
aminations of hyperon spectra. So the quark masses, that will be needed to fit the 
hyperon splittings in Sect. 6.3, will be compared with the ratio given by eq. (6.39). 
Furthermore the kaon decay constant fx equals /jr in this approximation and for 
the constituent quark masses one obtains M^ = M^ = Ms- 

6.3. Collective Quantization of the SU(3)-Soliton 

In sect. 4.2. we have dealt already the quantization of rotational modes in 
the case of SU(2) symmetry. Now in the case of the SU(3) group, the procedure 
is a little bit more involved. This is a refiection of the fact that SU(3) is now 
a rank 2 group, whereas SU(2) has rank 1 and that the configuration space for 
the SU(3) rotations is now restricted due to some trivial embedding of the SU(2) 
isospin subgroup into SU(3) (Mazur et al. 1984). Especially this embedding, 
which was proposed first by Witten (1983b), ensures in the end that only the 
physical representations emerge as the lowest possible ones (Balachandran et al. 
1985, Chemtob 1985, Mazur et al. 1984). 

But first we concentrate on deriving the left and right generators of the group. 
According to Witten, we make for the SU(3) chiral field the following ansatz : 

/ U2{x) 0\ 
U{x)=i^ Q J (6.40) 

where U2{x) = (o"(2)+^757rr)//7r and the SU(2) sigma field a/2) is defined according 
to a/2) = ^Ov 3 + ^^8- Then the aQ and ag fields themselves are constrained by 

"^OV 3 7q^8 = /tt? cis long as there is no fiavour mixing interaction present. This 

gives the 1 in the lower right corner of U{x). However this constraint is not really 
necessary and recently solitons were found (Kato et al. 1993), where it is released. 
In order to define quantization rules, we again introduce the time dependence by 
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writing 

U{x, t) = A{t)U{x)A\t) (6.41) 

where A{t) is now a time dependent matrix of SU(3). The theory remains invariant 
if one performs on A{t) the symmetry operations 

A{t) ^ A{t)G,paU G,pat e SU{2) (6.42) 

and 

Mt) ^ Gji,,A{i), Gfi,, e SU{3). (6.43) 

The right multiphcation of A{t) can be identified as spatial rotation because any 
isospin rotation of the chiral field can be undone by a suitable space rotation. 
This is due to the symmetric hedgehog ansatz. Note that because U{x) commutes 
with the hypercharge group t/y(l), with generators expzggAg, the configuration 
space for the generalized coordinates of the SU(3) matrix A(t) is restricted to 
SU{2)j Cg> t/y(l) (Balachandran et al. 1985). However, following the work of 
Balachandran et al. (1985), Chemtob (1985) and Praszalowicz (1985), one can 
regard the symmetry of the collective coordinates as SU(3) and treat the freedom 
for right multiplication Uji{l) as a constraint for the states. 

Then one can perform the time dependent rotation of the effective action as 
described in sect. 4.2 and obtain 

S^ff = -Splog (dr + H + At(t)i(t) - z74At(t)mA(t)) (6.44) 

with 

H = -i-f4^ {idaz - mU{x)) (6.45) 

Because of the antihermiticity of A'{t)A{t)j any expansion oi S^ff in odd powers 
of this quantity constitutes a contribution of the imaginary part of the effective 
Euclidean action and is therefore a finite quantity. As in the case of the baryon 
number in sect. 3.2, it is reasonable not to regularize these quantities. This 
is because they are connected with topological indices, which otherwise are not 
exact integers. We will see this, when we consider e.g. the right hypercharge in 
this model. To this aim, we write the Maurer Cartan form A){t)A{t) as 

A^{t)A{t) = qaA^daA = '-qaC^^XA = '^^A^A (6-46) 
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where the qa are the coordinates of SU(3) and the C(/^ are the vielbeins, which 
fulfill 

C^''c\ = Si (c/)"'=(7"^ (6.47) 

and the important Maurer Cartan identity, following from the definition of the 
structure constant /abC of ^^e A-matrices, 

C^Bd5C\ - C\dsC\ = -fBAEC\ (6.48) 

The effective Lagrangian after this SU(3) rotation is given by (in the chiral limit): 

^M = I^aIaB^B - ^B{U)n8 = lqa9apQp + Z'^Qa (6.49) 

where B{U) is the baryon number of the system with chiral field U(x) and the 
metric 

^a/3 = C^/^ASC'/, Z" = -^C,8 (6.50) 

Note that in contrast to SU(2), there is now a term linear in the angular frequency 
Og. This is due to the imaginary part of the effective Euclidean action, which is 
vanishing for SU(2). As mentioned above such a term poses a problem for the 
quantization, because the corresponding generators are constraint in this case (see 
e.g. Balachandran et al. 1985 for a discussion of this). 

Note that this Lagrangian also describes a particle moving in a monopole 
gauge field (Jackiw 1983), given by the gauge-variant Z^. As a result, it leads to a 
quantization condition for electric and magnetic charge. In the present model it is 
the Wess-Zumino term, which is the leading order term of the gradient expanded 
baryon current B{U) (2.38), that plays the role of the monopole. The consequence 
of this gauge variant term in the present quark model is a quantization condition 
for the number of colors (see App. D for details). 

But let us proceed step by step. From eq. (6.49)we can define canonical mo- 
menta by 

dL 

TTa = 77- 6.51 

Oqa 

and right generators by the symmetrized form (Toyota 1987) 

Ra = -l{na,C^a}- (6-52) 
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In the same way one can define left generators by using the form A(t)A^{t) with 
some vielbeins i?"^ and 

La = -l{7Va,E^,}. (6.53) 

From the Maurer Cartan identity and imposing the canonical quantization by 

ka, /] = -iS^: ka, TT^] = [fo, q/]] = (6.54) 

one can derive 

[Ra, Rb\ = -ifabcRc, [La, ^d = ^fabcLc (6.55) 

and the action on SU{3) transformation matrices as 

[Ra, A] = -A^Xa, [La. A] = ~XaA (6.56) 

Now we see from eq. (6.55) that the left generators La obviously obey correct 
commutation relations for SU{3) and can be identified with left fiavour rotations, 
because of the remarks at the beginning of this paragraph and eq. (6.56). The 
generators Ra can be related to some right transformations according to eq. (6.56) 
and act on the space of spatial rotations. Because the eighth component is con- 
nected with the baryon number, the space is also called the spin-baryon number 
space. 

We set the baryon wave functions according to Salam and Strathdee (1982) 
as 

*(A) = r]Vn < i I L'(")(A) | j > (6.57) 

where n is the dimension of the representation and r] is a suitable phase factor. 
The basis states < i \=< /, 13, Y \ and | j >=\ J, J3, Yji > carry the fiavour and 
spin quantum numbers. The Wigner function D^^'{A) is a matrix of the adjoint 
representation of SU(3) and is given by 

D^^liA) = ^tr AhAAXB (6.58) 

Using the commutation relations eq. (6.56) one obtains the action of the generators 
on the wave function 

[Rc, Dab] = ^IcdbDad, [Lc, D^b] = ifcADDoB (6-59) 
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From these one can deduce the action of the generators i?^, L^ on the wave func- 
tions. Furthermore, one can identify left indices of -D^^j as flavour quantum num- 
bers, chosen as hypercharge Y and isospin /, I^. The right indices therefore cor- 
respond to the right hypercharge Yji and spin J, J3. Now we come back to our 
rotated Lagrangian L^j^^. Using the properties of the hedgehog ansatz, we find 

(Ii6ab forA,B=l,2,3 
Iab = h^AB for A,B=4,5,6,7 (6.60) 

i for A,B=8 



where Ij^b is given by 



Nc f duj 



-^A- — rrr^B 



(6.61) 



.iuj + H iuj + H 

and the separation into valence and sea-part is analogue to Sect. 4.2.. The right 
generators get explicitly 

- -IiQa, A=1,2,3 

< -h^A, A=4,..,7 (6.62) 

JV3 = ^y^ A=8. 



Ra- ^''^ 



dflA 
so that we can go from L to the hamiltonian H by 



H = -J2Ra^A-L (6.63) 

A 



and obtain 



= Md + ^C2(5t/(3)^/i) + (^ - i^^C2{SU{2)r) (6.64) 

where the (72(5*^/(3)^/^) is the Casimir operator of SU(3) and C2{SU{2)jij is the 
corresponding one for the right SU(2). The eigenstates of these operators are the 
irreducible representations, which are usually labelled by the quantum numbers 
p and q, according to the rank 2 of SU(3). From the construction of states, 
which was done e.g. by de Swart, it is clear that Yji = 1 restricts the possible 
representations to those, which contain a basis state with Y = 1. Therefore only 
triality zero states survive, for which (p — q)m.od 3 = 0. As lowest possible one 
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these are the octet {8} = (p = 1, g = 1) and the decuplet {10} = (p = 3, g = 0). 
With formula eq. (6.64) at hand, we can determine the masses for the center of 
the octet and decuplet states. In order to remove the degeneracy of the states, 
within the multiplet, we have to switch on a finite symmetry breaking via a finite 
current quark mass. In this way, we obtain 



where S is the SU{2) sigma commutator, defined in sect. 2.1 and -ftT^^j are the 
so called anomalous moments of inertia, given by (compare with Park and Rho, 
1988, in chiral bag models) 



Kau = —i — - / — tr 
^^ 4 i 27r 



1 X 1 X 



(6.66) 



_iuj + H iuj + H 
and obey a similiar structure like the Tab'- 

(Ki6ab forA,B=l,2,3 
Kab = K26AB for A,B=4,5,6,7 (6.67) 

{ for A,B=8 

Because they originate from the imaginary part of the effective Euclidean action, 
they need no regularization. The quantization condition changes to 

(-{I^n^-^KiD^A), A=l,2,3 
Ra = -^ = -ih^A - '-^K2D8a\ A=4,..,7 (6.68) 

so that one obtains for the symmetry breaking part of the collective hamiltonian 



rsb _ _ ^2^ 2m, 



'^ch-rr.jfY ^ 



(^-f)J<(^)^.. 



™» .Efl-I)g'(A)U^^I)<«'(^). (6.69) 



3 ruu + rrifi \ »» y 3/2 

AT 

where we have used the relation X)^ D^a^-A = Lg = -^Y. The different values for 
the moments of inertia /i, I2, Ki and K2 from the selfconsistent solitonic solutions 
can be found in Tab. 6.1. for two values of the constituent quark mass M. 

Now we are in the position to evaluate the mass splittings within the hyperon 
multiplets. This will be done in the next section by evaluating the expectation 
value of the collective hamiltonian in a given baryonic state. 



6.4. Static properties of SU(3) Baryons 

After writing down the collective hamiltonian, we can evaluate it in different 
ways. First we can set baryon wave functions as described in eq. (6.57), which 
are eigenfunctions of the symmetric hamiltonian H Jj^ . All we have to do is 
then to sandwich the D-functions of the symmetry breaking part of the collective 
hamiltonian H^j'^ between baryon wavefunctions (Adkins et al. 1983). On the 
other hand, one can parametrize the right generators by differential operators 
of the Euler angles of the given representation (Yabu and Ando 1988, Park et 
al. 1991, Park and Weigel 1992). Then it is possible to evaluate the hamiltonian 
between the exact wave functions, which can be perturbatively expressed as a series 
in the symmetry breaking quark mass nis- After we have clarified the different 
procedures, we will comment on the numerical evidence of the exact treatment. 

The Perturbative Treatment of the Collective Hamitonian 

Having collected the leading terms in view of a I/Nq expansion for the collective 
hamiltonian, we have to sandwich this operator between the wave-functions of 

the baryons. Because of the symmetry breaking terms in -ff^^;;, the functions 

(n) 
Dj^^{A) in octet and decuplet representation are not exact eigenstates of the full 

hamiltonian. In principle, one has to write down a perturbative expansion 

\B>^\B,R>+ Z ^^'''^^''"Jf-''~'' \B,N>+... (6.70) 

Ny^R ^R ~ ^N 

where R is the lowest representation in which the baryon B can be found and the 
summation goes over all higher dimensional representations N, which have non- 
vanishing overlap with the ground state in R. So it is clear that the proton, e.g., 
is a sum of a proton in an octet representation and corresponding states in 10 
and 27. Furthermore these states enters in eq. (6.70) linearly with the symmetry 
breaking parameter nis- Because only the symmetry breaking terms in the collec- 
tive hamiltonian have a no n- vanishing overlap with these higher representations, 
they appear in the mass for the first time in the second order of nis- So one can 
define a perturbative treatment of m^ in first order, in which the first order nis- 
correction of the wave-function (6.70) do not appear. The effects of these terms 
will be discussed in the next section. 

In order to evaluate H^^^ in the baryon states | B,R>, we need the integral 



89 



over 3 D-functions. These integrals are well known (Blotz et al. 1993b) and are 
summarized in Tab. 6.2 (cf. De Swart ,1963, for a general discussion). As a result, 
we can express the mass splittings in terms of the 2 quantities 

A = |^^E + ™i2^-3^) (6.71) 

and 

S = ms— i. (6.72) 

h 

This is little bit astonishing, since we have three operators Y, Dgg and J^A^l 2 3 -^S^-^A- 

However due to a non-trivial group theoretical property (Blotz et al. 1993b), the 

two quantities A and 6 are sufficient. Then we obtain for the splitting between 

the center of the multiplets (experimentally ^ 230MeV) : 

Ag-io = ^ (6.73) 

which is solely given by the difference of the Casimir operators C2{SU{3)j^/]^) 
and C2{SU {2) ji) in octet and decuplet representation. One should note that this 
formula coincides with A^ — A splitting in SU(2) theory. This means that the 
same expression has to be considered as a different physical quantity depending 
on the use of the SU(2) or SU(3) collective quantization. Then we can write for 
the masses of the strange baryons relative to the mass of the E*, which we fix for 
the moment to the experimental value, as^ J 

Am^r = -—A -5- Ag-io (6.74a) 

AmA = -—A - Ag-io (6.746) 

Ams = Yo^ ~ ^^"^" ^^'^"^""^ 

Am= = \/\ + d- Ag-io (6.74d) 

^ 5 

AmA = — A - 5 (6.74e) 

8 

Ams* = (6.74/) 



^ J Semi classical quantization is known to yield always values for the masses which 
are by several hundred MeV too high. Remedies to this by means of band-head 
and centre-of-mass corrections will be used below. 
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Ams* = -A + S (6.74fif) 

8 

Am^ = jA + 26 (6.74/i) 

The reason that we fixed the E* to the experimental value is that the quantization 
always yields too large values for the masses if one does not employ band head 
and center of mass corrections as discussed below. 

From the formulas eqs. (6.74) above follow the mass relation of Guadagnini 

(1984) 

m^* — m^* + m]\[ = -{llm/y — Sm^;). (6.75) 

" 8 

It is interesting to note that it is obtained in the pseudoscalar Skyrme model 
only by introducing the hypercharge operator and several coefficients (Guadagnini 
1984) by hand. In our approach this term arise naturally in the theory and the 
coefficients in front of all of them are completly determined by the selfconsistent 
soliton solution. 

Furthermore one obtains the Gell-Mann Okubo relations (Gell-Mann 1962, 
Okubo 1962) 

2(m7v + ?^s) = 3mA + "^S (6.76) 

and 

rufi — m^* = m^* — mj^* = m^* — m/\ (6.77) 

automatically in this approach by means of eq. (6.74). One should note that 
these relations are mass sum rules which rely basically on the fact, that the SU(3) 
flavour symmetry breaking part of the strong interaction can be treated in 1st 
order perturbation theory (Cheng and Li 1984). Their validity shows therefore the 
consistency of the whole quantization procedure. On the other hand the values 
of the hyperon masses themselves depend on the dynamics of the model which is 
considered. In Tab. 6.1 we have collected the difference of the theoretical mass 
from the experimental one for two values of the constituent quark mass M. As one 
can see from the table for the value nis = IbOMeV, which correspond to eq. (6.39) 
for Too = 6.1MeV, the splitting prediction is better than 50MeV, except for the 
O, which drops out with 80MeV. However it is instructive to increase nis to a 
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value nis = 200MeV, i.e. this corresponds to enlarge nij^ to ~ 570MeV (exp. 
nix = 496MeV, which gives a very good agreement with experiment up to 20MeV 
for all baryons under consideration. 

Comment on Validity of Chiral Expansion 

In Sect. 6.1, following (6.39), we mentioned the failure of chiral perturbation 
theory for some vacuum parameters. Now we are in a position to answer this 
question for the baryonic sector and it will be done here for the case of the total 
masses as an example. Other observables are considered in the literature (Blotz 
et al. 1993a, 1993d). These two sectors need not to have the same behaviour, 
because the vacuum sector is dominated by different integral equations compared 
to the soliton. Actually we can write the classical soliton mass in perturbation 
theory for mg as 

Ma = M(o) + nisM^^^ + mlM^2) + 0{ml) (6.78) 

with the result within the present model: 

= 1250 + 426.9 -26.0Mey (6.79) 

for the preferred values of M = 418MeV, E = bQMeV , Nq = 0.668fm and ms = 
186MeV. So one can say that at least for the action the validity of the chiral 
expansion seems to be more reliable for the baryon sector than for the vacuum 
(Hatsuda 1990). 

The Yabu-Ando Diagonalization Method 

In contrast to the former method of evaluating the collective hamiltonian be- 
tween the symmetric eigenfunctions of the Casimir operators, Yabu and Ando 
(1988) developed a method for treating the baryonic wave functions to all orders 
in nis- The idea is to the express the rotation matrix A by its eight Euler angles 
and derive also for the right generators i?^ an explicit differential operator form 
in terms of these Euler angles. 

If we adopt the same parametrization of the rotation matrix A as discussed 
by Yabu and Ando (1988), we can write 

A = R{a, /3, 7) exp {-iuX4)R{a', /?', 7') exp (-zpAg/ Vs) (6.80) 
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where the /?(«, /3, 7) is the Euler-angle rotation matrices of SU(2)-isospin. Because 
of the hedgehog ansatz for the chiral field, isospin rotations and space rotations 
are intimately connected in a way, that any isospin rotation can be undone by 
a corresponding rotation in coordinate space. Therefore the R{a' , /?', 7') is the 
SU(2) spin rotation matrix. Now we would like to obtain an expression for the 
right generators Ra in terms of these eight Euler angles, conveniently written as 
ctfl. Following the work of Nelson (1967) and Park et al. (1991), we make the 

ansatz 

d 

Ra = idab{a)-^, (6.81) 

such that Ra are a linear differential operator of the aa- The unknown matrix d^b 
can be determined by inserting this ansatz into the definition of the Ra 

ARbA^ = AXbA^ (6.82) 

Now both sides of (6.82) can be explicitly worked out and give with eqs. (6.80, 6.81) 

XaEac{a)dcbia) = Cab{a)K 

with the some known matrices Eac{,(y) and (7^5(0;). Comparison yields dab{(y) = 
i?Q^ (a)C(,5(a). Their rather lengthy forms are summarized by Park and Weigel 
(1992). Then one can directly act with these operators on the wave functions, 
which are themselves D-functions and therefore also given in terms of A. 

The result of this method can be seen in Fig. 6.1, where the deviation of the 
theoretical prediction from the experimental one is shown in dependence of the 
strange current quark mass. Although mg is fixed to be ~ 150MeV from the 
mesonic sector, it is instructive to see that although the predictions for rus = 
150MeV are already good, a slight increase to nis — 187MeV leads to an almost 
perfect agreement (see also Tab. 6.3). The difference to the former perturbative 
treatment is in addition illustrated in Fig. 6.2, where the Yabu-Ando method is 
compared to the perturbative treatment for two particles of the multiplets. As it 
is clear from these curves, the Yabu-Ando method gives a significant contribution 
to the masses for current masses above 100 MeV. Similiar calculations have been 
performed by Weigel et al. (1992b). The formalism of these authors differs from 
the present one in the fact that they do not follow stricly a perturbative expansion 
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in 1/Nc and nis- Hence they also have different constituent quark masses for the 
strange and non-strange quarks. The resulting numbers are quite close to those 
presented here. 

Furthermore Praszalowicz et al. (1993a) made a refined calculation in the 
lines of the preceding sections with the inclusion of isospin breaking current quark 
masses. As a result they found that one can make a prediction for the hadronic part 
of the isospin splitting within the multiplets, which agrees with the experimental 
data (Gasser and Leutwyler 1982) with a high accuracy. In this work it turned 
out that the anomalous moments of inertia play again a crucial role in order to 
predict the splittings such as the neutron proton mass difference, which could 
be explained neither in the former SU(2) pseudoscalar Skyrme model, where it 
vanishes, nor in a U(2) or SU(3) extension (Jain et al. 1989) of the model. In 
Fig. 6.3 the hadronic part of the isospin splitting within the octet is shown and 
compared with the experimental predictions from Gasser and Leutwyler (1982). 
In addition it could be shown by Praszalowicz et al. (1993) that the complete 
agreement between theory and experiment within the experimental error bars is 
achieved for that m^ — niu difference, which follows from the mesonic sector from 
the masses of the charged kaons and pions. This can be viewed as a reflection of 
the fact, that perturbation theory works quite well for the small isospin breaking 
parts. In addition it was shown by Blotz et al. (1993e), that the recently measured 
Gottfried sum, which was considered in the present model for the case of SU(2) 
by Wakamatsu (1992b), comes out in the SU(3) version quite close to the value of 
the NMC measurements (NMC, Amaudruz et al. 1991). In addition the E-term 
comes out to be E ~ 45MeV, which is again close to the experimental value of 
Gasser et al. (1991) 

Total Masses and Zero-Point Energy Corrections 



By now we have discussed the mass splittings within the baryon octet and 
decuplet and the splitting between the octet and decuplet. As mentioned already 
the absolute energies of particles come out by several hundred MeV too high in 
the semiclassical approximation no matter which chiral model is used. Subtraction 
mechanism, as they are discussed by Pobylitsa et al. (1992) in SU(2) and by Blotz 
et al. (1993b) in SU(3) (compare with Jain et al. 1988b for a discussion in the 
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Skyrme model) can in principle cure the situation. In order to obtain reasonable 
values for the absolute masses of the hyperons one has to subtract the spurious 
zero mode energies similiar to the way described in sect. 4.4.. The corresponding 
terms of the rotational zero modes in SU(3) and of the translational one read 

AM-* = ^(C2(5t/(3))) + l[^,~Jl) <^2(^t/(2))) 
_ 1 7 19 

^^iransl = ^j^^J (6-83) 

where we used the values of the Casimir operators for the fundamental represen- 
tation. These are {C2{SU{3))) = 4Nc/3 and {C2{SU{2))) = 3Nc/4. The transla- 
tional zero mode subtraction coincide for SU(3) with the SU(2) result AM 
of (4.66). Subtracting these terms from the classical soliton mass M^i gives for the 
center of the octet and decuplet 

Mg = M^i + ^^ + ~- AM'^«* - AM**^""*^ (6.84) 

/2 4 /i 8 

MiQ = M^l + ^^ + ^—- AM™* - /^M^-ransl (ggg) 

h^ /i 8 
If one performs these corrections to the classical energy analogous to (4.65) in the 
present case of SU(3) (Blotz et al. 1993b) the absolute mass of the E*-particle 
changes from 2262 MeV to 1494 MeV (experimentally m^* = 1385Mey). Thus 
with these corrections the masses of all the hyperons of the octet and decuplet are 
reproduced with a constant shift of ^ lOOMeV, being in fact a quite impressive 
result. 

Expectation Value of Axial Currents 

In a similiar way to chap. 5 one can evaluate the expectation value of the SU(3) 
axial vector current operator. However, as it was mentioned in the SU(2) case 
already, there are important 1/Nc corrections to vector and axial vector currents. 
These are entirely due to the fact, that the collective operators O and A^X^A in 
general do not commute and have to be explicitly time-ordered (Christov et al. 
1993b, Blotz et al. 1993d). Without going into details we will nevertheless give 
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here a shortcut derivation. From the path-integral one obtains for A^(x), a 
0,3,8: 

< Alix) >= — --[Sp(to)log {d^ + H + iVLe - i^^A+mA + i s(x)747^75A+AM 

(6.86) 
where Sp/^^^A means now that the time dependent operators (A(t),0) within the 
trace have to be time-ordered, before performing the trace in time direction. This 
however disagrees with the usual definition of the trace mapping and its origin 
is the fact, that we are not going to perform the path-integral over the rotation 
matrices A{t) in the functional sense (Blotz et al. 1993d, Christov et al. 1993b). 
Instead one usually switches to the operator formalism, as it was done throughout 
this work (cf. Dyakonov et al. 1988), which involves the evaluation of time-ordered 
products of operators. Because of this, the trace in eq. (6.86) has to be modified 
in a way, that respects an explicit time-ordering of the operators (cf. Blotz et al. 
1993d). To obtain a c-number value for the current, finally one has to sandwich 
the operator < vl^(x) > between suitable baryon wave- functions. 

In addition to the SU(2) case, eq. (6.86) contains the symmetry breaking term 
from the strange current quark mass. Therefore the expression for axial current 
up to the first order in the rotational matrices iQj^ and the strange current quark 
mass m is given by 



< A^'^ix) > = SP(to) ^ .rr mi^lbA^IaA 



84 + H 

- SP(to)^^^^747M75^^/aA^^^(5if + 0{SH^) (6.87) 

where the perturbation 6H is now given by 

SH = iVtE - i-fiA+mA (6.88) 

Evaluating this yields for the a = 3, 8 part of the axial vector coupling constant 
(Blotz et al. 1993d) 

AMaa 

h 
2iQi2 2z(545 

, 2M83 4m.sKi 
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Atyt 

' NssDasil - Dss) (6.89) 



and for the singlet part (Blotz et al. 1993a) 



.0_2v|43^3_,„,^48,(| 



Sa = , Rs - 4m,D>°' -J-M83 - iV83 (6.90) 



There we used the foUowing definitions for the various moments of inertia. For 
the anomalous moment M^^, we find 

^bc = ^ XI < ^ I ^3-^6 \m>< m\ Xc\n> TZj^{En, Em) (6.91) 

with 

^ /p P N 1 sign(£;n - li)- sign{Em - y) ,„^^. 

and where the chemical potential y lies always between the valence level and 
positive continuum of states, in order to describe a baryon number B = 1 system. 
For the proper time regularized normal moments A^^^. = A^^^ val + ^bc sea from the 
symmetry breaking, we find 

.J Nc^<n\ asXh \v><v\ Ac70 I ^ > /« oq^ 



z ^ £/|^ rjy 



and 

Nu. „„„ = - 
4 



^bc.sea = -rY. <^\ ^Sh I TO >< m I Ac70 I ^ > ^/^(^n, ^m) (6.94) 

n.m 



with 

^/3(^-^-) = ^I ^At) r J, 7^ ] (6.95) 

2^y7V JO ^t En - Em 

and for the antisymmetric moments Q^^, from the explicit time-reordering we define 

W6c ^bCjVal ' ^bc,sea- 

Qbc,val = ^ Z^ F^ ^ signE^ (6.96) 

z n -^n ~ -l^v 



and 



Qbc.sea = ^ H < "^ I (^3h | n >< n | Ac | TO > TlQiEn, Em) (6.97) 
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In the case that the regularization function (^(t) is given by (p{t) = J2i CiO{l — 1/A^) 
the IZq^Eu, Em) takes the simple form 

^ (J, P ^_ /■! da ajEn + E^) - Em exp {-[aEl + (1 - a)El]/K]) 
^^ "' "^^ " "' Jo 2n y^(IT^ aEl + (1 - a)El 

(6.98) 
In the infinite cutoff limit it again reduces to 

TleiEn, E^) = l ^ig"(g^.-M)-y"(gn-M) (6.99) 

Note that the regularization function TZQ{En, Em) in eq. (6.98) is now antisym- 
metric with respect to the states m and n in contrast to lZj^{En, Em) in eq. (6.92). 
This is a reflection of the fact that the matrix elements of eq. (6.97) are now an- 
tisymmetric with respect to m and n. 

Numerical Results 

Similiar to the case of SU(2) in Sect. 5.2., the terms with the antisymmet- 
ric Qfj(, expressions also serve here as large corrections to the lowest order terms 
(O^) and pushes especially the value of gj^, which was 50% too low without these 
terms, a little bit beyond its experimental value (see Tab. 6.4). Together with 
the values of g^ and g^, which has the interpretation of being the spin of the 
proton that is carried by the quarks, the contributions from the lowest order O 
and from the symmetric and antisymmetric rotational corrections are presented 
in Tab. 6.4 and compared with the experimental numbers from recent EMC and 
SMC experiments. One should note that (7^ should be interpreted as the spin of 
the proton (normalized to unity), which is carried by the quark. As it is clear from 
these numbers, the antisymmetric contributions from the time-reordering play a 
significant role. At the present stage of the work the NJL numbers are very good. 
The corresponding expressions emerge from the real part of the Euclidean effective 
action, though they turn out to be finite and need no regularization. Furthermore, 
and this is most important, they are not present in the usual Skyrme model ap- 
proach. So this is in a sense another refiection of the explicit quark degrees of 
freedom, which obviously cannot be described perturbatively in the chiral fields 
and which provides therefore a clear distinction between Skyrme type models and 
Nambu-Jona-Lasinio type models. 
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7. The NJL-model with vector mesons 



In this chapter we consider the influence of vector couplings on the NJL model. 
After giving some motivation, the bosonization and regularization of the extended 
model are presented. The vacuum and mesonic sectors are analyzed both assuming 
the mesons to be on-shell or off-shell. Finally a system with baryon number one 
is constructed and its solitonic solutions are analyzed. 

7.1. The Effective NJL Action with Vector Mesons 

Why Vector Mesons ? 

Since their discovery (Nambu 1957, Frazer and Fulco 1959, 1960), vector 
mesons have been intimately linked to the internal structure of the nucleon. In 
fact, the success of the vector dominance hypothesis (Sakurai 1960; Gell-Mann 
and Zachariasen 1961) and its later fleld theoretical realization through current- 
fleld identities (Kroll et al. 1967) in many hadronic reactions is obvious (for a 
review see e.g. Gourdin 1974). In addition, vector mesons play a crucial role in 
low energy nuclear physics, since they are responsible for the medium and short 
range NN interaction (Machleidt et al. 1987). It seems natural to ask how vector 
mesons may be implemented in a NJL type model. 

Up to now, we have referred to the NJL model with scalar-isoscalar and 
pseudoscalar-isovector couplings (cryr-version of NJL). The introduction of other 
couplings makes possible to describe a wider meson spectroscopy, since it is known 
that the cttt version of the model does not account for vector or axial vector degrees 
of freedom, i.e. the corresponding correlation functions do not posses poles. 

If additional vector couplings are included, the NJL-model indeed incorporates 
in a natural way important phenomenological principles found long before the 
advent of QCD: Sakurai's universality (Sakurai 1966, 1969) and vector meson 
dominance realized by means of current-fleld identities (Kroll et a. 1967). It also 
provides a relationship to effective low energy Lagrangians with vector mesons 
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(Kleinert 1978, Dhar et al. 1984,1985; Ebert and Reinhardt 1986; Wakamatsu 
and Weise 1988) both in massive Yang-Mills (Lee and Nieh 1967 ; Gasiorowicz 
and Geffen 1969; Meissner UG 1988) or in hidden symmetry (Bando et al. 1984, 
1988) representation. Furthermore, it includes the gauged Wess-Zumino term 
(Witten 1983) with the vector mesons interpreted as internal gauge fields (Dhar 
et al. 1984, 1985 ; Ebert and Reinhardt 1986; Wakamatsu 1989). In addition, 
almost all attempts claiming any parentage of the NJL model to QCD require the 
explicit inclusion of vector mesons (Dhar et al. 1984,1985; Cahill and Roberts 
1985; Schaden et al. 1990; Chanfray et al. 1991). An exception to this is the 
instanton liquid model (Diakonov and Petrov 1986). 

Coupling of Vector Mesons - Bosonization 

The extended SU{2)ji (g) SU{2)i ^U{1)b classically invariant NJL- mo del in- 
cluding scalar, pseudoscalar, vector and axial-vector couplings reads 



^NJL =q{^^ - mo)q + —[{qqf + [qmrqf] 



(7.1) 



It is important to mention that the terms multiplying the coupling constants Gi, 
G2 and G3 are themselves chirally invariant. Thus at the classical level the con- 
servations laws given in chapter 2. remain valid. Furthermore, since the coupling 
constants have dimensions of inverse mass squared it is convenient to choose them 
as follows 

n'^ 0^ n^ 

Gi = ^ , G2 = - — 2 ' ^3 = -^ (7.2J 

where (7,^, gp and (7^ are dimensionless coupling constants and nip and rM^j are the 
vector meson masses. The parameter n has mass dimension and it will turn out 
to be identical to the // introduced in chap. 2. 

Similarly to the an case the corresponding generating functional is made well 
defined by performing a Wick rotation. In order to keep track of Lorentz invariance 
the time components of the vector fields are rotated as well (cf. app. A) 

X -^ ix , uj {xQ,x)^iuj (x4,x), u}{xQ,x) -^ u}{x/^,x) (7.3) 
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Following the steps of chapter 2 we multiply in addition by the Gaussian factor 

An important point is that the Wick rotation makes the exponent in the Gaussian 
factor to have a well defined negative sign. Similarly to chapter 2 Euclidean 
indices will be understood unless otherwise stated. The Gaussian factor is chirally 
symmetric if the fields are assumed to transform as 



^Ai —"^PL 



Pn ^ Pn + a X p^ + (3 X a^ (7.5) 

a^ ^ a^ + a X d^ + /3 X p^ 
under the global chiral SU{2)^ (g) SU{2)i group. Performing the shifts 

"t'UJ 

Pi^^P^^- ^Qlf^Q^^Q (7.6) 

rwi 2 

one gets the semibosonized Lagrangian 

C'-^-^ix) = q(-i^ + fi-TT (o- + Z75r7f) +gp{j) + ^75) - - g^^ + mo)? 

+ y(^ + tt ) + -fip^ + fl/J + —U^^j, 

At the classical level, the quark part of this Lagrangian is invariant under local 
chiral transformations. We will see later that this local symmetry is broken by the 
regularization thus leading to a chiral anomaly. After integration of the quarks 
the extended action reads 

Sejf[(^. TT, w, p, a] = -Splog(-z^ + fi-TT (cr + ^75T'7f) +gp{j> + ^75) 7T " 9io^ + ^0 



2 

+ ^ / d^x{a^ + -') + ^l dUpl + 4) + ^ / d'^4 

(7.8) 

Again we treat the theory in the stationary phase approximation. The correspond- 
ing classical equations of motion give the current field identities: 

il f/ _^P-. jA_'^f 

9iu ^ gp 9p 



iB_'fri-u W-'^Pn r^-'^Pd f7ql 
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with X) , jj^ and Jrf being the baryon-, vector- and axial vector currents 

j]i = Qlf^2^Q^ J^ = Ql^lb^rq, J^ = q-ff^q (7.10) 

The parameters m^, gp, rriuj, g^j will be considered later. Due to the Wick rotation 
the effective action is a complex number which can be separated into real and 
imaginary parts 

S = ReS + ilraS , Re5 = ^(5 + 5t) , IraS = ^{S - S^) (7.11) 

On the basis of the combined symmetry operation Q = G75 with G the usual 
G-parity (charge conjugation plus rotation of 180 degrees around the y isospin 
axis) it has been shown (Doering et al. 1992) that the real part of the action is 
an even function of the uj field whereas the imaginary part is an odd function. 
In particular, if the uj field vanishes the Euclidean action becomes a real number. 
Notice that this argument does not apply anymore if one considers SU(3) fiavour. 

Regularization of the Effective Action - Chiral Anomaly 

As we have said the quark contribution to the semibosonized Lagrangian is 
formally invariant under chiral local transformations. On the other hand, a regu- 
larization has to be introduced to make the effective action finite. However, there 
is no regularization which preserves both vector and axial gauge symmetries si- 
multaneously. Thus there appears a chiral anomaly (Adler 1969; Bell and Jackiw 
1969). The subject of chiral anomalies is rather involved and we refer the reader to 
the recent review of Ball (1989) for a detailed discussion. (For a more elementary 
level see e.g. Petersen 1985). We just mention here that there is some mathemat- 
ical ambiguity as to which symmetry should be conserved and which should be 
destroyed. This freedom depends on the particular physical application. The fact 
that QCD is the fundamental theory underlying the NJL model suggests shifting 
the anomaly to the pure axial sector, i.e. to make use of a vector gauge invariant 
regularization. 

Dhar et al. (1984,1985) have proven that the proper prescription to achieve 
vector current conservation is to regularize the real part and not to regularize the 
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imaginary part, so that the effective action becomes 

1 roo dj- 7-1 1 7-1 1 / + \ 

S^ff = -- / — 0(r, A)Spe"^^ + -(Splog(zL') - Splog(-zL'T)] + mass terms 
■' ■' 2 Jo T 2^ ' 

(7.12) 

Then the real part of the one loop contribution is both vector and axial gauge 

invariant while the imaginary part exhibits an axial anomaly. This prescription 

leads to the original Bardeen's form of the chiral anomaly (Bardeen 1969), however 

with the vector and axial fields interpreted as dynamical degrees of freedom. 

In the case of the NJL model with dynamical vector mesons one gets the 
following anomalous Ward identities (Wess and Zumino 1971, Wess 1972), which 
are generalizations of the Noether theorem in quantum field theory, in Minkowski 
space 

d^j^ = 

d^^jj = 

d^^J^ = 2zmog75^g + j^9.9pe^,apd^u;''{ld''pf^ + gp{fP x p^ - a'' x d^)) 

(7.13) 
As we see, the anomalous contribution to the divergence of the axial current sur- 
vives, even in the chiral limit niQ = and in the absence of external fields. This 
is called an internal anomaly which does not even disappear if vector mesons are 
integrated out. In fact it has been proven by Wakamatsu (1989), Ruiz-Arriola and 
Salcedo (1993a), that this anomaly breaks some low energy theorems of QCD. 

7.2. Fixing of the Parameters in the Vector Mesonic Sector 

As it has been done in previous cases, the parameters have to be fixed by 
looking at different mesonic properties. As input parameters we use the pion weak 
decay constant /tt = 93MeV and the meson masses 171-^ = l39MeV, rup = 770MeV 
and mcj = 783MeV. 

The main new point in the meson sector is the occurrence of a mixing between 
the pion and the axial meson very similar to the one found in the early massive 
Yang-Mills approach (Lee and Nieh 1967). This causes after redefinition of the 
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physical axial field a finite renormalization of the pion kinetic energy (Kleinert 
1978; Ebert and Reinhardt 1986). As a consequence the corresponding cutoff 
increases with respect to the case without vector mesons. In addition, the axial 
mass acquires a substantial contribution from a partial Higgs mechanism. For 
simplicity we will work on the chiral circle. 

Real Part - Massive Yang-Mills and Hidden Symmetry Approach 

In the second order heat kernel approximation the real part of the effective 
action reduces to the following expression (Ebert and Reinhardt 1986; Wakamatsu 
and Weise 1988) 

- ll2{9ltr[{V^.f + {V^,f] +glu%} (7.14) 

where the regularization dependent integral /2(M) has been defined in eq. (2.28). 
The constituent quark mass M = Qtj- f-K has been also introduced. The covariant 
derivatives and field strength tensors read 

D^U = DJ^U - UDJ! = d^U - igpiV^U - UV^) 
where 



(7.15) 



U = —{a + if-n) 

V/ = ^(p, + a,)-r (7.16) 

The former expression for C can be identified with the old massive Yang-Mills 
Lagrangian (Lee and Nieh 1967) if one demands 
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^ ^ iS - 3^^ f^-^^' 



and also that rup and muj are the physical vector meson masses. Also it has 
been shown (Ball 1987; Wakamatsu and Weise 1988) that if one performs a field 
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dependent chiral rotation the hidden symmetry Lagrangian of Bando et al. (1984, 
1988) is obtained. As we can see at (7.14) the chiral symmetry of the Lagrangian 
suggests that the chiral partners p and a have the same mass. However, due to 
the spontaneous breaking of chiral symmetry there appears a term of the form 
Ma^, ■ d^n which would imply an unphysical decay process. This problem may 
be solved by introducing a new physical axial field a'^ = a^ + ^9^7? and fixing the 
parameter ^ in a way that the mixing for the new field disappears. This feature 
is also present in the hidden symmetry approach. In any case the solution to the 
A — n mixing leads to the following algebraic conditions (Ebert and Reinhardt 
1986; Wakamatsu and Weise 1988) 



.2 ^l ,,j2t 2 M2 6m^ 2 1 2 .7 1 ox 



and the axial meson mass is given by mj^ = rUp + 6M . After this the parameters 
may be fixed in the way described below. For a given constituent quark mass M 
the cutoff is adjusted to reproduce the pion decay constant f-^ = 93 MeV and 
the vector meson mass nip = 770 MeV. Then the rest of the parameters as for 
instance gp and m^ are determined uniquely. It is interesting to note that in 
the limit nip -^ oo the above conditions become identical to those of sect. 1.6. 
For Trip = \^M and if a constituent quark mass of M = 315MeV is chosen, 
the KSFR relation 2gpf^ = nip (1967) and the Weinberg sum rule ni\ = 2nip 
(Weinberg 1967) are fulfilled simultaneously. Let us finally mention that the fact 
that the massive Yang-Mills Lagrangian and further relations are obtained even 
for a finite cutoff indeed requires the continuation of vector fields into Euclidean 
space as it is done in this review. 

Imaginary Part - Gauged Wess-Zumino Term 

As we have already mentioned vector mesons generate an imaginary part for 
the Euclidean action. This has been computed in the low momentum limit (Dhar 
et al. 1984, 1985; Ebert and Reinhardt 1986), reproducing the vector gauged Wess- 
Zumino term (Wess and Zumino 1971; Witten 1983) after rotation to Minkowski 
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space 

+ 3i9pd^[daUUW^ - U'fdaUV^ + tgpiuW^UV^ - V^V^]) 

which is SU{2)y (g) U{1)b gauge invariant. It must be said that this part of the 
Lagrangian saturates the anomalous Ward identity (7.13). Hence higher order 
terms ought to be chiraUy gauge invariant. Another interesting point is that al- 
though such an effective Lagrangian preserves vector gauge invariance it breaks 
global chiral symmetry and hence breaks some low energy theorems such as the 
amplitude for the decay 7 ^ Stt (Wakamatsu 1989, Ruiz Arriola and Salcedo 
1993a). The main reason can be found in the a — n mixing which produces sys- 
tematic corrections in any vertex with external pions. This is a clear drawback 
of the model, also present in the topological soliton model. Nonetheless it repro- 
duces the correct result for the neutral pion decay via intermediate neutral vector 
mesons n -^ up -^ 77 and the experimental value for the strong decay u ^ Sir 
(Gomm et al. 1984; Kaymakcalan et al. 1985). 

Meson Propagators 

The calculation of on-shell meson propagators proceeds similarly as sketched 
in sect. 2.6. for the scalar-pseudoscalar case and has been treated in detail in 
several works. Here we will follow Jaminon et al. (1992). As in the heat kernel 
expansion, there appears an a — tt mixing term, which can be diagonalized after 
a proper redefinition of the axial field. In summary, the following conditions are 
obtained 

2 _ , ,2 4:NcF{-ml) 2 _ 1 2 _ 1 



(7.20) 



with the proper time regularized functions F(q ) and S{q^ 



S(,^) = J-j J\ f (1 - u^) /r ^e-[«=+i(l-"^),V ,7.22) 
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It should be mentioned that the approximation S{—mp) ~ 5'(0) = ^/2 and 
F{—m^) ~ F(0) = I2 corresponds to the heat kernel approximation as can be 
seen by comparison of formulas (7.20) with (7.18). Moreover, in the limit oi p — u 
degeneracy the relation gp = 2(7^ holds. Jaminon et al. (1992) have found that 
for constituent quark masses lower than M = 385 MeV the p meson is no longer 
a bound state. This is a direct consequence of the absence of confinement in 
the model. For these masses various remedies have been proposed. Takizawa et 
al. (1991) suggest to look for the complex solutions of the Bethe-Salpeter equa- 
tion, where the corresponding imaginary part represents the decay width into free 
quark-antiquark pairs. Jaminon et al. (1992) propose rather to take advantage of 
the approximate linear behaviour of the inverse propagator in the time-like region 
and to use a linear extrapolation into the space like region. For constituent quark 
masses around 500 MeV they find the approximate formula 

m\r^ml + 3.1M^ (7.23) 

which differs from the corresponding result in the heat kernel approximation (7.18). 
Finally, it has been found (Schueren et al. 1993) that these results do not depend 
strongly on the particular regularization scheme employed. 

Numerical Results 

Vector mesons do influence the vacuum properties since their masses and the 
A — n mixing enter explicitly the determination of the cutoff. In general the cutoff 
increases as compared to the case without vector mesons. This change is more 
dramatic if the parameters are fixed off-shell rather than on-shell (Schueren et 
al. 1993). The interesting result is that only in the latter method the results 
for the vacuum parameters seem to be reasonable and fairly independent of the 
regularization. For illustration we quote the results for M = 350 MeV and in the 
proper time regularization 

Q 1 

< qq >= -(271MeV)'^ mi = -(mu + mj) = 8.2MeV 

The effect of vector mesons on pure pionic properties, i.e. pionic radii and threshold 
parameters for nn scattering, has been investigated (Ruiz Arriola 1991a; Schueren 
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et al. 1993). A very interesting point is that vector mesons generate in a natural 
way an axial coupling constant for the constituent quarks gj^ lower than one as 
given by the expression (Vogl et al. 1990; Ruiz Arriola 1991a) 

^? = 1 - ^ (7.24) 

due to intermediate axial vector meson contributions to the axial current even in 
the leading order of the large Nc expansion. A value smaller than one has been also 
predicted in somehow different approaches and models (Peris 1991, 1992; Weinberg 
1992; Blotz and Goeke 1992) as subleading large Nc corrections. 

The main result found by Schueren et al. (1993) is that vector mesons do not 
change noticeably the calculated pionic properties, at least if the p meson mass 
is fixed to its experimental value. As in the (a, tt) case a strong dependence of 
the pion scalar and isovector radii on the constituent quark mass is found. In fact 
the precise numbers do not differ very much if vector mesons are introduced. In 
particular, fitting the pionic radii would require a quark mass of 250 MeV. As we 
have said, for this value of the mass the p and uj mesons lie in the continuum, and 
as it will be seen below no solitons have been found. 

More recently it has been found that a g^ 7^ 1 in the Nambu-Jona-Lasinio- 
model breaks the proper QCD anomalous structure. This can be only obtained 
for g^ = 1 and correspondingly gp = (Ruiz Arriola and Salcedo 1993). The 
problem arises whether vector mesons can be described in a NJL type model 
without breaking the QCD anomaly. 



7.3. SoLiTONic Solutions and Nucleon Observables 

The solitonic sector of the NJL model has been studied with quark couplings 
corresponding to p- mesons (Alkofer and Reinhardt 1990), p and A (Doering et al. 

1992 ; Alkofer et al. 1992), u (Schueren et al. 1992a; Watabe and Toki 1992; 
Alkofer et al. 1993) and p, A and uj mesons (Doering et al. 1993 ; Schueren et al. 

1993 ; Ruiz Arriola et al. 1993; Zueckert et al. 1993). The main difficulty arising 
in the calculations including the w-meson is the fact that the Euclidean static 
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energy is complex, due to the uj meson. In this section we wiU present a way from 
the complex valued Euclidean action to the description of selfconsistent solitonic 
solutions for hedgehog field configurations (Schueren et al. 1992a; Doering et al. 
1993 ; Ruiz Arriola et al. 1993). At the end of this section we comment on other 
approaches (Watabe and Toki 1992 ; Alkofer et al. 1993) 

Statement of the Problem 

To understand the nature of the problem related to the Wick rotation let H 
and if' denote the single particle Dirac hamiltonian and its hermitean conjugate 
given by 



-9toOJ4. + 9pl^{pA + Mlb) 



-9loUJ4: + 9pi^{pA + ^475) 



(7.25) 



H = h + i 
H^ = h-i 
respectively. The hermitean part h of the hamiltonians H and H' has the form 



h = -iaiVi + l3gT:{a + i75f ■ tt) + a, ( -gujUJi + gp-{pi + AjTs) J (7.26) 

Due to the Wick rotation H is in general non-normal, i.e. [H^ W] 7^ 0, and 
therefore H and W may not be diagonalized simultaneously. For time independent 
configurations the real part of the action reads after regularization 

ReSf = -^iV, r ^ r ^^(r. A) tr e(-+^)(-+^^) (7.27) 

■' 2 J-00 27r JO T 

and the imaginary part 

lmSf = -Nc —\trlog{iu + H)+trlog{-iu + W)} (7.28) 

2 J-00 27r '- J 

Now one would expect to proceed similiar to perturbation theory, i.e. one should 
evaluate the trace, compute the //-integral and rotate back to Minkowski space. 
Unfortunately, since the one particle Hamiltonian is non-normal, the u integration 
in the real part cannot be done analytically. Thus, the eigenvalues in the exponent 
should be computed for any value of the variable u. After that, the analytical 
continuation to Minkowski space should be done numerically. Needless to say that 
such an approach is not feasible in practice. 
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Total Energy of the B=l Soliton 



The calculation of the functional trace for baryonic systems is conceptually 
involved and has been discussed in full detail by Schueren et al. (1992b, 1993). The 
interesting feature is that the analytical behaviour of the spectrum is such that 
the rotation back to Minkowski space can be performed by solving the eigenvalue 
problems 

H^i;ti^)=4i^ti^) (7.29) 



where 



H 



± 



h± 



T 



9LuUJo-gp-{po + Ao-f^) 



(7.30) 



We just quote the final result. The total energy for a system with baryon number 
equal to one has been found to be 



E = ^val + Eji + Ej + E. 



mes 



(7.31) 



where the mesonic contribution 
.2 



E^es = ^j dM-' + -^ - fn) + ^ / dMpl + 4) + ^ / d'-4 (7.32) 



the valence contribution 



Eyal = ^{^val)^ 



val 



the sea real 



iV, 



Er-'^I: [eaR{ea, A) - e^^R{el A) 



and the sea imaginary contributions 

Ei = -— ^5Isign(ea)(e+-e^) 

^ a 

have been defined. We have also introduced the averaged eigenvalues 

_ _ 1 + _ 



(7.33) 



(7.34) 



(7.35) 



(7.36) 



where e^ and e^ are the eigenvalues of the Hamiltonians H^ and H^ respectively 
with the condition that for guj = they coincide ej = e^. It should be mentioned 
that the expression for the total energy has been checked asymptotically against 
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the heat kernel expansion in the pure a, n and u case for fixed meson profiles 
(Schueren et al. 1992b). From the former equations it can be seen that the valence 
part contributes to the total energy of the system as long as the averaged eigenvalue 
iyai is non-negative. We will see below that this is in perfect agreement with the 
way the valence quarks contribute to the baryon number. Since the valence part 
of the energy is discontinuous for e^g^i = 0, the question arises whether the sum of 
sea and valence contributions behave continuously. It can be shown analytically 
and also numerically (Schueren et al. 1992b) that continuity is fulfilled due to the 
contribution of the imaginary part. The remarkable effect is the strong evidence 
for a vanishing imaginary contribution if e^^; > 0. 

Let us remind that in the long wavelength limit the imaginary part of the 
action reduces to the gauged Wess-Zumino term eq. (7.19). Therefore the valence 
quarks coupled to all degrees of freedom seem to be complementary to the gauged 
Wess-Zumino term as long as e^^i becomes negative. 

Another consistency check is based on the infinite cutoff limit A ^ oo. In this 
case it is possible to perform the whole calculations also in Minkowski space. The 
result found coincides exactly with formula in the infinite cutoff limit. 

Currents of the B=l Soliton 

From the expression for the total action the corresponding mean field quark 
densities may be obtained as functional derivatives of the quark contribution with 
respect to the relevant fields. The expectation value of a bilinear in the soliton 
background is given by 

< q{x)rq{x) > = ^(e^;a;)^+^;(f)r^+;(x) 

- ^llsign(ea) ['!/'«(:?) r'i/'+(x) - tlj~{x)rt(j~{x)] 

^ " (7.37) 

- ^ E[^2(ea, A) + ea7^2(ea, A)]sign(ea) 

['0+(x)rV'+(:r) +t(ja{x)riij^{x)] 

with F = l,i75r, 7^. Here 7^2(^5^) and 7^2(^5^) represent the regularization 
function and its derivative as used in sect. 3.3. As we see there are three distinct 
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contributions to any mean field quark density: 1) a valence one provided the 
averaged eigenvalue Cy^i is positive, 2) a real sea contribution which depends on 
the cutoff and 3) an imaginary part contribution which is not regularized. This 
also true in particular for the baryon number density. Moreover, if vector fields are 
switched off the baryon density coincides with the non-regularized baryon density 
of chapter 3. 

Self- Consistent Solitonic Solutions 

The equations of motion for general static configurations can be obtained by 
varying the total energy of the soliton with respect to all the meson fields. This 
yields 

a(x) = —^ < q{x)q{x) > 



(7.38) 



7f(a;) = 


: 2 < Q{x)i-f5Tq{x) > 


f^{x) - 


= 2% < '"(->^''^^'(-' > 


a^(x) = 


2m'2 ^ ^(^)^^^5^9(a^) > 


u^ix) : 


= ^ < q{xh^q{x) > 



The first two equations of motion are formally identical to the ones for the (ttt case 
discussed in chapter 3. For vector mesons the equations of motion reproduce the 
current field identities, i.e. the w, p and a fields in the nucleon are proportional 
to the expectation values of the baryon, vector and axial currents in the soliton 
background respectively (see (7.9)). 



For practical calculations we use the hedgehog ansatz 

a = U<^{r)cosQ{r), Wa = XaUHr)smQ{r), p\ = 2e,kaX^ p{r), p^ = 0, 



c^^ = 0, uj^ = u{r), o^^ = 25],As{r) + 2xax\As{r)-\AT{T)), 4 = 

(7.39) 
Notice that the chiral circle condition a {x) + tt (x) = /^ is achieved if the polar 
field $(r) is taken to be a constant equal to one. It can be checked that this 
variational ansatz is a true solution of the equations of motion provided the Dirac 
sea contribution only involves closed shells in grand spin quantum numbers. 



112 



Numerical Results 

The mean field equations of motion have been solved iteratively for the non- 
linear case ($(r) = 1) by Schueren et al. (1993) and in the linear case ($(r) 7^ 1) 
by Ruiz Arriola et al. (1993) in a similar way as it was done in chapter 3. To 
simplify the discussion we will refer mainly to the full anup a system where the 
parameters have been fixed by computing the on-shell meson propagators and in 
the proper-time regularization. We remind that fitting the meson masses and 
the pion decay leaves the constituent quark mass M as the only free parameter. 
An interesting result is that vector mesons prevent the soliton to collapse if the 
chiral circle constraint $(r) = 1 is relaxed (Ruiz Arriola et al. 1993). This is 
a consequence of the repulsive nature of the u meson. For the simpler ayr case 
the collapse has been discussed in sect. 3.5.. The results for various mean field 
observables both in the non-linear and in the linear case are shown in Tab. 7.1, 7.2 
for different values of the constituent quark mass. 

As it can be seen, the total soliton mass is around 200 MeV higher than without 
vector mesons (see chap. 3). Unfortunately (see Tab. 7.2 the isoscalar radius and 
the axial coupling constant are 25% too high and 60% too low corresponding 
to experimental values respectively being nearly independent of the particular 
constituent quark mass M. Actually in the physical region the computed quantities 
do not depend strongly on the chiral circle condition, although the deviations from 
unity at the origin are significant in the linear model. Those become larger if the 
CO coupling constant is reduced and lead eventually to a coUapase of the soliton for 
values much lower than the ones nedeed to fit the uj meson mass. For instance, 
for M = 340 MeV a (7^ = 2.24 is required to fit the u meson mass, whereas 
collapse takes place below g^ = 0.8. The experimental isoscalar nucleon mean 
square radius can only be obtained for (7^ = 1.60. For this last value the averaged 
quark eigenvalue has been found to be positive. 

For illustration we show in Fig. 7.1 the self-consistent vector meson fields 
u), p, Ag and Aj^ in the non-linear case for the particular values M = 340MeV, 
Qp = 4.61 and Quj = 2.24. In principle, they resemble the corresponding fields 
found in other soliton models like e.g. the Skyrme model (Meissner UG 1988) and 
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the quark- meson model (Broniowski and Banerjee 1986), and their shapes do not 
depend strongly on the regularization scheme employed (Schueren et al. 1993). 

As we have discussed in chapter 3 one of the virtues of the NJL model is that 
no assumption is made a priori whether the valence quark picture is correct or not. 
This is in fact decided upon the results in the calculation of physical observables. 
This means if a calculation reproduces the relevant experimental data and if then 
a clearly separated bound single quark state exists in the energy gap between 
positive and negative continuum, then the valence picture is valid. In SU(2) and 
SU(3) calculations with no n- vector mesonic couplings all calculations by now yield 
this picture. 

When vector mesons are present the validity or invalidity of the valence quark 
picture depends on the sign of the averaged quark eigenvalue e^^;. From Tab. 7.2 we 
see that this is positive both in the linear and in the non-linear model. However in 
this context all calculations involving vector mesons performed so far do not allow 
to reach final conclusions because the nucleon observables are not yet reproduced 
well enough. In particular the sign of Cyg^i depends very much on the particular 
vector mesons included and the way the parameters are fixed. Schueren et al. 1993 
have shown that in the non-linear model and as long as the nucleon isoscalar radius 
is used as nucleon observables to be reproduced this question cannot be decided 
unambiguously. In other words, it is possible to adjust the radius but with different 
signs for e^^;. In this sense the calculation of other nucleon observables is desirable 
to reach definite conclusions. 

Alternative Approaches 



For the cr, tt, cu, p, Ai-system there are other approaches to the soliton sector 
of the NJL model with vector mesons. Alkofer et al. (1993) propose to make the 
whole numerical calculation in Euclidean space without rotating the single particle 
hamiltonian back to Minkowski space where the problem is originally formulated. 
Instead they use a particular operational prescription to obtain the mass of the 
soliton. The parameters are fixed by means of a heat kernel expansion although 
the masses of the vector mesons are not negligible. These authors find solitons for 
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the whole system. Within the same scheme Zueckert et al. 1993 have also found 
more recently that vector mesons stabilize the soliton of the a, tt, cu, p, Ai-system in 
the linear case. Following their calculations it is possible to fit the nucleon radius 
with a solution having a negative valence quark eigenvalue. Watabe and Toki 
(1992) made calculations in the a, tt, cu-system without continuation of the omega 
field into Euclidean space, which avoids the appearance of an imaginary part of the 
Euclidean action. Besides an explicit breaking of vector gauge symmetry (which 
dissappears in the infinite cutoff limit), their prescription corresponds to have a 
regularized baryon density. In addition they release the chiral circle condition. 
However neither bound nor stable solitons exist in this scheme. 

As we have mentioned at the begining of this section, the introduction of vec- 
tor mesons in the soliton sector is by no means a trivial problem from a conceptual 
point of view. As an unfortunate consequence the various prescriptions proposed 
by different authors quoted at the beginning of sect. 7.3. lead to different results. 
This makes the situation unsatisfactory. In this sense it would be highly inter- 
esting to study this problem more deeply and to see which or even if any of the 
prescriptions suggested is the correct one. 

In particular the problem of regularization has to be addressed since none of 
the above authors really starts from a regularized action as basis for all further 
developments. 
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8. The Nambu— Jona-Lasinio-model and Other Effective Chiral Models 

Finally we want to compare shortly the N JL soliton consisting of valence quarks 
and the polarized Dirac sea with other effective chiral models which have been dis- 
cussed in the literature so far. Both the Skyrme model (sect. 8.1.) and the chiral 
sigma model of Gell-Mann and Levi (sect. 8. 2.) can be formally related to the 
NJL model via gradient and heat-kernel expansions. Whereas the Skyrme model 
contains no explicit valence quark degrees of freedom but relates the baryon num- 
ber to the topological winding number of the Goldstone fields, the valence quarks 
are present in the Gell-Mann-Levi model from the very beginning. As we have 
discussed in chapter 3 both pictures can be incorporated in the NJL model and 
it is therefore interesting to see how these models are formally and numerically 
related. The renormalized chiral sigma model containing polarized Dirac sea as 
well as explicitly the kinetic energy of the mesons shows a translational not in- 
variant ground state and is described in sect. 8. 3. In sect. 8. 4 we consider recent 
extension of the NJL model containing dilaton fields. General aspects concerning 
the connection between those various effective models are given in sect. 8. 5. 

8.1. The Skyrme model 

The Skyrme Model Without Vector Mesons 

The SU{2) Skyrme Lagrangian in the scalar sector without vector mesons 
reads (Skyrme 1961,1962, Adkins et al.l983, for a review cf. Zahed and Brown 
1986, Holzwarth and Schwesinger 1987) 



4- 

+ 



CsKY = \f>r[{^^,u'f)i^^'U) 

^Tr,|[(a^t/tc»,t/)(c»^t/tc»^t/)] - [{d^U^){dm) 



(8.11 



where U{r) = e^^''''''"'' with ^(oo) = denotes the hedgehog chiral field on the 
chiral circle and e is a model parameter, which can be chosen to reproduce certain 
observables. The topological winding number n of the chiral field, which is defined 
by ^(0) = — TiTT, has to be an integer number in order to obtain a finite energy 
^SKY = / dt{—CsKY)- It can be shown that the 4-divergence of the corresponding 
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current, which is the topological or anomalous current (Goldstone and Wilczek 
1981, Witten 1983) 

(-1) 



^ton(^) — 



247r2 



^/iJ/ll^2i^3'JY 



{d^,U){d^,,U){d^,,U)u'f 



1.2) 



vanishes and n is identified with the baryon number of the system. As we have 
already stated (cf. sect. 3. 2, app.B), jf(j„{x) is obtained as the A — th order gradient 
expansion of the 1-loop (sea) part of the baryon current {b^{x)), whereas the 
second order vanishes. If we look at he 4th order gradient or heat-kernel expansion 
of the fermion determinant S^^AU) of the effective action in the NJL (eq. (2.16)) 
itself we find (Aitchison and Frazer 1984,1985a,1985b, Dhar et al.l985, Ebert and 
Reinhardt 1986): 
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(8.4) 
where we have used the proper time scheme for regularizing the fermion determi- 
nant and performed a heat-kernel expansion (cf. sect. 2.5.). T[a,x) denotes the 
incomplete F-function: 

F(a,x) = r dte-^e-"" (8.5) 

Jx 

Indeed, ^cj^y ^^^ ^^^ form of the Skyrme stabilizing term in eq. (8.1) but in 
addition a 'destabilizing' term £^ appears (Aitchison et al.l985). Meissner Th 

et al.(1990) showed that for small constituent quark masses M < 700MeV the 

(4) 
total energy in 4th order heat kernel expansion is indeed unstable, because Ejj 

corresponding to Cjj goes to — oo for meson profiles with small size R (cf. eq. (3.7), 

Fig. 3.2). Due to the behavior of the incomplete F-functions Ejj changes sign for 

larger M and the theory in 4th order heat kernel expansion gets stable. Because 

for very high constituent quark masses the valence quark is part of the negative 
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spectrum one therefore recovers at least the philosophy of the Skyrme model for 
M — ^ oo. On the other hand we know from sect. 3. 3 that the physics of the nucleon 
demands M ^ 400MeV, from which we conclude, that the Skyrme model in the 
Goldstone sector is not a good approximation to the NJL approach. 

The Skyrme Model With Vector Mesons 

In the mesonic sector the heat kernel approximation of the NJL Lagrangian 
has been found to resemble the phenomenological chiral Lagrangian of Gomm et 
al.(1984), where the imaginary part of the effective action of the NJL corresponds 
to the gauged Wess-Zumino term in the Skyrme like theory, in which Gomm et 
al. (1984) adjusted the parameters to reproduce the hadronic processes p ^ tttt 
and Al -^ pn. Using an on shell mass renormalization (cf. sect. 7.2.) the NJL 
prediction for these decays is rather good for constituent quark masses around 300 
and 350MeV. 

Reinhardt and Dang (1989) have claimed that the NJL model with vector 
mesons reduces exactly to the original Skyrme Lagrangian (8.1) if the vector 
mesons are integrated out, and the nonlinear constraint of the chiral circle a + 
TV = f^ is assumed. This indeed is true in a second order heat kernel expan- 
sion. However, in the 4th order additional terms appear, whose magnitude are not 
known. Hence there arises an ambiguity about the proper way to compare with 
the original Skyrme Lagrangian. It is important to notice that the way vector 
mesons are presently described in the Nambu-Jona-Lasinio-model turns out to be 
in full agreement with the old massive Yang-Mills idea (Lee and Nieh 1968) that 
vector mesons are the gauge bosons of chiral symmetry. Such an ideology has been 
widely used within the Skyrme model approach (Meisser UG 1988). Not surpris- 
ingly, the gradient expansion produces a massive Yang-Mills type Lagrangian for 
the non-anomalous part and a gauged Wess-Zumino term for the anomalous part, 
as it is often used in the Skyrme model (Zahed and Meissner UG 1986, Meissner 
UG and Zahed 1987). 

Finally we want to state that up to now it remains an open question, if the 
solution of the NJL including u, p and Ai mesons as degrees of freedom supports 
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a bosonized, Skyrme like picture and moreover if the explicit inclusion of vector 
mesons in the NJL is reasonable at all (cf. sect. 7.3.). In fact it is known that 
the way they are introduced in the Nambu-Jona-Lasinio-model through direct 
vector and axial couplings breaks the proper chiral anomaly (Wakamatsu 1988, 
Ruiz Arriola and Salcedo 1993a). 

8.2. The Chiral Quark Meson Model without Vector Mesons 

The Chiral Sigma model (CSM) is based on a Lagrangian of Gell-Mann and 
Levi (1960) and has been quite successful in the last years for a detailed description 
of nucleon and A observables and form factors. The Lagrangian in SU{2) and for a 
and 7? fields reads in the chiral limit (rMjr = 0) and on the chiral circle a +7i = /^: 

C = qi^q - gq{a + iT^-i^)q + - (d^'adna + d^ndn^) (8.6) 



In practical calculations (Birse and Banerjee 1984, Birse 1985, Goeke et al 1985, 
Cohen and Broniowski 1986, Fiolhais et al.l987, Fiolhais et al.l988, Alberto et 
al.l988, Meissner Th et al.l989, Birse 1990, Neuber and Goeke 1992) the Dirac 
sea has been ignored and the q was taken to be a spinor of the Nc = 3 valence 
quarks. The a and n are first of all classical static fields. In this approximation 
the total energy of a soliton reads: 

EcSM = Nce.al + \jd^^ [(Va)2 + (W)^] (8.7) 

where e^^; is the valence quark energy. The e^^; and a{x) and ti{x) are obtained 
easily from the equations of motion: 

^EcSM _ Q 



5q{x) 
da{x) 
6tx{x) 



(8.8) 





in the hedgehog approximation. Solving these equations yields selfconsistent va- 
lence quark, a and tt fields of the Lagrangian eq. (8.6). Observables can be cal- 
culated in the usual way by coupling of Cq^j^ to external currents and inserting 
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the selfconsistent fields into the corresponding expressions, where the quark and 
meson fields have to be properly quantized, so that the system carries good spin, 
isospin or momentum quantum numbers, respectively. This quantization has been 
done in the semiclassical cranking way (Cohen and Broniowski 1986). Another 
way consists in assuming coherent Fock states for the sigma and pion field and 
to use Peierls-Yoccoz projection techniques (Birse 1985, Fiolhais et al. 1987,1988, 
Alberto et al. 1988, Neuber and Goeke 1992). 

Actually the E(jgf^ of eq. (8.7) can be directly obtained in the gradient 
or heat-kernel expansion of the fermion determinant S^rAa^ii) of the effective 
NJL action eq. (2.16) (Meissner Th et al. 1988, 1990). If one truncates the ex- 
pansion after the first nonvanishing order one obtains directly the kinetic energy 
^jd X (Vcr) + (Vtt) of the meson fields (cf. sect. 2.5.), which, by adding 
^c^vab yields E(jg]\^. This means that the kinetic terms of the mesons in the 
CSM have been generated from the polarized Dirac sea of the NJL by gradient 
expansion. A comparison between the two models will be performed now very 
easily: We once solve the NJL model in the 1 quark loop approximation and eval- 
uate then observables. This is called NJL (exact) in Tab. 8.1. Then we perform 
a gradient expansion of the energy and several observables up to the first non- 
vanishing order. We insert the selfconsistent NJL solutions into those gradient 
expanded expressions yielding NJL (grad.). These numbers are compared with 
the selfconsistent solutions and the corresponding observables from the CSM. Us- 
ing in all three cases the semiclassical quantization the results of some relevant 
observables are compared in Tab. 8.1, where the contributions from the Dirac sea 
and of the meson cloud are treated on a equal footing. Apparently the total val- 
ues of the observables agree within 15%, while the separated contributions from 
valence quarks and mesons (sea quarks) differ. The gj^ is a clear exception, there 
the gradient or heat kernel expansion seems to converge badly. The gj^ from the 
Dirac sea of NJL basically vanishes, while the mesonic contribution in the CSM is 
comparable with the valence quark contribution. Thus a general statement on a 
comparison between these models is impossible and can only be done separately 
for each observable. 
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Finally we want to note that also the full linear version of the CSM, which 

2 

contains in addition to (8.6) the mexican hat potential reading — n^(o" +7f — /^) 
(if m-n; = 0) can be obtained from the NJL by gradient or heat kernel expansion, 
where the a-mass conies out to be m^ = 4:(gfj^) . In contrast to the NJL, where 
without constraining the meson fields to the chiral circle the soliton collapses (cf. 
sect. 3.5), the mexican hat potential arising in second order gradient expansion 
is highly stabilizing and therefore the linear CSM shows perfectly stable solitonic 
solutions. Moreover it turns out that the nucleon observables calculated from those 
solutions in the CSM do not differ very much from the ones which are obtained in 
the nonlinear version. 

This fact and the behavior of gj\ as described above show that the gradient 
expansion cannot be regarded as unreliable in the solitonic sector as a generally 
reliable approximation. It has to be used with care. 

8.3. The Chiral Quark Meson Model without Vector Mesons 

The chiral quark meson model with vector mesons proposed by Broniowski and 
Banerjee (1985,1986) is based on the massive Yang-Mills Lagrangian of Lee and 
Nieh (1968) supplemented with valence quarks coupled to all mesons {awivp a). 
Many practical calculations (Broniowski and Banerjee 1985,1986; Broniowski and 
Cohen 1986; Ruiz Arriola et al. 1990; Alberto et al. 1990b; Ruiz Arriola et al. 
1993b) neglect the effects of the Dirac sea. Similiarly to the case without vector 
mesons, they may be considered to be included in an approximate way in the 
kinetic energy of the mesons. It should be also mentioned that this type of models 
do not include a gauged Wess-Zumino term as it is done in Skyrme type models. 
From the point of view of the Nambu-Jona-Lasinio-model this is not a problem 
since the gauged Wess-Zumino terms correspond to the valence quarks (Schueren 
et al. 1993) (see also the discussions in sect. 7.3.) The dynamical question, 
whether a chiral quark model with vector mesons is supported by the NJL model, 
remains an open problem. 
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8.4. The Renormalized Chiral Sigma Model - Vacuum Instability 

Some authors (Soni 1987, Ripka and Kahana 1987, Li et al.l989) considered 
a renormalized CSM, i.e. they solved the Lagrangian (8.6) in the 1-quark loop 
approximation and subtracted the divergent parts as local counterterms. Such an 
approach does not contain an UV cutoff by construction. In our language the total 
energy of the renormalized CSM on the chiral circle reads: 

E'JSM = Nce^al + hld\ [(Va)2 + (VTr)^" 

A /• (8-9) 

+ hm {i?,ea(A) " ^^9^1^)1 / d\ \{Vaf + (W)^ } 

where /2 (M) is the logarithmically divergent integral from eq. (2.28). Due to 
the subtraction of the counterterm {AN^g l2)\ j d r (Va) + (Vtt) from the 
divergent sea energy Egea the whole expression remains finite in the limit A ^ oo. 



^2 , (^r?\2 



In the renormalized CSM the kinetic energy of the mesons \ J d r ( Va) + (Vtt 
is added explicitly rather than generated by the polarized Dirac sea as in the NJL 
where the cutoff' A is properly adjusted. Actually one can easily see, that the 
total energy of the NJL and the renormalized CSM formally only differ by the 
UV convergent terms. In the NJL those are evaluated with a finite UV cutoff A 
reproducing the pion decay constant, whereas in the renormalized CSM they are 
treated with infinite A. Though this small difference has a tremendous effect. The 
authors quoted above showed that the sea energy Egea gets negative for meson 
profiles with finite size R. This means that the full vacuum is not translationally 
invariant. Moreover it turns out that Egea gets more and more negative, if one 
increases the topological winding number n of the pion field and therefore neither 
a stable vacuum nor a soliton solution with finite baryon number B exists. 

Recently it has been suggested that this problem does not occur if dynamical 
vector mesons are included in the CSM (Kahana and Ripka 1992) in a massive 
Yang-Mills scheme and considering the corresponding quantum corrections. 
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8.5. Scale Invariance and Dilaton Fields 

Ripka and Jaminon (1992) extended the NJL to a model which in addition to 
the spontaneously broken chiral symmetry also exhibits the anomalous breaking of 
scale invariance in QCD. Following the work of Jain et al. (1987) they introduced 
to this end a scalar-isoscalar dilaton field X, which is coupled to the effective action 
by: 



oo 



2 



1/(A2X2) 

(8.10) 
where a is a dimensionless parameter. The first two terms in eq. (8.10) are 
scale invariant, the scale breaking being related to the X-field Lagrangian C{X). 
One chooses C{X) in such a way that the divergence of the dilation current s^ 
is d^s^^ = X and therefore the vacuum value Xy can be related to the gluon 
condensate (Gj^j^). It has been shown (Jain et al.l987, Ripka and Jaminon 1992) 
that, when a vacuum value Xy ~ 350MeV is used, as deduced from the QCD sum 
rule estimate (Shifman et al.l979), the X field stays almost constant throughout 
the chiral phase transition, in which the chiral symmetry is restored at high baryon 
densities or high temperatures. 

Motivated by this fact Meissner Th et al.(1993) and Weiss et al.(1993) used a 
constant dilaton field X = Xy by solving the theory (8.10) in the solitonic sector. 
The scale invariant form of the mesonic interaction in eq. (8.10), which differs from 
the mesonic mass term in the original NJL (2.16) 

,2 . 

/ d^x{a^ + 7f2 - af ) 



2 



has the effect, that solitonic solutions of eq. (8.10) exist even in the case, when 
both a and tt are not restricted to the chiral circle a + n = ay (non-linear 
model), but both a and yf degrees of freedom are fully allowed (linear model). 
This was not the case in the original version of the NJL, where the constraint of 
the chiral circle turned out to be necessary, because otherwise the soliton collapses 
to a configuration with zero size and zero energy but baryon number B = 1 (cf. 
sect. 3.5). One can easily convince oneself, that, with the notation of sect. 3. 5, the 
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meson energy in eq. (8.10) gets now proportional to U R oc R . The soliton 
energy can then only fall to zero as R ^ when a < |. However, as it was shown 
in sect. 3. 5, the soliton can only maintain a baryon number B = 1 while its total 
energy tends to zero if ct > 1. For this reason the B = 1 soliton does not collapse 
when it is calculated with the scale invariant action (8.10). 

Furthermore Meissner Th et al.(1993) and Weiss et al.(1993) have shown that 
for the physical relevant region of the constituent mass M ^ 400MeV the devia- 
tion of a and n obtained as solitonic solutions of (8.10) from the chiral circle is 
numerically very small. Moreover all calculated observables are nearly identical to 
the ones which are obtained from the original version of the NJL, where a and tt 
are restricted to the chiral circle from the very beginning. 

The nonlinear version of the NJL, as it has always been used in the past and 
on which the present review is based on, can therefore be justified by a model 
which implements the anomalous breaking of scale invariance in QCD. 

8.6. General Aspects 

If one looks from a bird's view on section 8.1 ,8.2 ,6.4. and 7.3., the following 
qualitative evaluation can be done: The simple Skyrme model with a pion field 
only cannot be compared with NJL. In the heat kernel expansion a destabilizing 
term occurs and it also turns out that the result for semiclassical quantized ob- 
servables are different. The inclusion of vector mesons in the Skyrme Lagrangian, 
however, allows a better comparison with the NJL model in the scalar and pseu- 
doscalar sector. The observables are similar and in the collective hamiltonian of 
SU{3) structurally the same terms appear and the numerical results are close as 
well. Thus it seems that the role of the valence quarks in the NJL is qualitatively 
played by the vector mesons in the Skyrme model. Presently one cannot judge 
the role of vector mesons in NJL. Altogether it seems, and this includes the chiral 
sigma model as well, that the degrees of freedom relevant for the low energy regime 
of QCD can effectively be parameterized in rather different ways still reproducing 
the low energy hadronic phenomena with similar accuracy. One, however, is quite 
clear: 
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The Goldstone boson fields or Goldstone bosonic quark-quark interactions are play- 
ing the dominant role. 
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9. Summary 

The present review article deals with effective chiral symmetric models of 
quarks and non-dynamical mesons for the description of mesonic and baryonic 
ground states. In general effective theories are trying to model low energy QCD 
phenomena by identifying the appropriate degrees of freedom and the relevant 
symmetries neccessary for a qualitative and if possible, quantitative description 
of low energy hadronic phenomena. Experience over the last 20 years has cum- 
mulated in the common opinion that for baryonic ground states spontaneously 
broken chiral symmetry and the incorporation of the corresponding Goldstone- 
boson fields are the important ingredients for a proper effective approach. 

The review concentrates on the use of the Nambu-Jona-Lasinio-model and re- 
lated approaches in the sector with baryon number B = 1. In its basic form it is 
defined by being the simplest pure quark theory with local biquadratic interactions 
of Goldstone character exhibiting spontaneous breakdown of chiral symmetry. One 
should note that the Nambu-Jona-Lasinio-model plays a central role because it 
allows to build conceptual and numerical bridges to other basic approaches, such 
as fully bosonized approaches of Skyrme type and theories, where valence quarks 
are coupled to dynamical meson fields. 

The Nambu-Jona-Lasinio-model is treated explicitly in the path integral for- 
malism performing a suitable bosonization procedure. Then the theory is described 
in terms of quarks and composite, non-dynamical meson fields. Actually the model 
is solved in the zero-boson and one-quark-loop approximation, such that it can be 
uniquely defined by a so called fermion determinant. In this form and neglecting 
the scalar fiuctuations around the vacuum it reduces to the solutions of the chiral 
quark loop and agrees with the structure of the Chiral Quark Model. Performing 
an expansion of the determinant under the assumption of slowly varying pionic 
fields gives the structure of the Skyrme model, including the whole anomalous 
sector being described by the Wess-Zumino term. 

Actually in the present model and for time-independent chiral fields on the 
chiral circle the fermion determinant is evaluated exactly and provides a full treat- 
ment of the polarization of the Dirac sea. The latter is caused by the presence 
of a discrete and localized valence level within the single-particle spectrum. The 
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procedure used in the literature consists of three steps: First, a regularization 
scheme is chosen and the parameters of the model are fixed in the mesonic sector 
to reproduce PCAC, meson masses and decay constants. Second it is checked 
that the vacuum condensates and the current quark masses are reasonable. Third, 
without changing the parameters solitonic solutions are selfconsistently obtained 
in the sector with baryon number B = 1. 

In order to describe baryons with the quantum numbers of the spin 1/2 and 
spin 3/2 multiplets, the hedgehog based quark field ansatz has to be quantized in 
the collective subspace. Therefore a time dependent rotation in the direction of 
the symmetry of the model is performed and canonical commutation relations are 
imposed on the collective coordinates and the canonical momenta. This results in 
proper commutation relations for the generators of the group, which in the case 
of SU(2) are the components of the total angular momentum and in SU(3) are 
complemented by some generators acting on the spin-baryon number space. The 
hedgehog ansatz for the SU(2) fields, which manifests the relation between spatial 
and isospin rotations, results after the quantization in the constraint that the ab- 
solute value of spin always equals isospin. However the value of the spin itself can 
be integer or half integer and therefore the soliton can be quantized as fermion as 
well as boson. Only the generalization to SU(3) gives a constraint to fermionic 
solitons as long as the number of colors is odd. 

The numerical results for the solitonic sector of the Nambu-Jona-Lasinio- 
model are basically independent on the regularization scheme used as long as 
the cut-off is treated as a parameter of the system, which is fixed in the meson 
sector. The calculations of the literature, all being based on hedgehog structure 
and physical values for the meson sector, can be summarized as follows (status 
autumn 1993): 

i) SU(2), a and n field, chiral circle: For a constituent quark mass M = 
420MeV the following observables within 15% are reproduced. Isoscalar and 
isovector charge squared radius of the nucleon, < r >x^q ^^'^ < ^ ^r==l' 
nucleon sigma term, E, nucleon axial coupling constant and magnetic mo- 
ments of proton and neutron (if the recent I/Nq corrections are included), 
(7^, Up and Hn, the nucleon-delta splitting, M^ — -^tVj the q-dependence 
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of the form factor G^^{q^), G^j^{q'^), G'^{q'^), QAio^)- The nucleon energy- 
comes out at about SOOMeV if rotational and translational zero modes cor- 
rections are included. The neutron squared radius, being a very sensitive 
quantity, is twice as large as the experimental value, the G^{q ) is generally 
by a factor of two too large at finite q . 

ii) SU(3), a, tt, K and r] fields, chiral circle for sigma and pion, trivial em- 
bedding of SU(2) into SU(3), perturbative treatment of iris up to second 
order or Yabu-Ando approach: Splitting between and within spin 1/2 and 
spin 3/2 baryons is reproduced within few MeV. The splitting within isospin 
multiplets is reproduced and a common value for rriu — rrid is found. If the 
zero point corrections to translation and the SU(3)-rotations are included 
the energies of all octet and decuplet baryons are too large by a constant 
shift of approximately 100 MeV. If recent 1/Nc corrections are included the 
9a^ 9 a ^"^^ 9 A °^ ^^^ nucleon are reproduced slightly outside the experimen- 
tal errors. The nucleon E-term is reproduced with a strangeness content of 
15% and the recently measured Gottfried sum comes out as well. 

iii) SU(2), o", TT, p, Ai and uj mesons, chiral circle for sigma and pion. None of 
the present approaches has been developed to an extent that one can judge 
the infiuences or the necessity of vector mesonic coulings. The observables 
calculated by now are the isoscalar charge radius of the nucleon and (7^, both 
deviating noticebly from experiment. There seems to be also conceptual 
problems with treating the omega meson. On the other hand the repulsive 
character of the omega guarantees a stable soliton and one does not require 
any more the restriction to the chiral circle. 

The problem encountered by now can be summarized as following: The semi- 
classical quantization procedure is not perfectly understood yet, although it is 
common usage. Its relationship to well-established quantum-mechanical theories 
like Peierls-Yoccoz projection approaches is yet unknown. If one looks at g/^ and 
magnetic moments of the nucleon its convergence in the rotational frequency O 
seems to be slow. Furthermore it has been shown that the PCAC relation is re- 
produced within the cranking scheme if rotational corrections in linear order of O 
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are included. 

The quantum corrections due to zero modes of translation and rotation add 
up to 30 — 40% of the resulting baryon mass and hence are by far too large for a 
correction. The corresponding corrections of non-zero modes are of order 0{Nc) 
as well and they are not known yet. There are indications in the Skyrme model 
that they are small, but a clear statement in the Nambu-Jona-Lasinio-model is 
missing. 

The treatment of the w-meson is not quite settled yet. A consistent scheme, 
which start from the very beginning with an action, the real part of which is reg- 
ularized in Euclidean space, is still missing. Furthermore one needs some more 
observables to be evaluated. Hence the question if vector mesons in such a quark 
theory are indeed necessary to be included as quark-quark couplings is not yet 
settled. 

Altogether one can summarize: Quark models like the Nambu-Jona-Lasinio- 
model with a polarized Dirac sea and involving Goldstone type couplings (i.e. 
non-dynamic mesons), have been by now quite successful in describing the ground 
states of baryons in the octet and decuplet. 

If one ignores vector mesons, whose relevance in the present model is still under 
debate, the calculations support a clear picture of the baryons consisting of three 
localized valence quarks interacting with a moderately polarized Dirac sea. This is 
in contrast to Skyrme-like models, which are governed by meson fields, exhibiting 
certain topological properties, and quark-meson models, where valence quarks are 
coupled to dynamical meson fields. 

The problem encountered concern the semiclassical quantization method. It 
has certain merits, but it deserves improvement. Last but not least, the NJL- 
model lacks confinement which should be necessary for orbitally excited states of 
baryons. 
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Appendix A. Minkowski- and Euclidean Space-Time: Notation and 
Convention 

Tab. A.l compares the notation for 4- vectors A^, the metric tensor gr^jy, the 
scalar product between two 4-vectors as well as the Dirac 7-matrices in Minkowski 
and Euclidean space. 

Generally a Minkowski 4-vector is transformed to an Euclidean 4-vector by 
performing a Wick rotation: 

Ae^ = {+t)AM^ 

AE^ = -AE'^ = {-i)AM^ = {-i)AMo ^^ ^ 

Ae' = Am' 

Asi = Ami 

Hereby both Ae and A^ are assumed to be real numbers, which means that any 
function /(A^^) depending on the real variable Aj^ g^ts analytically continued 
to the complex plane. For the 4-space vector and the 4-momentum vector we have 
especially: 

XM^ = {t, rl XE^ = (r, r) 

r = {+i)t J d xe = i+i) J d XM d^^^^^>dl 



and 



PM^ = {^^P) PE^ = {u,P) 



u = i+i) V J d'^pE = (+^) J d^M g^ = (-^)^ 
with the notations: 

t: Minkowski time r: Euclidean time 

u: Minkowski frequency u: Euclidean frequency 

Further examples handled in this way are the cranking frequency O (cf. chap. 
3), which couples like the time component of an isovector vector meson (p meson), 
or the LJ meson, which is time component of an isoscalar vector meson (cf. chap. 
7). 
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The transformation of a scalar product reads: 

(Am ■ Bm) = 

^m^Bm^ = Am Amo + Aj^Ami = Am Am - Am^Am^ 
- Ae^Ae^ - Ae'Ae' = Ae^Ae^ + Ae'Ae, = Ae^Be^ = 
{Ae ■ Be) 



I.e. 



{Am ■ Bm) = {Ae ■ Be) (AA) 

Furthermore we notate: 

(Am)^ = {Am ■ Am) 

{Ae)^ = Ae^'Ae^ = -{Ae-Ae) 
giving: 



(A.5) 



(Am)^ = -(Ae)^ iA.6) 



Fig. A A shows the contours in the complex plane arising for the calculation of 
static observables of the soliton in Minkowski space: 

27r Y « - CA 
and in Euclidean space: 

2n^ tu-ex 
if the one neglects regularization (cf.eq. (3.16)). 

In order to obtain for a certain expression the same results in both cases one 
has to make the central assumption, that during rotating the contour-line from 
Minkowski into Euclidean space no singularity in the complex plane is hit. This 
is clearly true for the non-interacting case, but has to be postulated in general 

{Euclidean Postulate). 

A system with baryon number B = 1 can be constructed either by extending 
the corresponding contours as shown in Fig. A A (dashed lines) or by introducing 
a thermochemical potential //, which causes a shift: 

ex^ex-^i {A.9) 
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in the single particle orbitals and therefore increases the number of singularities 
in the closed contour line by 1, namely by the valence orbit (cf. sect. 3. 2). It is 
clear that this prescription only works, if ^ is both in the Minkowski space and 
the Euclidean space is a real number and therefore no analytic continuation (Wick 
rotation), as it was posted in eq. (A.l), may be performed. The treatment of the 
thermochemical potential /u, in Euclidean space is therefore different from that of 
the time component of a physical 4-vector, as it is the lj meson. 

Also the Dirac 7-matrices can be formally transformed into Euclidean space: 

1E^ = (+^)7M° 1e' = 1m' (^-10) 

Though the 7-matrices are no physical objects, which transform under Lorentz- 
transformations as a Lorentz 4-vector. The eq. (A. 10) is therefore only a simple 
redefinition, which means, that 7^; is antihermitian as it is 7*. The reason for 
doing this is just because the rule for transforming a scalar product {AA) can now 
also be applied to the expression yl: 



^M = 4^E 



(All) 



For the matrices (3 and 75 we will not define Euclidean values (as it is sometimes 
done especially for 75), but always use the same matrices, namely: 



and 



/9 = 7M =H)7£;^ (A12) 



75 = 7^ = (+«)7m°7mSm^7m^ = lE^lE^lE^lE^ {A.l?,) 



Generally we will use the Bjorken-Drell notation (1965): 



„_/l2 \ ^,(0a' 



-« -U _ij -«-U.' oj <^-"' 



with the Pauli matrices: 



^-Ci) ^<:) Mil.) <-' 
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and especially: 

as well as: 

7* = Pa' {A.17) 



Appendix B. Gradient Expansion of the Baryon Current 

Because of the importance for our model we show explicitly the gradient expan- 
sion of the baryon current b^{x) and especially the baryon number B = J d xbQ{x), 
as it has been performed by Goldstone and Wilczek (1981) and discussed in detail 
by Witten (1983) in context of the topological quantization of B. 

From the general treatment given in sect. 3. 2 we have for the expectation value 
of 6^: 

(b^'ix)) = —{x \Tr^Xc [(«^)~S'i I x) - vac.contr. (S.l) 

Hereby the U denotes a SU{n) chiral field which is parameterized in terms of the 
Goldstone field ^(x) 

U{x) = e*^(^) (S.2) 

$,{x) is hermitian and traceless, i.e. 

ax) = E eix)- (S.3) 

a—l 

where the A*^ being the generators of SU{n). iD is the quark Dirac operator with 
the chiral field U: 






-1 .. .. (5-4) 



with: 



t^5 = ^(t/ + t/t) + l75(t/-t/t) (S.5) 
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We therefore find: 

{x \Tt^Xc [l^idr^d"" + M^)-\i^ + M)] I x)} 
In this expression we can now perform the power expansion: 



x)- 



(S.6) 



{A + B 



| — 1 _ ^— 1 /| — 1d/|— 1 



A-^BA-^ + ... = A-^ J2 i-nBA 

171=0 



— Into 



(S.7) 



with: 



A := {d^d" + M^) 
B := i^U^ 



G- 



{B.8) 



Because of Tr;^^ = we find that aU terms in eq. (-B.8) up to order m — 2 vanish. 
For m = 3 we get: 

{b^'{x)) - {b^'{x))v = i^M'^Tt^xi^ \ [7^G(^t/5)G(^t/5)G(^t/5)Gt/5t] \x) + ... 

(S.9) 
The G is diagonal in the momentum space basis | pe)- In the lowest order gradient 
expansion the ^U^ does not vary with the coordinate x (cf. sect. 2. 5). After 
inserting complete sets of | x^) and \pe)j respectively, one obtains: 



mx)) - r{x))v 

,3^4 / d^PE 



Tr 



7A 



7 



and finally after performing the 7 trace 



^(^t/5(x))(^t/5(x))(^t/5(x))t/5t(x) 



+ ... 
(S.IO) 



{bf'ix)) - mx))v 



(-1) 

247r2' 



fiiyiiy2V3r^j. 



X 



{d^,U){d„^U){d„^U)U^ (S.ll) 



which is the anomalous baryon current of Goldstone and Wilczek (1980). 

The same result can be obtained, if one starts from the imaginary part of the 
Euclidean effective action and performs a gradient expansion (cf.sect.2.4), which 
leads to the famous Wess-Zumino-Witten term (2.37). The algebraic manipulation 
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for doing this is very similar to the one, which was used above for gradient ex- 
panding the baryon current itself. 't/(l)-gauging' the Wess-Zumino-Witten term 
(2.37), as it was shown by Witten (1983), gives also the expression in eq. (-B.11) 
as conserved current. 

In the special case of the 5't/ (2) -hedgehog form for U: 

U{r) = cos e{r) + i{Tr) sin e{r) (S.12) 

one finds for the time component of eq. (B.ll): 



1 dOsin^e 
and after integration J d r and using the fact, that 6{oo) = 0: 



{b^{x)) - {b^ix))v = ir^^-^ + ■■■ (5-13) 

2n^ dr r^ 



27r n 

For ^(0) = —mr one recognizes that in case of integer n the gradient expanded 
baryon number (B) coincides with n and therefore the gradient expansion is valid 
exactly. On the other hand for non integer n the gradient expanded expression 
(S.14) is also non integer in contrast to the exactly calculated (B) and therefore 
the gradient expansion does not converge in this case. 

A necessary condition for the validity of the expansion (-B.7) as well as the 
neglegibility of the space dependence of ^U^ in (-B.9) is: 

\dU\<^M (S.15) 

which has to hold at any point r in the radial integral {B) = J d r{bQ{r)). If the 
meson profile U{f) is characterized by some size parameter R (as e.g. in Fig. 3.1, 3.2) 
one finds therefore: 

|n|7r <M-i? (S.16) 

which means that the gradient expansion is good either for large profile sizes or 
also in case of selfconsistent solutions for high constituent masses M. 
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Appendix C. Lorentz Transformation of the Mean Field Soliton 
Energy 

In this appendix we want to show explicitly that the soliton mean field energy 
behaves under Lorentz boost transformations like the time component of a Lorentz 
4- vector with vanishing space components, and therefore prove formula eq. (4.43), 
which is the main ingredient in the pushing approach. 



Let us to this end consider a Lorentz boost in xi-direction with velocity v. 



X' 



f^^x'^ = Ai^V) 



X 



(C.l) 



where: 

uj = arctanh(t;) ■ v 

v= 
and the contra- and covariant Lorentz transformation read: 



A^.( 



UJ 



K'^'i^) 



( cosh u — sinh u \ 

— sinhcu cosh a; 

10 

V 1/ 
^ cosh a; sinh a; 0^ 

sinh a; cosh a; 

10 

V 1/ 



respectively. The generator of this transformation is the unitary matrix: 



(C.2) 



(C.3) 



E{Q) = e-^^ 



(C.4) 



where the infinitesimal boost generators K^ are defined by: 



K' 



t d j . d .a* 

-7-^ — X 1-— + I — 
i ox^ ot 2 



(C.5) 
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One can easily convince oneself that: 

and therefore derive the action of B{uj) on the Dirac operator iD: 

B^'^iu) Z7^a^-5f((j(x) +zf7f(x)75) B{u)= i'f^d^ - g (^a{x') + iTn{x')'y5^ 

(C.7) 

In the 1-quark loop approximation the fermionic part of the energy of a system 
with meson field configuration {a{x)jw{x)} is defined as the expectation value of 
the 1-particle hamiltonian h = ^j- + gl3[a{x) + zf7f(x)75] due to eq. (3.14): 

J4, 



dt{h{x)) 



VqVqJd'^x [q{x)(3hq{x)] e*/'^ XEq{x){tD[a(x),7r{x)])q{a 



VqVqe^ / d'^xq{x){iD[a{x),Tr{x)])q{x) 



d 



(as) 



« oYSpln[«9t - h{x) - T/i(x)]x=o 



In case of time independent fields a{x) and 7t{x) this reduces to 

/ dt{h),tat = (-^r) / ^ E -^ = ^ E ^A = r^;^ 

J JCb ^TT ^ u- ex ^^^^ 



iC.9) 



which is of course the same as the unregularized expression obtained in Euclidean 
space (eq. (3.116)), if the integration contour C^ is chosen in such a way that 
the appropriate number of orbitals necessary for obtaining a system with baryon 
number B gets occupied. In Fig. A.l this has been indicated for the cases 5 = 
and B = 1. 

On the other hand in the boosted system (') we have the expectation value of 
the boosted hamiltonian hcj{x) = ^ + g/3[a{x'{x)) + if7f{x'{x))'y^] 

dt{haj{x (x))) = 



T=0 



d r 

« oy Spin [i'j^d^ - g{a{x'{x)) + zf7f(x'(x))75) - T/3huj{x'{x)) 

i^Spln {B-^-CJ) [z7^a^ - g{a{x'{x)) + irnix'ix)^^) - T/3/i^(x'(x))] B{-CJ)}^^^ 

(C.IO) 
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Applying eq. (C.6) one finds (Betz and Goldfiam 1983): 



dt{huj{x'{x))) = T{u)E{u) = 

dt{—idt + h{x)) 
+ (Tcosh(a;)) ■ [cosh(a;)(/i(x)) — sinh(a;)((aa))/i(x)) 
+ sinh(a;)((pa))/i(x)) — tanh(a;) sinh(a;)((acD)(pa)))] 



(ail) 



For time independent meson field it holds: 

{-idt + h{x)) =0 

(ah) = (p) =0 
Furthermore for fixed index i (e.g.z = 1 in case of eq. (C.2)) one has: 

i^iPi) = {xi ■ [gl3d^{a{x) + iTn{x)'yr^)]) (CIS) 

which can easily be shown to vanish if the meson fields fulfill the classical mean 
field equations of motion (3.51) refiecting a virial theorem (Rafelski 1977). 

Therefore we end up finally with: 

f dt{hu;{x'{x))) = T{u)E{u) = {Tcosh{uj)){EMFCosh{uj)) (C.U) 

Because the time interval T as itself is boosted to 7^ = cosh(a;)T, we obtain the 
desired relation (4.43) between the energy in the boosted (moving) system i?^ and 
the static mean field energy E^p: 

E{u) = cosh{u)EMF (C-15) 
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Appendix D. Quantization condition for Nc and why solitons are 
fermions 

We know from the action of right generators Rj^ of Sect. 6 on the baryonic 
wave- functions '^{A), that (cf. Balachandran et al. 1984) 

expieai?a^(^) = ^(Aexp(iea^Aa)) {D.l) 

For the special case of a rotation around the 8t/i axis, we have to take into account 
that i?8 = NcB/{2\^) is constraint. Using the fact exp (i27r\/3A8) = 1, which 
corresponds to ©g = 4:n^/3 in eq. (-D.1), one obtains 

expi47T^BNc^{A) = ^{A) {D.2) 

so that B X Nc is integer. The number of colors therefore is an integer quantity. 
This corresponds to the quantization condition for electric and magnetic charges, 
when we interpret the Lagrangean L^ as describing the moment of a charged 
particle in a monopole gauge field (Jackiw 1983). But we can deduce another 
quantization result. Under a rotation of 2tc around the 3th axis, we observe that 

exp (iOshh) = ^^P (^©S^Ag) {D.3) 

for Gg = 27rv^ and G3 = 2tt. Using eq. (-D.1) again one gets on the one hand 

exp (zegi?g)^(A) = exp (z7r)^(A) = -^(A) (DA) 

and on the other hand using eq. (-D.3) for Nc = 3 

exp(iegi?g)vl/(A) = vl/(Aexp(ie3^A3)) {D.5) 

Comparison of eqs. {DA, D.5) gives that the wave functions "^{A) changes sign 
under a space rotation of 27r. Therefore in the actual case of Nc = 3, and in fact 
for an arbitrary odd value of the number of colors, we have quantized the soliton 
as a fermion (and as a boson in the case of an even number of colors). It was 
noted already by Finkelstein and Rubinstein (1968) and Witten (1983b), that in 
the case of SU(2), where the corresponding Wess-Zumino action - which leads to 
the constraint for the eighth generator - is vanishing, the soliton can be quantized 
as fermion as well as boson. 
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FIGURE CAPTIONS 

1.1. Free Dirac spectrum for current quarks {chiral symmetry spontanously 
unbroken) 

1.2. Free Dirac spectrum for constituent quarks with mass M ( chiral symme- 
try spontanously broken) 

1.3. Dirac spectrum for the sohton. The bound valence quarks e^;^; > 
polarize the Dirac sea weakly [valence quark picture). 

1.4. Dirac spectrum for the soliton. The bound valence quarks e^^^ < 
polarize the Dirac sea strongly [fully bosonized picture). 

2.1. The cutoff A for the the non-covariant (0(3)), covariant (0(4)) and single 
(PT) step proper time regularization as a function of the constituent 
quark mass. 

3.1. The bound state spectrum of the 1-particle hamiltonian h for an expo- 
nential profile form 9[r) = —mie~^ with n = 1. e;^ and R are given 
in scaled units [^ and R • M). The valence orbit (0^") comes down 
from the positive into the negative spectrum (Kahana and Ripka 1984). 
Furthermore a value for the thermochemical potential ^ is shown. 

3.2. Sea energy Egea^ valence energy e^^; and total energy Etot as well as 
the corresponding baryon numbers Bgea^ By^i and Bfot = Bgea + -S^a^ 
respectively, for a linear profile form 



r 

e[r) 

10 iir> R 



-n7r(l-^) iir<R 



with n = 1 at a constituent mass of M = 565MeV. The total energy 
Etot[R) shows a local minimum dX R = 0.8fm. B^^j^i and Bsea both 
jump if e^g^i gets negative [R ~ 1.3fm), whereas Btot = 1 in any case. 
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(Diakonov et al.l988, Meissner Th et al. 1989) 

3.3. The mean field energy Ej^p of tfie selfconsistent solution in the non-linear 
model (tott = 0) (full line) splitted in valence (dashed-dotted) and sea 
(dotted) part in dependence of the constituent quark mass M. (Meissner 
Th et al.l989) 

3.4. The quadratic isoscalar electric radius < Br >=< R >p + < R >n of 
the selfconsistent solution in the non-linear model (tott = 0) (full line) 
splitted in valence (dashed-dotted) and sea (dotted) part in dependence 
of the constituent quark mass M. (Meissner Th et al.l989) 

3.5. The selfconsistent a field (normalized to f-n-) in the non-linear model 
with TOtt = for 3 different constituent quark masses M (Meissner Th 
and Goeke 1991) . 

3.6. The selfconsistent n field (normalized to /t^) in the non-linear model 
with m-jr = for 3 different constituent quark masses M (Meissner Th 
and Goeke 1991) . 

3.7. The baryon densities bQ{r) of the selfconsistent solutions in the non-linear 
model with m-n- = for 2 different constituent quark masses M (full 
lines). The valence contributions (dashed-dotted) are explicitly shown as 
well (Meissner Th and Goeke 1991). The upper two curves correspond 
to M == 725MeV and the lower two to M = 363MeV. 

4.1. The rotational moment of inertia (right y-axis), its valence and sea-quark 
contribution, and the nucleon-delta-splitting (left y-axis) is given in de- 
pendence of the constituent quark mass M. The pion mass is chosen to 
be 139MeV. 

4.2. The mean field energy Ej^jp (HEDGEHOG) as well as the energies of 
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nucleon E^q (NUCLEON) and delta ^a (DELTA) and the translational 
and rotational corrections of the selfconsistent solution in the non-linear 
model (tott = 139MeV) in dependence of the constituent quark mass M 
(Pobylitsa et al. 1992) . 

5.1. The quadratic radii of the proton and the neutron are given in dependence 
of the constituent quark mass M. The experimental values are indicated. 

5.2. The proton electric form factor Gj^{Q) for 4 different values of the con- 
stituent mass M (Gorski et al.l992). 

5.3. The neutron electric form factor G^{Q ) for 4 different values of the 
constituent mass M (Gorski et al.l992). 

5.4. The electric charge distribution of the proton for a constituent mass of 
M = 465MeV (full line) splitted in valence (short dashed) and sea (long- 
short dashed) (Wakamatsu 1991, Gorski et al.l992). 

5.5. The electric charge distribution of the neutron for a constituent mass 
of M = 465MeV (full line) splitted in valence (short dashed) and sea 
(long-short dashed) (Wakamatsu 1991, Gorski et al.l992). 

5.6. The proton magnetic form factor G^{Q ) for 4 different values of the 
constituent mass M (Christov et al. 1993a). 

5.7. The neutron magnetic form factor G^j{Q ) for 4 different values of the 
constituent mass M (Christov et al. 1993a). 

5.8. The magnetic moments for neutron and proton including contributions 
up to first order in the rotational frequency O are presented in dependence 
of the constituent quark mass M. The calculations are performed with 
rriTT = 139MeV. 
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5.9. The axial coupling constant (7^ including contributions up to first order in 
the rotational frequency O are presented in dependence of the constituent 
quark mass M. The calculations are performed with TOtt = 139MeV. 

6.1. The deviation of the theoretical mass from the experimental one is shown 
for the hyperon spectrum and for the Yabu-Ando treatment and a con- 
stituent quark mass of M = 419MeV (Blotz et al. 1993b) 

6.2. The deviation of the theoretical mass from the experimental one is shown 
for the E and the A, comparing the perturbative and the Yabu-Ando 
method for M = 419Mey (Blotz et al. 1993b). 

6.3. The hadronic part of the isospin splitting for the octet baryons from our 
theory (Praszalowics et al. 1993) compared with the experimental ranges 
for these splittings according to Gasser and Leutwyler (1982). 

7.1. The selfconsistent vector and axial vector fields u, p, Ag and Aj^ for a con- 
stituent quark mass M = 340MeV. The parameters are fixed according 
to the on-shell definition for the two-point functions, i.e. A = 877MeV, 
gp = 4.61 and g^u = 2.24. 



7.2. The polar field A(r) = T^-\/c^(r) +n'^{r) (which equals one in the non 
linear model) of the linear model is shown for different constituent quark 



masses. 



A.l. Contour integration for evaluating static observables of a system with 
baryon number B = (closed lines) and B = 1 (dashed lines) in Minkowski 
and Euclidean space. 
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TABLE CAPTIONS 

1.1. The current quark masses and corresponding charges from the QCD vac- 
uum (Gasser and Leutwyler 1982) 

2.1. The cutoff A, the n parameter, the original NJL coupHng G, the quark 
condensate uu, the current quark mass mo and the vacuum energy density 
for the non-covariant (0(3)), covariant (0(4)), single (PTl) and double 
(PT2) step proper time and the Pauli-Villars (PV) regularization, defined 
in Sect. 2 and Sect. 6. 



3.1. The critical values of the constituent quark mass Mcr and the corre- 

Ma, 



sponding values for the coupling constant gcr = —p- c^s well as the ratio 



Acr = mZ^ ^°^ "^TT ~ ^ (chiral limit) and ra-,^ = 139MeV. Solitonic 



^cr 



solutions of the non-linear model exist if M > Mq^ {g > gcr, A < Acr) 
(Meissner Th and Goeke 1991). 

5.1. The isoscalar n and the isovektor n magnetic moment for 4 values 
of the constituent quark mass M. The value in brackets denotes the 
contribution of the sea quarks. (Wakamatsu and Yoshiki 1990). The 
calculations are performed in the lowest non-vanishing order in the the 
cranking frequency, i.e. up to 0{fl^ ') for the isoscalar quantity and up 
to 0{Q,^ ') for the isovector magnetic moment. 

5.2. The pion nucleon coupling constant g-jj-NN ^^ shell q = m^ and at q = 
in the lowest order of the cranking frequency (O ) for 3 values of the 
constituent quark mass M. For the calculation the nucleon masses M^ 
from Fig. 4. 2 have been used. The experimental value of g^NNiQ = 0) 
is extracted from the Goldberger-Treiman relation gj^NNi^) = ~r^9A- 
(Meissner Th and Goeke 1991). 

6.1. The contribution of the valence and the sea part of moments of inertia 
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for M = 391.5MeV and M = 418.5MeV. 

6.2. Values of < Y,T,Ts \ DgS I YR^J.-h > and of < Y,T,T^ \ Dg^J^ \ 
Yr, J, - J3 >• 

6.3. The deviation of the theoretical mass from the experimental value is 
shown for the perturbative treatment with a constituent quark mass M = 
390.6MeV and the Yabu-Ando method with a constituent quark mass 
M = 418.5MeV. For two values of the strange current quark mass, these 
deviations are shown together with the absolute experimental mass. The 
theoretical value of the E*-mass is adjusted to the experimental one. 

6.4. The axial vector coupling constants (7^^, g"^ and g^ for M = 423.5Mey, 
given in the lowest order contribution (O ) and with linear rotational (O ) 
and additional strange quark mass (ml) corrections from the effective 
action, evaluated for wave functions, which contain linear mg corrections. 
These are compared with 'experimental' numbers from recent EMC and 
SMC measurements (EMC, Ashman et al. 1988,1989; SMC, Adeva et al 
1993) 

7.1. A comparison of the mean field energies ( in MeV ), the axial nucleon 
coupling constant (7^, the isoscalar mean squared radius < r > ( in fm 

) and the value of the polar field at the origin $(0) = y- y o"2(0) + 7r(0)^. 
For a constituent quark mass M = 340MeV the u coupling constant g^j 
is increased up to the correct value g^ = 2.24, which is in accordance 
with the mesonic sector. The linear model is used for the calculation. 

7.2. The same as Tab. 7.1 in the linear and non linear version of the Nambu- 
Jona-Lasinio-model with anpu and Ai mesons for different constituent 
quark masses M. 

8.1. Various observables (energy mean field energy Ej^jp, isoscalar electric 
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mean square radius {R ) jf^^ moment of inertia G, axial vector cou- 
pling gj\) for the self consistently solved NJL (valence part and sea part) 
(1. column), the selfconsistently solved NJL, where the sea part of the 
observables has been approximated by the 2nd order heat kernel or gra- 
dient expansion (valence part and mesonic approximation of the sea part) 
(2. column) and the selfconsistently solved CSM (valence part and mesonic 
part) (3. column). 

A.l. Notation for 4- vectors A^, the metric tensor g^^y, the scalar product of two 
4-vectors and the Dirac 7-matrices in Minkowski and Euclidean space. 
Hereby Greek indices run from — 3 in Minkowski space and from 1—4 
in Euclidean space, whereas Latin indices generally run from 1 — 3. 
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